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EXECUTIVE SUMMARY 


Objectives 


This contract dealt with the use c e the emissive or thermal infra- 
red band (10.5 to 12.5 pm) of the Heat Capacity Mapping Mission (HCIM) 
for specific agricultural applications: freeze damage assessment, 
planting date advisory, crop water stress, and temperature as affected 
by vegetative cover. More general or synoptic applications dealt with, 
for example, temperature gradients inland from the coast and the differ- 
ence between HCMU indicated and screen height air temperature maxima 
(THCMM-Tmax, daytime scenes) or minima (THCMM-Tmin, nighttime scenes) 
across the 220-km wide test site in extreme south Texas. In addition, 
the 0.55 to 1.1 pm band, which we refer to as the reflective band!/ 
was investigated as it helped edit clouds from the data, indicate major 
soil associations and as it related to vegetative cover. 

Scope of Work 

Computer compatible digital magnetic tapes (CCT) of 10 daytime 
emissive ( I R) and reflective (VIS) scenes and six nighttime emissive 
scenes were obtained, preprocessed and analyzed. The CCT were retro- 
spectively ordered after viewing hundreds of photo-products supplied 
to us on a standing order basis. The scenes obtained covered the time 
span 27 May 78 to 26 Feb 79. 

The test site we used was encompassed by longitudes 97°00' to 
99°30‘ W and latitudes 25°3Q' to 27 ° 00 ‘ il. The generalized soil map of 
the area was digitized, and the soil association number for each 20- 
seconds of latitude and longitude was entered into computer memory. 

For more detailed analyses of smaller ground sites, 44 sites represent- 
ing typical natural range plant communities, cities, water, and irrigated 
and dryland crop sites were chosen by reference to high altitude photog- 
raphy and tandsat images. Twenty-four weather station sites from which 
daily maximum and minimum air temperatures were reported during the 
winter months ware used to relate HCsfi indicated surface temperature 
to reported minimum temperature on cold nights. Ground truth was 
obtained weekly from 01 Jun to 15 Jul 78, for over 100 fields that 
constituted 22 intensive sites where soil water stress and ground cover 
were to be related to HC.UM indicated surface temperature. 

The "representative”, weather station, and intensive sites mentioned 
above constituted 94 ground locations in the test site area that were 
located in the HCUi! data by latitude and longitude of the center pixel. 
The temperature and albedo of the center pixel (1x1 matrix)— which re- 
presented 35 -ha on the ground—as well as the temperature and albedo 


— This band includes a portion of the visible (3.4 to 0.7 pm) and 
reflective infrared (0.75 to 1.35 pm) bands. Therefore, it is more 
accurate to refer to it as the reflective band than a visible (VIS) 
band, as is done in the HCMM User's Guide. 






of the 3x3, 5x5, 7x7 and 9x9 pixel matrices centered on it were deter- 
mined. Standard deviations of temperatures and albedos and the dominant 
and the next most dominant soil association for the 3x3 and larger 
matrices 'were determined. On 15 Jun and 12 Jul 78, Landsat data 
were used to calculate vegetation indexes for most of the 94 sites. 

The vegetation greenness indexes were used in conjunction with ground 
truth, weather data, calculations of root zone soil water depletion 
and actual to potential evapotranspiration ratios from the USDA-ARS 
Irrigation Scheduling Program to relate HCMM radiometric temperatures 
to plant cover and crop water stress. 

The above investigations were supported by three other facilities or 
efforts. One was a central weather station where numerous weather 
parameters were continuously monitored and their hourly integrated 
values were logged. Another was two 6-m square fallow plots, one of 
which was irrigated weekly and the other of which received only rainfall. 
In these plots, soil temperature was monitored at the 0.5, 5, 10, 20, 50, 
and 100-cm soil depths and recorded hourly on a strip chart recorder. 
These data helped to investigate the thermal response of the shallow 
soil layers to meteorological variables and the difference in thermal 
response of wet and dry soil. The third effort was daily maximum and min 
imum soil temperature at the 10-cm depth at four monitoring sites, during 
the period 01 Jan through 31 Mar 79, which the NWS-N0AA used for planting 
advisories. 

The above efforts provided the data for the contract objectives. 

In practice, however, we found that the HCMM scenes of our test site 
typically contained some clouds and that subvisible cirrus layers that 
affected the thermal data were also present over a part of some 
scenes. Consequently, considerable effort was expended in devising 
techniques for editing cloud-contaminated pixels from the data (see 
report section 2.2.1). Vie settled on what we called "modified cluster • 
screening". In this procedure, scatter diagrams of the reflective and 
emissive pixel values were displayed by the printer-plotter. Clouds had 
random combinations of values in the two bands whereas land features 
and water had values that clustered. The analyst examined the scatter 
diagrams and subjectively positioned the line segments that separate 
cloud from noncloud pixels. The line segments were described mathemat- 
ically and used to machine edit the data. It was found for eight summer 
overpasses that, on the average, each 10'i cloud contamination reduced 
the radiometric temperature of land features by 0.9 C. Because temper- 
ature differences between surface features are small, cloud-contaminated 
data ware edited out before most analyses were run. However, exemplary 
analyses are included that used data before and after Cloud-contaminated 
pixels had been removed to illustrate the degradation in relationships 
observed. 

Two other considerations in interpreting radiometric temperatures 
from HCMM were encountered. One was the low radiometric surface tem- 
perature determined (below ambient air temperature) under humid atmos- 
pheric conditions unless the radiative transfer model provided by 
Or. John Price, HCMM Project Scientist, was used to make the "surface 
to satellite" temperature correction based on radiosonde observations 


made by NWS-NOAA within the experiment test area (Brownsville). The 
correction for dry winter days (1 ci<< preci pi table water in the atmo- 
sphere) was nil whereas it was up to 9 C for humid summer days with up 
to 4.4 cm precioitable water in the atmosphere (see Table 1.3). In 
this report, anal'/ses involving temoerature differences between dates 
or between the air and the surface usually have incoroorated the radi- 
ative transfer model temperature correction whereas within date or 
among date regression analyses have not because the correction would be 
the same for all ground sites in the test area based on a single radio- 
sonde observation. 

The other consideration was number of HC M M pixels used to represent 
a given ground test site (field or larger area). In multiple regression 
analyses the multi Die correlation coefficient, R, increased slightly 
as number of HCMM pixels in the sample increased from one (central pixel) 
to 81 (matrix of 9x9 pixels centered on the central one), but the im- 
provement was not statistically significant at the p = 0.05 level (see 
report section 2. 3. 1.1). The standard deviation of the mean decreased 
but the value of the mean was not much affected. This was important 
in our study because the IFOV of HOT1 is 36 -ha— larger than all but a 
few percent of the crop fields. 

Results and Significance 


Synoptic analyses . A test site 220 -km E to W and 80-km N to S at 
the southern tip of Texas was used for several synoptic thermal analyses. 
Isotherm lines at 2 C intervals adequately displayed the thermal patterns 
within the test site. The HC‘ ? M radiometric temperatures by soil associa- 
tions from 10 daytime overpasses were typically highly correlated. 

Summer nighttime temperatures tended to be uniform across the test area 
and close to the dewnoint temperature which is reached each night. On 
winter nights with strong radiational cooling and dewpoints lower than 
the air reaches, there is a temperature gradient across the test site 
that is dominated by distance from the coast. Warm land areas adjacent 
to inland water bodies were observed in both summer and winter niahttime 
thermal data. 

An "aridity index" defined as the difference between the daytime 
HCMM radiometric temperature and maximum air temperature for the day at 
the site of weather stations (THC.MM - T“’ax) increased with distance from 
the coast, except for a short time following general rains. The differ- 
ence between the HCMM nighttime radiometric temperature and minimum air 
temperature (THC.’IM - Tmi n ) showed that THC M M was within 2 C of the dew- 
point temperature on nights when the air reaches the dewpoint 
temperature. 

The HCMM data contain useful direct synoptic information such as 
landscape temperatures durinci freezes, tenoerature gradients across 
scene areas and with respect to water bodies and estimates of minimum 
or dewpoint temperature. Day minus night radiometric temperatures, or 
radiometric and air temperature differences also appear to be useful 
aridity indexes. For the latter, the radiometric temperatures for the 
land surface need to be corrected for atmospheric path radiance. 


x 








Freeze hazard . Areas of least freeze hazard in our south Texas 
test area were identified by computer processing of HCMM data from cold 
nights. First a temperature frequency distribution was made of the pix- 
els from 21 Oct 73 that comprised our test area, and identification was 
made of all pixel locations whose temperature was in the top third* The 
same procedure was used with the cold night data of 26 Feb 79 (Had us- 
able data from other cold nights been available, they would have been 
processed to determine the locations of warmest pixels on those nights 
also.). Then the locations of warm pixels from both dates were super- 
imposed to form a map showing the areas that had least freeze hazard. 

A bold symbol via s used to identify locations that were warm on both 
dates. This approach offers promise as a useful method of classifying 
land areas according to freeze hazard. Because of variations in temper- 
ature patterns from date to date, information from several freeze or near 
freeze events should be combined in determining the safest areas to grow 
temperature-sensitive crops. 

HCMM nighttime data from the serious freeze of 03 Jan 79 could not 
be validly analyzed for landscape temperatures because an absorbing-- 
emitting atmospheric layer, probably subvisible cirrus (SCi), obscured 
surface emissions from view of the satellite. The SCi layer, which was 
not apparent to ground observers, retarded the loss of nocturnal terres- 
trial radiation to space and slowed the cooling of the landscape and air 
that night. The ability to detect an atmospheric layer that blocks the 
loss of terrestrial radiation should have important application in real- 
time freeze forecasting. Knowledge of. the presence of such a blocking 
layer would be justification for raising upward the minimum temperature 
forecast for a cold night and could save growers from undertaking expen- 
sive freeze protection measures unnecessarily. 

On the coldest SCi -free night for which we have HCMM data (26 Feb 
79) the HCMM indicated surface temperatures surrounding weather stations 
approximated the minimum air temperatures observed at the stations. 

Since the HCMM overpass occurred at 0209 CST and the minimum air temper- 
atures occurred about dawn, the HCMM radiometric temperatures were below 
air temperature at overpass time by an amount that approximately equaled 
the air temperature decline between overpass time and dawn. The results 
suggest that the satellite radiometric temperatures can be used to esti- 
mate minimum nighttime air temperature synoptically or at particular 
locations where air temperature is needed but ground observations are 
unavailable. 

Planting advisory . HCMM radiometric temperatures on a clear night 
early in the 1979 planting season (26 Feb), when there was little atmos- 
pheric attenuation, were very close to mi iimum air temperature, 3.1 C. 

The radiometric temperature of the sandy upland was 1.4 C and differed 
from that of the loams and clays (3.3 C) that had the same centroid dis- 
tance from the coast. Three environmental variables— average daily air 
temperature, insolation at ground level and extraterrestrial radiation- 
accounted for from 89 to 93;' of the variation in maximum daily, minimum 
daily and average dailv soil temperature at the 5-cn depth in dry plots 
and from 88 to 91% of the variation for wet plots. Standard errors of 
the estimate of maximum, minimum, and average daily soil temperatures 
from the three environmental variables ranged from 1.5 to 1.9 C. 


Estimation of planting depth soil temperatures from radiometric 
temperature! appears feasible, provided the relation between planting 
depth and surface temperature is ynderstood. The estimates should be 
superior to esti rates currently made from air temperature observations. 

Crop water Stress . HCMM temperatures for sorghum, sugarcane and 
citrus crop sites were related to soil water depletion (DEPL) since 
the last rain or irrigation calculated from the USDA-ARS Irrigation 
Scneduling Program (i'SP), to days since irrigation or rainfall (DSR) 
and to the ratio of estimated actual to potential evapotranspiration 
(ET/ETp) the day of and the dsy prior to HCMM overpasses. For a com- 
bined crops data set involving three HCMM overpasses, radiometric 
temperatures (means of 3x3 pixel arrays, T3) were highly significant- 
ly negatively correlated with ET/ETp (r - -0.704**) and significantly 
positively correlated SfTfch DSR (r = 0,484*). Temperatures were also 
highly significantly correlated with insolation on the day of the pass 
(I OOP) r = 0.812**; the vegetation greenness index (PVI6) calculated 
from Landsat for the corresponding dates (r = -G.642**) and maximum 
air temperature (Tmax) tha .day of the HCMM pass (r = 0,619**). The 
linear correlation coefficients between HCMM radiometric surface 
temperature and variables important in evapotranspiration for seven 
dates (27 May through 08 Jul 78) consisting of 51 observations were: 
insolation the hour of the HCMM overpass (IHROP), r = 0.668**; IOOP, 
r - 0.660**; ET/ETp, r = -0.581**; DSR, r = -0.477**; Tmax, r = 0.378** 
and percent vegetation cover (PC), r -0..277**. For this data set, 

Tmax was highly significantly correlated with DSR, DEPL and ET/ETp. 

The data demonstrate the strong correlation between variables 
related to evapotranspiration and HCMM surface temperatures, as well 
as the strong coupling of maximum daily air temperature to availability 
of water for evapotranspiration. In general, the findings indicate that 
surface temperatures observed by satellites, such as HCMM and the NOAA- 
§ to -9 series, are responding to shifts in the energy balance of the 
Earth's surface that are strongly affected by vegetative cover and 
availability of water for evapotranspiration. Such interpretations 
could have applications under the interagency AgRISTARS program as they 
relate to crop condition and yield assessments. 

Plant cover. HCMM temperatures for five crops (sorghum, citrus, 
sugarcane, cotton and grass pastures) as one data set were significantly 
correlated with the vegetation indexes GVI, PVI apd PVI6 calculated 
from Landsat multi spectral scanner data acquired wi this four days of the 
HCMM overpasses but not with percent cover (PC) or plant height (PH). 
Weather (insolation, saturation deficit of the air and daily maximum air 
temperature on HCMM overpass dates), geographic (distance of sample sites 
from the coast) and time (days since irrigation or rainfall) variables 
influenced HCMM temperatures as much or more than vegetative conditions 
when crops that differed in stage of maturity, planting configuration 
and PC vs. PH relation were combined in one data set. For cotton and 
sugarcane which were increasing in ground cover during the observation 
period (27 May to 15 Aug 78), the cooling effect of increasing ground 
cover was about equally offset by the seasonal increase in ambient air 
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and agricultural landscape temperature. In general, the spread in HCMM 
temperatures among crops was narrow and the unconnected temperatures 
were below ambient air temperatures at the time of the HCMM data acqui- 
sitions, 1400 to 1416 CST, 

From these studies we anticipate that surface temperatures observed 
by the polar orbiting operational meteorological satellites (flOAA-6 to 
-9 series) will correlate with vegetation indexes calculated from their 
visible and near-infrared bands and be interpretable in the AgRISTARS 
program for synoptic seasonal crop development and drought occurrence 
applications. 

Recommendations 


Radiometric temperatures require interpretations along energy 
balance lines. Insolation (solar irradiance) needs to be known and 
available for the interpretations as does soil water availability for 
evapotranspi ration. Therefore NASA should strongly support the efforts 
of NOAA to make insolation data available operationally from the GOES. 

Soil water balance models for the near surface soil layer and root zone 
are needed that can be readily implemented over a wide range of environ- 
mental conditions. It appears also that microwave data should be made 
more widely available for interpreting thermal data, and vice-versa . 

Our experience in these studies is that radiometric temperatures 
from satellite altitudes need to be routinely adjusted for radiative 
transfer losses through the atmosphere when preci pi table water in the 
atmosphere exceeds approximately 1 -cm. The radiosonde data for such 
corrections are available only for a limited number of sites in the 
U. S\ and those available were not taken at the time of HCMM overpasses. 
Consequently, provision of alternative sources of the data for such 
corrections is in order. For example, atmospheric sounders already or- 
biting may be able to supply the needed parameters. Or it may be that 
atmospheric sounders need to be placed on future earth observation 
satellites alongside the multi spectral or pushbroom scanners. 

Agriculturalists tend to associate specific observations with 
individual crop fields. However, IFOV's the size of those of HCMM (36- 
na) or the N0AA-6 to -9 series of satellites (l-km2) result in smoothing 
of the data across field boundaries so that it will be necessary to think 
in terms of larger sample sites on the ground that are representative of 
particular physiographic areas or ecological environments and less about 
individual cropped or fallow fields. 

In these studies the vegetation indexes calculated from visible and 
reflective infrared bands of Landsat were significantly correlated with 
the thermal band data from HCMM. That is, the vegetation indexes provided 
useful vegetation cover and condition information for intepreting the 
thermal data. Therefore, in spite of the disappointing experience with 
the thermal band on Landsat- 3, wa recommend that NASA persevere in in- 
cluding thermal sensors on future earth observation systems. 
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1.0 INTRODUCTION 


1 .1 Ob jecti ves 

The objectives of our HCMM investigation took advantage of the 
optimum day and night timing of the HCMM passes. Agriculturally 
related studies were conducted in the subtropical Lower Rio Grande Valley 
of south Texas. Rainfall decreases by l-cm/4.5-km inland from the coast 
at the test site giving a natural gradient in plant species and soil 
water available for evapotranspiration. An important irrigated area of 
260,000 ha provided fields that vary in irrigation date (degree of water 
stress), soil type (crop rooting depth), and crop species. In addition, 
freezes are a hazard to high value/ha citrus, sugarcane, and vegetables. 
Our research was organized for seasonal emphasis: 


1. Synoptic studies (year around) 

a. Document surface temperature gradients inland from coast 
and around inland lakes 


b. Use HCMM to develop synoptic aridity index, THCMM - Tmax 

c. Relate nighttime HCMM temperatures to ambient air dew- 
point temperatures 

d. Relate HCMM surface temperatures to ground observed day 
and night air temperatures 

2. Freeze hazard (Dec to Feb) 

a. Identify areas of least freeze hazard for growing tempera- 
ture-sensitive crops 

3. Planting date advisory (Jan to Apr) 

a. Relate nighttime soil temperatures to major soil texture 
variations over the Lower Rio Grande Valley 

b. Correlate HCMM soil surface temperatures with ground 
temperatures reported by NOAA in agricultural advisories, 
or with those directly measured at 5 -cm planting depth 

4. Plant cover (year around but primarily Mar to Jun) 

a. Document effect of percent ground cover on HCMM measured 
surface temperatures 

b. Identify bare fields 

5. Crop water stress (Apr to Sep) 


Relate HCMM temperatures to soil water depletion and days 
since irrigation or rainfall 
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1.2 Data Sources 


1.2.1 HCMM images and computer compatible tapes (CCT) 

The 700 -km wide swath of HCMM covered our test site (97°0G' to 
99°30*W longitude and 25°30' to 27°00'H latitude) partially or totally 
in the daytime on reference days 3, 8, and 14 whose NNW-bound orbits 
centered about over Brownsville, Laredo, and Houston, Texas, respec- 
tively. Nighttime overpasses, which occurred on the SSW-bound orbits 
of reference days 6, 11, and 1, had subsatellite cities, for relative 
positioning of the orbits, of McAllen, Houston, and Laredo (Falcon 
Lake), Texas. The dates of overpasses for the above reference days 
and ground identifying cities are summarized in Table 1.1. 

Because we had a standing order with the Information Processing 
Facility (IPF), NASA, Goddard we received single 9x9-inch (230x230-mm) 
black and white prints of each coverage of our site as relayed to the 
Merritt Island, Florida, receiving station during the period May 78 
through May 79. (Our site was on the very fringe of the range for 
both the Greenbelt, Maryland, and Goldstone, California stations; con- 
sequently they received little data for our site). 

We received hundreds of prints to screen for data quality and cloud 
cover. From these, we selected 10 daytime emissive and reflective band 
scenes, and six nighttime emissive scenes for which we ordered 9 track 
800 bpi CCT, single negative transparencies, and single positive prints 
for use in analyses and reporting. The dates of the overpasses, their 
ID, type of scenes (day or night), and average cloud cover over the 
test site are given in Table 1.2. For analyses purposes, the scenes are 
grouped chronologically as indicated in report section 1.1. 

Because of our low latitude, consecutive day /night or night/day 
(12 hour difference) overpasses were not possible. However, day/night 
passes 36 and 60 hours apart were possible. We. could retrospectively 
order day/night temperature difference (AT) and thermal inertia differ- 
ence (ATI) products from Goddard. We did retrospect!' . ..ly order the AT 
and ATI processing of the 08 Jul 78 day and 06 Jul 78 night scenes. 

\ 

1.2.2 Intensive test sites 


For crop water stress and plant cover studies, 18 intensive test 
sites were selected using aircraft support photography, Landsat imagery, 
and visits to field sites that contained contiguous fields of the same 
crop or very large single fields --defined as land parcels managed alike 
(same crop species, planting date, cultivations, etc.) that ware not 
dissected by field access roads, public roads, canals, fences, etc. 

The 18 segments were composed of 99 fields that were ground truthed. 

In four segments, large fields of two different crops occurred so that 
the data set consisted of 22 fields. The observed sites ranged in area 
from 1.03 km2 (108 ha) to 12.40 km2 (1240 ha). Representative fields 
ranged from 18 ha to 1240 ha and averaged 248 ha in size. These sizes 
were the basis of examining the HCMM temperature data for the center 
pixel of each site and the 3x3 pixel matrix (325 ha area) centered on 
the center pixel. 
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We waited until HCMM was launched and operational before we began 
weekly ground truthing individual fields (no ground truth data was 
obtained when muddy roads prevented field access) from 30 May to 13 Jul 
78. The information obtained included crop species, plant height (PH), 
percent crop cover (PC), percent weed cover, stage of crop maturity, 
and date of irrigation. Notes were taken about plant condition (abnor- 
mal defoliation, leaf curling or rolling, color), soil condition 
(crusted, wet or dry soil surface, tillage in progress or recently com- 
pleted, standing water), and obvious insect or disease incidence. 

The data were summarized for each field and filed on discs for 
merging with data from other sources. The crop, longitude and latitude, 
and site number for each of these sites is given in Table 1.4 and their 
distribution in the test area is shown on the site map of Figure 1.1. 

The sites were concentrated in Hidalgo County to minimize the time and 
mileage to visit them. 

1.2.3 Representative sites 

Forty-eight ground locations were chosen within the test site that 
represented typical water, urban areas, irrigated and dryland cropped 
sites, native and improved rangeland, wildlife refuges, and special 
interest sites such as barren Padre Island (sand) and a native oak 
grove on deep sand. These sites were chosen for internal homogeneity 
by using high altitude photography, Landsat images, and various maps of 
the area. They ranged in latitude from 40-km off the coast to 160-km 
inland. They were intended for synoptic analyses and to be large enough 
so that HCMM data in up to 9x9 pixel matrices could be examined for 
these sites. These sites are listed by number, longitude and latitude, 
and a brief description of each is given in Table 1,5. A map showing 
their locations and distributions is Figure 1.2. 

1.2.4 Supplementary Landsat data 

The earth observation satellite, Landsat series, with its 4-band 
multi spectral scanner (MSS) passed over the test site on two dates of 
low cloudiness, 15 Jun (Landsat-2) and 12 Jul 78 (Landsat-3) , when the 
fields were being ground truthed. The CCT for these two dates were 
ordered from the EROS Data Center, Sioux Falls, South Dakota, and the 
digital count data for each Landsat band were extracted for subsamples 
of the 22 intensive sites as well as for many of the "representative" 
sites. 

The digital samples were used to calculate Kauth transformations, 
greenness (GVI) and brightness (BRT) (Kauth and Thomas, 1976), and 
perpendicular vegetation indices (PVI, PVI 6) (Richardson and Wieaand, 
1977). 

These vegetation indices were needed to corroborate ground truth, 
and to provide spectral indicators of vegetation greenness and ground 
cover for additional ground sites that were not ground truthed. These 
data were studied in relation to HCMM indicated earth surface and crop 
canopy temperatures. 
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1.2.5 Climate and weather 


Supporting weather data were collected at a centrally located 
environmental station (Weslaco, Texas) where continuous measurements 
were made of solar radiation, total hemispherical radiation, wind 
movement, air temperature, and vapor pressure (dew point). These 
data were integrated for each hour of the day and archived as hourly 
means and daily totals. 

Daily precipitation and daily maximum and minimum temperatures 
for 15 reporting stations within the HCMM test site area that are 
regularly reported and published (NOAA-EDIS, 1978 and 1979) were 
summarized. Also, specifically for the HCMM study, a recording thermo- 
graph station was installed in an unrepresented rangeland area. When 
needed for analyses, data from the reporting station nearest our indi- 
vidual test sites were used to estimate their precipitation and air 
temperature. 

From Nov to Mar each year, flOAA-NWS uses volunteer observers to 
measure and report minimum daily temperature. Frequent observations 
are also made during periods of freezing or near freezing temperatures. 
These data were given to us by the NWS office at Brownsville, Texas. 

The data supplied had detailed air temperatures for interpreting results 
for 02 and 03 Jan 79, when a freeze occurred. The site number, coor- 
dinates, and names of such stations are given in Table 1.6, and their 
distribution over the test area is shown in Figure 1.3. 

Radiosonde atmospheric soundings of temperature, pressure and 
humidity (dew point depression) made by NOAA-f'WS at Brownsville on HCMM 
overpass days were used to correct the thermal satellite data by applying 
a radiative transfer model. Table 1.3 lists information from atmospheric 
soundings that were made nearest in time to the satellite data that we 
utilized. 

1.2.6 Aircraft support flights 

The contract proposal was submitted to NASA in Aug 75. At that 
time it was anticipated that HCMM would be launched in late 1976. We 
specified aircraft support flights on cold (the criterion used was that 
the predicted nighttime minimum be 5 C) clear calm nights in Jan and Feb 
77, and in May, Jun, and Jul 77. Five flight lines were used for daytime 
flights that totaled 310 nautical miles (574-km) . 

Flights were requested at low (3,000-m AGL) and medium (9,1 00 -ni) 
altitudes with the C-130 aircraft that was equipped with a portable radi- 
ation thermometer (2°F0V), MSS, and aerial cameras. The low altitude 
data were needed to examine differences and gradients within the fields 
and the medium altitude data were needed to make comparisons among fields 
High altitude (12,000 to 18,000-m AGL) data obtained with the 8,857 
airplane were needed for synoptic studies. 


-to: 






The Cl 30 aircraft overflew three flight lines at night on 20 Jan 
77 and part of the flight lines at 1,825-m and 7,900-ni on 26 Oul 77. 

The photography was valuable to help select the representative and 
intensive test sites. No overflights occurred during the period when 
HCMM data were being supplied. 

A locally-leased aircraft, equipped with a Zeiss RI!K 9-inch 
camera, overflew the intensive test sites on 09 Jun 78. The photographs 
were used to locate the test sites and individual fields in the line 
printer gray maps of the Landsat data, determine the size of individual 
fields by planimetry, verify land uses within the sites, document vari- 
ations in crop stands and vigor within individual fields, and document 
soil and topographic variations within test sites, 

1.2.7 Ancillary information 

Test fields and sites were located on U5GS 7 1/2 minute quadrangle 
maps using aerial photographs and features displayed on the maps as aids. 
Once located on the topographic maps, the center coordinate of each site 
was assigned a latitude and longitude to the nearest ±.001 degree. This 
became the center pixel that was searched for in the HCMM data, and we 
later refer to it as the lxl pixel matrix for the site (one pixel in 
latitudinal and 1 pixel in longitudinal extent). 

We carefully located the center pixel because each HCMM pixel repre- 
sented about 36-ha on the ground which corresponds to a displacement of 
adjacent pixels by 20 seconds of latitude and longitude, or .006°. That 
is, we did not want site location specification to be less precise than 
the retrievability of the HCMM data from the CCT. (The HCMM data were 
registered to the ground scene to within about ±1 pixel --see report sec- 
tion 2.1.1 for HCMM data preprocessing procedures). 

Soil association maps of the Lower Rio Grande Valley were provided 
in response to the HCMM study by the Soil Conservation Service (SCS) of 
the USDA. This information was digitized each 20 seconds of displacement 
to correspond to the HCMM pixel resolution. The identifying soil asso- 
ciation number at every 20 seconds of the latitude and longitude grid 
v/as stored on disc for use in interpreting the HCMM data. 

The SCS county soil surveys for Cameron and Starr Counties (Thompson 
et al. , 1972; Williams et al . , 1977) gave specific soil identifications 
for the intensive test sites in these counties. The soil data for test 
sites in Hidalgo and Mill acy Counties were either obtained from unpublished 
work sheets in the Raymondville and Edinburg SCS field offices or from 
the generalized map. 

1.2.8 Representative soil temperature, and experimental wet 
and dry plots 

Daily maximum and minimum soil temperatures measured at the 10-cm 
depth for four locations within the test area were obtained from the N0AA- 
NWS office in Brownsville for the period 01 Jan 79 to 31 Mar 79. These 
data are useful indicators of soil temperature for planting advisories 
and for deriving relations among environmental variables and soil temper- 
atures at the 19 -cm depth. 

5 





To provide additional data for direct comparison with HCMM- 
indicated temperatures and for predicting soil temperature at various 
depths from meteorological parameters, two 6.1 m square continuously 
fallow plots ware established. One plot received only rainfall while 
the other plot was also irrigated weekly. In each plot, thermocouples 
were installed at the 0.5, 5, 10, 20, and 100-cm soil depths, and their 
outputs were printed on a strip-chart recorder at hourly intervals. 

Data were recorded for the period 10 Aug to 29 Sep 73 and from 22 Nov 
78 to the present. The data provided information for relating HCMM 
data to existing ground and environmental conditions at the time of 
its acquisition. 
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Figure 1.3 Distribution of 24 weather stations for air temperature reference 
within the HCMM investigation test site. 
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Table 1.1 


Dates of HCMM overpasses by reference day orbits from launch 
through June 1979 that partially or totally covered our test 
site. A subsatellite city for relative positioning of the or- 
bits is given in parentheses for each reference day orbit. 



Davtime 



Nighttime 


Ref. 

Day 3 

(Brownsville) 

Ref. 

Day 8 
(Houston) 

Ref. 

Day 14 
( Laredo ) 

Ref. 

Day 1 
( Laredo) 

Ref. 

Day 6 
(McAllen) 

Ref. 

Day 11 
(Houston) 

1978 5/16 

5/21 

5/27 

5/14 

5/19 

5/24 

6/01 

6/06 

6/12 

5/30 

6/04 

6/09 

6/17 

6/22 

6/28 

6/15 

6/20 

6/25 

7/03 

7/08 

7/14 

7/01 

7/06 

7/11 

7/19 

7/24 

7/30 

7/17 

7/22 

7/27 

8/04 

8/09 

8/15 

8/02 

8/07 

8/12 

8/20 

8/25 

8/31 

8/18 

8/23 

8/28 

9/05 

9/10 

9/16 

9/03 

9/08 

9/13 

9/21 

9/26 

10/02 

9/19 

9/24 

9/29 

10/07 

10/12 

10/18 

10/05 

10/10 

10/15 

10/23 

10/28 

11/03 

10/21 

10/26 

10/31 

11/08 

11/13 

11/19 

11/06 

11/11 

11/16 

11/24 

11/29 

12/15 

11/22 

11/27 

12/02 

12/10 

12/15 

12/21 

12/08 

12/13 

12/18 

12/26 

12/31 


12/24 

12/29 


1979 


1/06 



1/03 

1/11 

1/16 

1/22 

1/09 

1/14 

1/19 

1/27 

2/01 

2/07 

1/25 

1/30 

2/04 

2/12 

2/17 

2/23 

2/10 

2/15 

2/20 

2/20 

3/05 

3/11 

2/26 

3/03 

3/08 

3/16 

3/21 

3/27 

3/14 

3/19 

3/24 

4/01 

4/06 

4/12 

3/30 

4/04 

4/09 

4/17 

4/22 

4/28 

4/15 

4/20 

4/25 

5/03 

5/08 

5/14 

5/01 

5/06 

5/11 

5/19 

5/24 

5/30 

5/17 

5/22 

5/27 

6/04 

6/09 

6/15 

6/02 

6/07 

6/12 

6/20 

6/25 

7/01 

6/18 

6/23 

6/28 

7/06 



7/04 

7/09 



Note: 


Orbit paths over Brownsville and McAllen came closest to being centered 
over our HCMM test site; however, passes over Houston and Laredo encom- 
passed it because of the wide swath. Ref. Day 3 (daytime) and Ref, Day 
1 (nighttime), and Ref. Day 8 (daytime) and Ref. Day 6 (nighttime) cov- 
erages, respectively, were 36-hours apart. Day and night coverages may 
occur each five (or six) days. The same orbit is repeated each 16 days. 




Table 1.2 Overpass dates, scene identification numbers, day (D) or night (N) 
scene, and percentage of clouds over our test site for HCMM scenes 
investigated. 


Overpass 

date 

Julian 

date 

Scene 

type 

Scene 

ID 

Percent 

clouds 

27 May 1978 

147 

D-VIS 1 

A-AOO 31-20160-1 

12.4 



IR 2 

A-A00 31-20160-2 


6 Jun. 1978 

157 

D-VIS 

A-A00 6 8-20060-1 

19.6 



IR 

A-A006 8-20060-2 


12 Jun. 1978 

163 

D-VIS 

A- A00 47-20150-1 




IR 

A-AQO'47 -20150-2 

11.3 

17 Jun. 1978 

168 

D-VIS 

A-A00 5 2-200 80-1 




IR 

A-A00 5 2-200 80-2 

25.8 

23 Jun. 1978 

174 

D-VIS 

A- A005 8-20200-1 




IR 

A-A00 5 8-20200 -2 

25.2 

3 Jul. 1978 

184 

D-VIS 

A- A00 41-20040-1 

18.8 



IR 

A-A0041-20040-2 


8 Jul. 1978 

189 

D-VIS 

A-A0073-20000-1 

3.4 



IR 

A-A00 73-20000-2 


15 Aug. 1973 

227 

D-VIS 

A-A0111- 20 040-1 

3.6 



IR 

A- A0111- 20040-2 


8 Nov. 1978 

312 

D-VIS 

A-A0196-19440-1 




IR 

A-A0196-19440-2 

0.1 

7 Feb. 1979 

38 

D-VIS 

A-A0237-19420-1 

0.2 



IR 

A-A0287-19420-2 


30 Jun. 1978 

181 

N- IR 

A-A0065-08260-3 

.... 

23 Aug. 197 8 

235 

N-IR 

A-A0119-08270-3 

- 

10 Oct. 1978 

283 

N- IR 

A- A0 16 7 -08240- 3 

- 

21 Oct. 1978 

294 

N-IR 

A- A0 178 -082 70- 3 

- ' 

3 Jan. 19 7 9 

003 

N-IR 

A- A02 52-0 8050- 3 

— . 

26 Feb. 1979 

057 

N-IR 

A- AO 306-0 8090- 3 

— 


1 

0.55 to 1.1 ym, or reflective band 

2 

10,5 to 12.5 ym, or emissive band 
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Table 1.3 Lowest level temperature (TLL) , lowest level pressure (PLL) , precipitable water (PW) , atmospheri 
optical depth (OD) , and surface to satellite temperature correction (IC) as determined by a radi 
tive transfer model for the HCMM scene dates investigated using morning (6am) and evening (6pm) 
radiosonde data. 


TLL 

(®K) 

PLL 

(mb) 

P m 

(cm) 

OD 

(--) 

TC 

( # c) 

297.5 

1012 

4.12 

1.221 

' 

3.5 

299.8 

1008 

4.53 

1.414 

5.6 

299.2 

1014 

3.83 

1.103 

3.79 | 

297.0 

1015 

3.91 

1.059 

3.97 ! 

i 

295.9 

1014 

3.63 

0.891 , 

j 

2.14 1 
• 

i 

299.2 

1010 

4.46 

1 . 174 

5.87 | 

| 

299.2 

1015 

3.65 

0.955 

3.99 | 

299.2 

1014 

4.32 

1.180 

4.45 I 


282.5 

1016 0.96 

0.125 

295.9 

1016 4.31 

1 . 108 

295.9 

1014 3.63 

0.891 

293.1 

1015 2.93 

0.795 

290.3 

1015 2.26 

0.560 

278.1 

1025 0.81 

0.092 



6pm 

' tLL 
(*K) 

PLL 

(mb) 

PW | 
(cm) 

OD 

(~) 

TC 

(°C) 

304.2 

| 

1101 

3.60 

0.914 ! 

6.45 

304.2 

1005 

4.35 

1.217 

6.83 

305.3 

1016 

3.88 

1.065 

7.21 

305.9 

1014 

3.94 

0.985 

8.66 

306.4 

1014 

1 

4.45 

1.186 

i 

9.77 

306.4 

1009 j 

3.41 

! 0.691 

5.85 



307.0 1014 

290.9 1017 


286.4 1020 


292.5 1026 


7.73 

0.88 | 0.133 0.68 

66 0.92 



.080 0.76 





Tabl'> 1.4 Site number, coordinates, and crops at 22 intensive test 
sites used in the HCMM investigations. 


Site No. 

Longitude 

Latitude 

Crop 


- - - - Degrees - - - - 


10101 

97.870 

26.368 

Sorghum 

30901 

98.245 

26.460 

Sorghum 

31001 

98.358 

26.404 

Sorghum 

31011 

98.354 

26.382 

Sorghum 

41341 

98.500 

26.333 

Sorghum 

41381 

97.865 

26.419 

Sorghum 

41421 

98.245 

26.372 

Sorghum 

10194 

97.885 

26.308 

Sugarcane 

10196 

97.875 

26.306 

Sugarcane 

30951 

97.902 

26.425 

Sugarcane 

41422 

98.250 

26.373 

Sugarcane 

41665 

97.884 

26 . 087 

Sugarcane 

51782 

98.120 

26.080 

Sugarcane 

10102 

97.872 

26.058 

Cotton 

10381 

97.920 

26.096 

Cotton 

20611 

97.984 

26.352 

Cotton 

30955 

97.898 

26.416 

Cotton 

20681 

98.134 

26.276 

Citrus 

20754 

98.402 

26.244 

Citrus 

31166 

98.262 

26.230 

Citrus 

20611 

97.970 

26.351 

Bermudagrass 

20651 

98.428 

26.321 

Buf f elgrass 
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Table 1.5 Site number, coordinates, and a brief description of each, land use 
at 48 representative sites used in the HCMM investigations. 


Site No. 

Longitude 

Latitude 

Site Description 


- - - Degrees - - - - 


001 

99.301 

26.854 

Falcon Lake (Inflow) 

101 

99.226 

26.969 

Rangeland, Zapaia 

111 

99.199 

26.610 

Falcon Lake (Main Storage) 

112 

99.010 

26 . 513 

Native Brushland, Nonsaline 

121 

99.034 

26.475 

Buffelgrass , Sandy Loam 

211 

98.972 

26.744 

Grain Sorghum (Dryland) 

212 

98.968 

26.659 

Saline Rangeland 

213 

99.000 

26.569 

Buffelgrass, Some Brush Reinfestation 

214 

98.901 

26.667 

Rangeland, Las Escobas Ranch 

221 

98.769 

26.415 

Rangeland, Ricaby Oil Field 

311 

98.674 

26.741 

Rangeland, LaGloria 

312 

98.688 

26.516 

Rangeland, Rincon 

401 

98.279 

26.777 

Rangeland (Red Soil) 

421 

98.435 

26.453 

Grain Sorghum (High Reflectance Soil) 

422 

98.384 

26.443 

Grain Sorghum (Low Reflectance Soil) 

423 

98.340 

26.334 

Mixed Crops, Mainly Citrus 

424 

98.268 

26.368 

Monte Christo Citrus Area 

425 

93.386 

26.405 

Citrus, Topographic Drain 

501 

98.194 

26.760 

Rangeland (Natural) 

521 

98.250 

26.366 

Sugarcane 

522 

98.136 

26.261 

Edinburg Citrus Area 

531 

98.231 

26.200 

City of McAllen 

532 

98.147 

26.075 

Santa Ana Wildlife Refuge 

533 

98.134 

26.219 

Mixed Crops , San Juan 

611 

97.879 

26.705 

Rangeland (Soil Coloration) 

612 

97.880 

26.604 

Rangeland (Numerous Swales) 

613 

97.855 

26.555 

Rangeland 

621 

97.947 

26.413 

Delta Lake (Irrigation Reservoir) 

622 

97.860 

26.376 

Mixed Cotton and Sorghum (Raymondville) 

623 

97.809 

26.356 

Mixed Cotton and Sorghum (Willacy-Raymondville) 

624 

97.881 

26.302 

Sugarcane 

631 

97.676 

26.212 

Mid-Valley Mixed Crops 

632 

97.789 

26.151 

Adams Gardens Citrus Area 

633 

97.775 

26.239 

Mixed Crops , Primera 

701 

97.623 

26.767 

Oak Trees ( Extensive Grove ) 

721 

97.555 

26.322 

Mixed Cotton and Sorghum (Willacy-Ramondville) 

* 722 

97.671 

26 . 466 

Sorghum and Cotton, Reflective Soil 

731 

97.699 

26.169 

City of Harlingen 

732 

97.572 

26.108 

"Homogeneous" Area 

801 

97.368 

26.827 

Padre Island Sand 

821 

97 . 441 

26.721 

Laguna Madre ( Shallow water ) 

831 

97.383 

26.233 

Laguna Atascosa 

832 

97.418 

26.067 

Grassland (High Water Table) 

841 

97.498 

25.905 

City of Brownsville 

901 

97.217 

26.750 

Gulf of Mexico - 12 km . E. of Padre Island 

931 

97.247 

26.169 

Laguna Madre (Deep Water) 

932 

97.201 

26.006 

Mud Flat 

941 

97.037 

25.787 

Gulf of Mexico - 12 km . E. of Washington 
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Table 1.6 Site number, coordinates, and station names of 24 weather 
station sites used in the HCMM investigations. 


Site No. 

Longitude 

Latitude 

Station 
Site Name 


- - - Degrees - - - 


1 

99.284 

26.868 

Zapata 

2 

98.865 

26.382 

Rio Grande City 

3 

98.386 

26.486 

McCook 

4 

98.408 

26.331 

Texan Gardens 

5 

98.396 

26.220 

Bates Power 

6 

98.117 

26.566 

San Manuel 

7 

98.158 

26.298 

Edinburg 

8 

98.159 

26.090 

Schuster Farm 

9 

97.988 

26.223 

Research Farm 

10 

97.962 

26.157 

Weslaco 

11 

97.791 

26.487 

Raymondville 

12 

97.869 

26.265 

Santa Rosa 

13 

97.777 

26.158 

Adams Gardens 

16 

97.211 

26.075 

Port Isabel 

17 

97.425 

25.913 

Brownsville Airport 

18 

97.968 

26.378 

Monte Alto 

19 

97.892 

26.064 

Mercedes 

20 

97.681 

26.209 

Harlingen 

21 

98.218 

26.202 

McAllen 

22 

98.604 

26.521 

Rincon 

23 

99.131 

26.566 

Falcon Dam 

24 

98.686 

26.311 

Garcias vi lie 

25 

98.901 

26.667 

Las Escobas Ranch 

26 

97.430 

26.558 

Port Mansfield 


2.0 METHODS, PROBLEMS, AMD SOURCES OF VARIATION 

2.1 Jata preparation and organization 

2.1.1 HC'1’1 whole test site (^100, 000 pixels) 

Once HCMM scenes (see Table 1.2 for a listing) were received from 
NASA Goddard, they were processed for analyses as follows where the 
name of the Weslaco USOA computer program is given in parentheses: 

a. Locate test site in magnetic tape, by referring to 
transparencies and/or Dicomed display (DSPLY). 

b. Determine digital count range of the test site (DIGN5). 

c. Prepare IR gray map on line printer-invert daytime data to 
match nighttime record sequence (because of onposite 
direction of satellite travel) (GRYN5 , DAYN5). 

d. Find georeference sites on IR data and produce location 
equations that relate pixel line and column to latitude 
and longitude. Cut and try to obtain minimum error of 
fit by discarding poorly identified georeference sites 
as indicated by computer analysis (Him). 

dl . Print IR gray map with georeference sites to check regres- 
sion equation coefficients that relate pixel line and 
column numbers to latitude and longitude. If this check 
is good then processing proceeds on to NTWKT. If not 
regression equation will need to be reworked (LATN5), 

e. Assign the IR value of the nearest pixel to each coordinate 
cell (20-seconds latitude x 20-seconds longitude). Store 
these data on disc (NTWK1). 

f. Repeat steps a through e with VIS data. 

g. Graphically represent six parallel IR (and VIS) transects 
that cross the test site (PROFL). 

h. Select and apply "threshold" values to the IR data so as 

to exclude all cloud contaminated data from further analysis. 
Determine threshold value by plot of individual IR pixel 
values vs corresponding VIS pixel values, and by examination 
of IP. (and VIS) transect graphs (50PFL, RATDF, TEMP 3, TEif>4). 

i. Prepare IR gray mao of test site using only retained, non- 
con tamina ted data (DISC2). 

j. Print out mini -maps of thermal environment surrounding all 
sites (Intensive, Representative, and Weather Sta. sites). 
Identify which pixels in the 9x9 cell scene were cloud 
contaminated (TEMP 3, TEMP 4) . 
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k. Calculate mean temperature and standard deviation of the 9x9 
pixel area surrounding (or comprising) the site. Include 
only cloud-free pixels in the calculations (TE;'!P3, TEMP4). 

l. Prepare histogram of temperatures in the 9x9 pixel areas, 
excluding cloud contaminated pixels (FREQ5, FREQ1). 

m. Prepare summary table of surface temperatures and albedos 
at all sites for lxl, 3x3, 5x5, 7x7, and 9x9 pixel areas. 
Tabulate mean values and number of cloud-free pixels in 
each size area (TEMP 3, TE'4P4) . 

n. Compare surface temperatures with air temperatures measured 
at the weather stations (TEMP3, TEMP4), 

o. Determine the mean temperatures and standard deviations for 
the 34 soil associations of the test area (exclude cloud- 
contaminated data from analyses) (S0IL2, S0IL3). 

p. (Required to be done only once) Determine the centroid of 
the area of each soil association. Information was used in 
analyzing for effect of proximity to the coast (CEfITR). 

q. Determine histograms of agricultural area and rangeland area 
(FREQ 2) . 

r. Using the frequency distribution of the agricultural area 
(excluding cloudy pixels), determine the temperature limits 
for Cold, Medium, and Harm areas, such that one-third of the 
area falls in each classification (FREQ2). 

s. Prepare map of test site showing areas of Harm, Medium, and 
Cold classification (FP.EQ2). 

t. Overlay selected nighttime IP, maps of test site that display 
the occurrence of pixels in the Warm category — of step r, 
above--by a. dark print symbol, pixels in the Medium category 
by a dot, and pixels in the Cold category by a blank space 
(CLOW*!). 

u. Using selected pairs of day/night IR scenes prepare a disc 
file of temperature differences, cell by cell. Process these 
"difference" data according to the procedures of steps g, i, 

j , k, 1 , m, n, q, r, and s (DIFF) . 

These processing steps are extensive but are necessary to derive 
the data to meet the various objectives of the contract. 
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2.1.2 HCMM temperature, albedo, and dominant soil by lxl, 

3x3, 5x5, 7x7, and 9x9 pixel cells for 94 ground 
sites 

The data were produced by the above processing steps j, k, 1, and 
m, The central point for each ground site of interest was identified 
on USGS 7Ja-minute quadrangle maps, and assigned latitude and longitude 
positions to the nearest 0.001 degree. Once this pixel was identified 
in the ground registered digital data, surrounding pixels were assigned 
to build the pixel configuration up to 3x3, 5x5, 7x7, and 9x9 pixel 
cells. 

The dominant soils could be determined because we digitized the 
generalized soil map of the four LRGV counties supplied by the SCS- 
USDA at 20-second intervals (0.006°) of latitude and longitude. Thus 
each HCMM pixel could be associated with the soil at that pixel location. 

2.1.3 Ground truth for 22 intensive sites 

The intensive test sites are described in section 1.2.2. A total 
of 99 fields were ground truthed weekly. The information recorded 
included Julian day and year of observation; sample stratum, segment, 
field and division within the field; crop; percent crop cover; percent 
weed cover; stage of crop development or maturity; plant height; unusual 
plant condition (defoliation, discoloration, etc.); soil surface condi- 
tion (wet or dry, cloddy, freshly cultivated or tilled, etc.); visible 
apparent disease or insect infestation; irrigation or rainfall amount 
(when determined); and, date. 

Aerial photographs of some of the sites were obtained with a 
locally based aircraft on 09 Jun, while photographs of the others were 
available from previous investigations dating back to 1972. The photo- 
graphs were planimetered to determine the hectarage of each field. The 
09 Jun photography was valuable for helping to locate fields of interest 
in the gray map printouts of Landsat data (see section 1.2.3) from which 
vegetation indexes were derived. The photography also verified within 
field variation in crop appearance and differences in maturity of adja- 
cent fields. 

A coding system for the ground truth was devised and the data were 
filed on discs for ready merger with data from other sources. 

2.1.4 Weather data 

Weather data described in section 1.2.5 were obtained from the 
central environmental station, 15 reporting stations in the test site 
area (N0AA-EDIS, 1978 and 1979), and N0AA-NWS volunteer observer reports 
of minimum daily temperature (from Nov to Mar). Twenty-four of the 
ground sites chosen for analyses of HCMM data were locations of NWS 
weather stations or volunteer observers (see Table 1.5), so that there 
would be a close correspondence bet ween reported air temperature and 
precipitation, and temperature deduced from HCMM data. Air temperature 
for the remainder of the 94 ground sites was taken from the nearest 
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reporting station. The temperature, rainfall, and other weather data 
were assembled for the HCMM overpass dates and times of interest, and 
were placed on disc files for ready access (WATST). 

Late in the investigation period, radiosonde data from Brownsville, 
Texas, was acquired for selected dates from Asheville, North Carolina, 
for use with a radiative transfer model provided by the HCMM Project 
Scientist. The objective was to find out the magnitude and direction 
of the temperature corrections it calculated for the HCMM data. 

2.1.5 Landsat data (15 Jun and 12 Jul 73) 

Landsat with its 4-band MSS passed over the test site on two dates 
of low cloudiness, 15 Jun (Landsat-2) and 12 Jul (Landsat-3). The 15 
Jun Landsat data were paired with the HCMM measured temperatures and 
albedos for the HCMM's 12 and 17 Jun 78, overpasses and the 12 Jul 
Landsat data were paired with the 03 Jul HCMM data. 

The Landsat scenes acquired from the EROS Data Center, Sioux Falls, 
South Dakota, were rectified to the ground scenes. The segments of the 
scenes, identified by latitude and longitude, over the ground sites of 
interest were extracted. These scene segments were used to produce gray 
map products of MSS 5 and 7. Use of USGS 7^-minute quadrangle maps and 
aerial photographs enabled the recognition of individual fields in the 
gray map printouts. Since the ground resolution of Landsat data is 
6.45 ha, samples of from 40 to 120 Landsat pixels from each band were 
used to represent a given ground site. Usually samples were from a 
single field that was representati ve of the larger area sensed by HCMM 
(35 ha/pixel or 325 ha/3x3 matrix of HCMM pixels). 

The Landsat digital samples were used to calculate the Kauth trans- 
formation ( GVI ) and brightness (BRT) (Kauth and Thomas, 1976), and the 
perpendicular vegetation indexes (PVI, PV 16) (Richardson and Wiegand, 
1977) (See chapter 7, section 7.3 for the equations that define the 
veqetation indexes,). 

The Landsat data were identified by ground site, field within site, 
crop, and date and placed on disc file for use with HCMM data. 

2.2 Special problems 

2.2.1 Coping with cumulus cloud contamination 
2. 2. 1.1 Daytime scenes 


Two of the 10 daytime HCMM scenes that we studied were clear. The 
test area in the other scenes was covered by from 3 to 26% clouds. These 
clouds were mainly of the "popcorn" cumulus type. If they were not 
densely aggregated we analyzed the scenes. To avoid erroneous conclusions 
from analyzing cloud contaminated data, it was necessary to develop 
methods of machine screening for clouds. The descriptions that follow 
enumerate the methods that we tried in our attempt to remove from 
analyses all pixels that were cloud contaminated. The problem of coping 
with contamination by s invisible cirrus (SCi ) layers not apparent in the 
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VIS images is discussed in section 2.2.2 (In this respect, we distin- 
guish between true clouds (mostly cumulus) that are apparent in both 
the VIS and IR images, and SCi layers that are apparent only in the 
IR.). 


2. 2. 1.2 Reflective (VIS) screening 

Screening reflective (VIS) data above a threshold value removed 
many of the clouds but there were areas of uncertainty. If the thresh- 
old was set low enough to detect thin clouds and pixels partly 
contaminated by clouds, the more reflective land areas such as fallow 
and dune areas were excluded. Also, the reflective criterion did not 
remove cloud shadows. 


2. 2. 1.3 Emissive (IR) screening 

Screening solely on the basis of emission was not feasible because 
water and cloud temperatures tended to overlap. 

2. 2. 1.4 Reflective/emissive ratio screening 

The ratioing of reflective (VIS) to emissive (IR) pixel values was 
tried for screening daytime clouds. This required that the data from 
the two wavelengths be registered so they could be overlaid and ratioed 
pixel by pixel. This was done by initially geometrically rectifying 
all HCM‘1 data to map coordinates as discussed in section 2.1.1 of this 
chapter. 

Screening for clouds by the reflective to emissive ratio was 
superior to either reflective or emissive screening alone. The ratio 
is useful because clouds are the brightest (highest digital count (DC)) 
scene component in the reflective channel and clouds are usually colder 
(lower digital count) than land features in the emissive data. Thus, 
the ratios, VIS/IR or IR/VIS, contrasted clouds very well versus other 
scene features for a given overpass date. 

A weakness of the ratio method is that the division between clouds 
and terrestial features is not a fixed ratio but varies from scene to 
scene and with season. 

2.2. 1.5 Standard deviation screening 

Tarpley (1979) described a process for eliminating cloud contamin- 
ated observations from GOES VISSR thermal data using standard deviations 
His reasoning was that the standard deviations for cloud-contaminated 
data would be more variable than noncloud contaminated data. However, 
data for water body-land interfaces such as for Falcon Reservoir, the 
Gulf of Mexico, and the bay (Laguna Madre), also have high standard devi 
ations; thus Tarpley 's scheme would incorrectly classify these water 
body-land interface areas as clouds. Also, the subvisible cirrus contam 
inated area does not yield a high standard deviation; thus, Tarpley's 
scheme would incorrectly identify the SCi contaminated areas as clear. 


21 


- . ,S§|£iii>*iai8a - ■»**-- • 


-Six***,* - 




It is possible that the Great Plains area, where Tarpley's scheme was 
applied, had a dry enough atmosphere and a dearth of water bodies so 
that neither was a problem. However, there are many other places in the 
world where this may be a problem. 

2 . 2 . 1 . 6 Movi ng -wi ndow screeni ng 

A cloud detection scheme was applied to emissive scenes that uses 
the ratio of the center pixel value (Xc) to the average pixel value 
(Xw) of a moving nxn pixel size window (Figure 2.1). Variations in 
Xc due to clouds are treated as high frequency image detail that is 
partially blocked by the lew pass filter, Xw data. Ratio values less 
than Xc/Xw-0.92 .we found experimentally to be due to a partial cloud 
and the corresponding center pixel values were censored. The value 
0.92 appears to be set by sensor NEAT and the other noises. 

This approach appears to distinguish between cloud values and 
surface and water body values satisfactorily. It does not censor 
vapor contamination. Unfortunately, it causes many individual pixels 
to be discarded in the vicinity of cloud edges because it is subject 
to the relative number of cloud-free and cloud-contaminated pixels 
in the nxn pixel window used as a filter. This may not be ojecti on- 
able where generalized results are desired. 

The same Xc/Xw=0.92 ratio value was applied to all test-site 
data for cloud detection. Figure 2.2 shows the effect of the 0.92 
ratio value for partly cloudy 27 .lay and 12 dun 73 flCMTl overpass dates, 
and clear days 21 Oct 73 and 07 Feb 79 on east-west temperature tran- 
sects across the test site. Mote the small variability on the right 
side of all scenes, the Gulf of lexico. 

2. 2. 1.7 Cluster screening 

In another cloud screening method, we took advantage of the fact 
that clouds are individualistic in scatter diagrams of reflective and 
emissive pixel values. Clouds have few locations where pixel values 
of both wavebands are the same as for the other scene features. That 
is, in two-dimensional spectral space, clouds do not cluster. Con- 
versely, pixel values from the Earth's surface tend to have numerous 
locations where both pixel values are the same as other locations. 

Thus, clusters are formed. The scatter in cloud data and the cluster- 
ing of clear areas were used to separate data with clouds from data 
for noncloudy areas. The scatter diagrams of Figure 2.3 show the 
clustering tendency of the data from land and water areas and the 
disperse nature of cloud data. 

Gray mans were prepared after discarding all reflective-emissive 
combinations that did not recur at least 17 times in the test area of 
approximately 100,909 pixels. Combinations that recurred at least 17 
times are identified by $ in the figures. 

Close examination of the gray maps prepared by the clustering 
criterion showed that some cloud shadows remained in the scene and 
certain land and water areas had been discarded. 
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2. 2. 1.8 Modified cluster screening 


Of the above methods of screening for clouds, clustering was the 
best. Most of the HCMM data discussed in this report were processed 
with a modification of the cluster method. 

The modification was applied by examining a scatter diagram of 
reflective and emissive pixel values (such as the samples in Figure 2.3). 
The analyst draws a subjectively positioned line (not necessarily 
straight) on the scatter diagram so as to separate cloud pixels from 
noncloud pixels. The line is then described mathematically and applied 
in machine processing of the data. All pixels having reflective- 
emissive combinations that fall on the cloud side of the line are 
discarded. 

A gray map prepared from the retained data was examined for un- 
desirable data (such as cloud shadows) or inappropriate exclusion of 
valid features (such as highly reflective land areas). If the gray 
map was unsatisfactory, the line separating the clouds from valid data 
on the scatter diagram was adjusted and the procedure repeated. Table 
2.1 lists the screening equations used in processing HCMM data. 

2. 2. 1.9 Nighttime scenes 

Consecutive day-night HCMM overpasses for thermal inertia analysis 
are not possible at our test site latitude, although pairs of scenes 
separated by several days have been examined. Thus, our main interest 
in nighttime emissive data has been confined to relationships with 
dew points and minimum air temperatures, and with cold-night surface 
temperature distributions. For these purposes, only clear night scenes 
have been analyzed, except that the freeze night of 02-03 Jan 79 was 
studied because of its economic importance, even though a subvisible 
cirrus layer was present. 

2.2.2 Coping with subvisible cirrus contamination 

A phenomenon observed in some HCMM scenes was illustrated in the 
emissive (IR) and reflective (VIS) wavelengths images for 03 Jul and 
15 Aug 78 (Figure 2.4). Visual comparison shows that cumulus clouds had 
the same pattern and extent in both wavelengths. However, there were 
additional patterns in the emissive scenes that had no counterparts in 
the reflective images. The more pronounced patterns seen in the emis- 
sive images were often observed in the reflective images, but the edges 
were very diffuse and effects were too subtle to be readily detected in 
automated processing of the digital reflective data, lie attributed 
the patterns readily discernabie in the emissive images, but not in the 
reflective images, to a subvisible cirrus (SCi) layer. 

Portions of the test area were under SCi layers on six of the 
10 daytime HCMM scenes considered in this report. These atmospheric 
layers occurred in the summer time. The 08 Nov 78 and 07 Feb 79 
scenes seemed to be free of SCi layers. 
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Identification of surface areas that were obscured from the sat- 
ellite sensor by a SCi layer was possible by simultaneously viewing 
the HCMM emissive (IR) scene and a location map on a zoom transfer 
scope. The SCi layers were clearly distinguishable from typical cumu- 
lus clouds because of their uniform and usually extensive nature. These 
characteristics contrast with the "popcorn" type cumulus clouds common 
in our test area. 


2.3 Sources of variation 

2.3.1 Small test sites (weather stations, representative 
sites, and intensive fields) 


We anticipated that many variables would affect the apparent Earth 
surface (agricultural landscape) temperatures and albedos as recorded 
by the emissive and reflective sensors of the HCMM. We designed a 
number of these (vegetative cover, plant water stress, ambient air 
temperature, distance from the coast, soil type, nighttime dewpoint 
temperature) into our objectives. Other variables also influenced our 
results, and it is the objective of this section to present an over- 
view of some of these sources of variation. 


The data used were pooled from the representative sites that were 
chosen especially for synoptic analyses (Chapter 3), weather station 
sites that were needed for freeze hazard analyses (Chapter 4), and 
intensive sites where ground truth were taken by individual fields for 
the water stress (Chapter 6) and plant cover (Chapter 7) studies. 


From the 94 possible sites for each HCMM overpass date, we deleted 
any that were urban (nonagri cultural) , those that contained water, and 
all sites on a given overpass date that contained any cloud or sub- 
visible ci rrus -contaminated pixels in the 9x9 matrix for the site. Con- 
sequently, if the 9x9 matrix was deleted so was the 7x7, 5x5, 3x3, and 
lxl matrix for that site and date. Eight summer dates, 27 May through 
15 Aug 78, were included. The restrictions put on the data resulted in 
168 observations for each of the following variables: HCMM tempera- 

tures (T1...T9) and albedo (RT...R9), respectively, for ground areas 
consisting of lxl, 3x3, 5x5, 7x7, and 9x9 pixels; saturation deficit 
of the air (SQEF.mb) measured at the central environmental station be- 
tween 1300 and 1400 hr CST; distance of each site from the Gulf coast 
(DIST.km); insolation the day of the pass (ID0P,lv day! ) and hour of 
the pass (IHROPJy); water holding capacity (WCAP) of the rootzone of 
the dominant soil type; maximum air temperature (TMAX,°C) the day of 
the HCMM overpass; time since rainfall or irrigation (DSR.days); and 
square root of days since rainfall or irrigation (DSR) a. Temperature 
was used as the dependent variable in step-wise multiple regression an- 
alyses. (A ridge regression analysis would have been better but that 
analysis capability was unavailable.) 
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2. 3. 1.1 Ground sample size 


O 

Figure 2.5 summarizes the percent of the sum of squares (R xlOO 
where R is the multiple correlation coefficient and is the coeffi- 
cient of determination) accounted for as the HCMM ground sample 
increased from the single central pixel (36 ha) to a sample of 81 
pixels (9x9 pixel sample=2900 ha=29 km2). The regression coefficients 
for the first five terms in the regression model were significantly 
different from zero (t-test) but coefficients for additional terms were 
not. To learn whether the R values were different as the sample size 
increased, we compared the R value for five variables in the regression, 
for the 81 pixel sample (0.936) with the R value for the central pixel 
sample (R=0.906). The test used was 

Z 1 - z 2 



where z=!iln(|~) 

and z-j identifies R for the 9x9 pixel sample, 

Zg identifies R for the lxl pixel sample, and 
n-j and n^ each equaled 164. 

Calculated t was significant at the 102 probability level. Thus, 
the same relationships seem to be discerned regardless of the HCMM 
sample size that was used. 

The values of IDOP, IHROP, SDEF, DIST, DSR, (DSP)’*, Tmax, and WCAP 
were common for all ground sample sizes. The means and standard devia- 
tions (S.D.) for these variables ware: 


Vari ables 



IDOP 

IHROP 

SDEF 

DIST 



ly day“^ 

ly 

mb 

km 


Mean 

408.6 

74.3 

28.8 

101.5 


S.O. 

46.6 

24.0 

7.1 

33.1 



DSR 

(QSRp 

TMAX 

WCAP 



days 

(days) 2 

°C 

cm 


Mean 

17.9 

3.74 

36.4 

16.2 


S.O. 

17.5 

1.96 

2.1 

1.4 
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The means and standard deviations of HCMM temperature and albedo for the 
various size samples were: 


Sample size 


Pixel configuration: 

lxl 

3x3 

5x5 

7x7 

9x9 

Number of pixels: 

1 

9 

25 

49 

81 

Mean (°C)^ 

33.00 

32.95 

32.90 

32.87 

32.84 

S.D. (°C) : 

4.60 

4.48 

4.41 

4.38 

4.37 

Mean R (%) : 

17.29 

17.30 

17.33 

17.31 

17.28 

S.D. R(%) : 

2.21 

2.01 

1.90 

1.84 

1.78 


The means and standard deviations presented showed that mean temper- 
ature and its standard deviation decreased slightly but steadily as number 
of pixels in the ground sample increased from 1, to 9, to 25*... to SI. 

The albedos went through a maximum for a sample of 25 pixels, but the 
standard deviation decreased slightly but steadily as the sample size 
increased. The decrease in standard deviations as sample size increased 
seems to explain why the R 2 values in Figure 2.5 increase as sample size 
increases . 

These results seem to indicate that a single HCMM pixel can represent 
ground conditions almost as well as a larger sample, and that the observed 
HCMM temperatures are influenced by variables such as insolation and air 
temperature to such an extent that these synoptic variables may be as 
strongly related to the observed temperatures as field to field variations 
in ground cover and water conditions. In Chapters 6 and 7 we used ground 
samples of 1 (central pixel) and 9 pixels because the ground truth and 
vegetation indexes represented individual fields or adjacent fields and 
we wanted the closest correspondence between sample sizes we could get in 
data from different observation sources. Inclusions of variables such as 
ground cover or vegetation greenness and evapotranspi ration rates in the 
analyses consistently raised the proportion of the sums of squares 
accounted for in the regressions by about 15?!. 

2. 3. 1.2 Saturation deficit, distance from the coast, 
insolation, and soil type 

The first five variables included in the multiple regressions were 
SDEF, DIST, IHROP, I OOP , and WCAP . Saturation deficit of the air and 
the two insolation variables are weather related, distance from the coast 
is geographic, and water holding capacity is a soil property. The four 
variables not included were Tmax, R, DSR, and (DSR) 1 *. 
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Table 2.2 presents the linear correlation coefficients between HCMM 
observed temperatures for each ground sample size and each of the varia- 
bles. The simple correlation between Tmax and HCMM temperatures was 
highly significant (r=0.686** to 0.738** for the various ground sample 
sizes). However, Tmax was highly significantly correlated with SDEF, 
DIST, and IHROP, and they rather than Tmax were included in the multi- 
ple regressions. Days since rainfall, DSR, and (DSR)^, were most 
highly correlated with SDEF and Tmax. The albedos (R) from HCMM were 
highly significantly correlated with insolation and distance from the 
coast, as well as the HCMM temperatures. 

2 . 3.1.3 Cloud contamination 

To determine the effect of percent cloud contamination (PCM) on the 
observed HCMM temperatures, we added the observations that contained 
clouds and vapor contamination to the cloud-free data of sections 2.3.1 . 1 
and 2. 3. 1.2 above for the 9x9 pixel samples. Mumerical values of cloud 
contamination were obtained by counting the number of pixels that were 
determined to be cloud-affected using man-machine interactive procedures, 
dividing that number by 81, and multiplying by 100 to express the propor- 
tion as a percentage. Then we analyzed the data for each of the 8 
summertime dates for which we had data. 

The regression equations established for each of the dates, the 
number of observations, and the r- values are summarized in Table 2.3. 

The average regression coefficent for the 8 linear equations is -0.093 
which indicates that each 10% increase in cloud contamination lowered the 
observed HCMM temperatures by 0.9 C. 

The distribution of data points for two of the dates (17 Jun and 
03 Jul) is illustrated in Figure 2.6. The temperature variability for 
the data points with low cloud cover was real and was associated with 
differences in surface cover, land use, distance from the coast, and 
days since rainfall included in the data set. Labeling of data points 
for the 27 May data by land use showed, e.g., that most high tempera- 
tures represented rangeland sites; these had low vegetative cover and 
were very dry. 

Generally, the cloud contamination effect on HCMM temperatures 
indicated that the weaker effects we were interested in would be found 
only by using noncontaminated pixels. Consequently, in later chapters 
only pixels judged to be unaffected by cumulus or cirrus clouds are 
included in the analyses. 
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Figure 2.2 The upper curves of each graph show the unfiltered and filtered daytime temperature profiles of east-west 
transects across the test site on latitude 26.55°N. Two partly cloudy and two clear dates are illustrated. 
The loxtfer curve of each graph shows the ratio of the temperature of the center pixel to the mean temperature 
of a 9 pixel by 9 pixel moving window. 
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Figure 2.3 Scatter diagrams of ref lective(VIS) vs emissive (IR) 
pixel values on 12 June 78 and 23 June 78* The wide scatter 
is attributable to the individualistic nature of clouds. 
Clustering over cloud-free land and water areas is shown by 
the $ symbols representing reflective/emissive combinations 
that recurred at least 17 times in the 100,000 pixel test area 
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HCMM daytime scenes of the test site on 15 Aug 78 (upper) and 03 Jul 78 
(lower). In the 15 Aug scene, scattered cumulus clouds paralleling the 
coast (right and lower center) are apparent in both VIS and IR images 
while an atmospheric absorbing-emitting layer (wispy dark areas left 
center to upper right) are apparent only in the IR band scene. In 03 
Jul scene, the dark globs and striations in the IR scene are clouds 
that are more apparent than in the VIS scene. 
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Figure 2.5 Percent of sum of squares (R xlOO) accounted for as variables 
were added to multiple regression analyses for ground samples 
of 1, 9, 25, 49 and 81 RCMM pixels. 
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Table 2.1 ^cations applied to reflective (VIS) and emissive (IR) digital counts to specify which data 
should be discarded as cloud contaminated. 

Overpass Date 


Equations 


27 Hay 78 
06 Jun 78 
12 Jun 78 
17 Jun 78 
23 Jun 78 
03 Jul 78 

08 Jul 70 
15 Aug 78 
00 Nov 78 


If IR<68, or VIS>13.9 + 0.518 IR 

If IR<75, or VIS>14. 4 to 0.552 IR 

If IR<63, or VIS>37.9 + 0.304 IR 

If IR£70, or VIS>-5.0 + 0.768 IR 

If IR<77, or VIS>44.2 + 0.360 IR 

If IR<6 6 , or IR<97 and VIS>_31.56 + 0.248 IR, or 
IR>97 and VIS>-20.28 + 0.776 IR 

If IR<41, or VIS>16 .52 + 0.41G IR 

If IR<69 , or VIS>-177. 80 + 2.645 IR 

If VIS>-52.174 + 1.304 IR 

If IR<50 and VIS>-84.94 + 2.359 IR, or 
IR>50 and VIS>23.76 + 0.208 IR 


07 Feb 78 



Table 2.2 Linear correlation coefficients between temperature and climatic, 
soil, geographic and time variables for HCMM sample sizes of 1, 9 
25, 49 and 81 pixels(Tl, T3, T5, T7, T9). 


Variable 

T1 

T3 

T5 

T7 

T9 







SDEF 

** 

.738 

kk 

.758 

r -J- — - “ 

ft* 

.763 

** 

.765 

** 

.768 

DIST 

** 

.721 

** 

.729 

** 

.725 

** 

.725 

** 

.721 

IHROP 

** 

.600 

** 

.613 

** 

.616 

** 

.614 

kk 

.610 

ID OP 

.128 

.126 

.115 

.109 

.101 

WCAP 

* 

-. 170 

* 

-.171 

k 

-.180 

* 

-.194 

-.205** 

TMAX 

** 

.686 

kk 

.712 

** 

.721 

kk 

.732 

.738** 

R-^ 

* 

.180 

kk 

.276 

** 

.296 

kk 

.276 

** 

.245 

DSR 

** 

.472 

kk 

.481 

** 

.485 

kk 

.492 

** 

.500 

(DSR)^ 

** 

.453 

kk 

.486 

kk 

.490 

** 

.498 

** 

.507 



27 May 

T9=31.45-0.073(PCN) 

60 

.214 

06 June 

T9=28. 40-0. 064 (PCN) 

54 

** 

.496 

12 June 

T9=30.21-0.105(PCN) 

64 

** 

• 53 &** 

.622 


T9=31 . 39-0 . 226 (PCN)+0 . 002 (PCN) ' 

64 

17 June 

T9-31. 72-0. 085 (PCN) 

46 

.498 

23 June 

T9=30. 86-0. 082 (PCN) 

73 

** 

.371 

03 July 

T9=32. 20-0. 105 (PCN) 

51 

** 

.505 

08 July 

T9=28 .05-0. 139(PCN) 

24 

.229 

15 August 

T9=33. 50-0. 101 (PCN) 

30 

.406* 



3.0 SYNOPTIC ANALYSES 

3.1 Abstract 

A test area 240-km across and 130-km N to S at the southern tip of 
Texas was used for various synoptic thermal investigations. Isotherm 
lines at 2 C intervals adequately displayed the thermal patterns within 
the test area. The temperature range -11 to +40 C generally encompassed 
the temperature range encountered for land features, water bodies and 
clouds. 

Digitization of the soil association map for a four county area 
each 20-seconds of latitude and longitude corresponding to the pixel 
resolution of HCMM enabled stratification of HCMM albedo and tempera- 
ture data by soil associations. Soil association temperatures from 
•10 daytim overpasses were typically highly correlated. Summer night- 
time temperatures tended to be uniform across the test area and close 
to the dewpoint temperature which is reached each night. On winter 
nights with strong radiational cooling and dew points lower than the 
air reaches, there was a temperature gradient across the area that was 
dominated by distance from the coast. Soil type and land use had much 
weaker influences. Warm areas adjacent to inland water bodies were 
observed in both summer and winter nighttime thermal data. 

The gradient of rainfall across the test area, that averages 1 cm/ 
4.5 km annually, prompted us to calculate an "aridity index" defined as 
the difference between daytime HCMM temperature and maximum air temper- 
ature for the day at the site of the weather stations. The aridity 
index for 10 dates increased with distance from the coast. The dif- 
ference between the HCMM-indicated nighttime surface temperature and 
the minimum air temperature showed that HCMM-indicated temperatures 
were within 2 C of dewpoint temperatures on nights when the air reaches 
the dewpoint temperature. 

The day temperature transect across the test area for 07 Feb 79 
and the night temperature transect for 26 Feb 79 were nearly mirror 
images of each other, and showed that ground sites that warmed most in 
the daytime also cooled the most at night. No rain fell between ob- 
servation dates, vegetation growth was controlled by low temperatures, 
and the atmosphere was dry on observation dates. The data indicated 
that when test area and atmospheric attenuation conditions were similar, 
day minus night temperature data separated in time by 2 weeks or more 
was meaningful. 

Overall, these studies demonstrated that data with approximately 
1-km ground resolution such as that from the N0AA-6...9 series of satel- 
lites contain useful information about synoptic thermal conditions. 

They can be used directly to measure thermal responses such as freezes, 
but they must be interpreted along with other relevant data for in- 
ferences about near surface soil water conditions, evapotranspiration, 
and crop condition. 
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3.2 Introduction and objectives 

The energy balance at the earth's surface is: 

Rn = H + LE + 5 

where Rn = net radiation flux received by the surface 
H = sensible heat flux to the air 
LE = latent heat flux (evapotranspi ration) 

G = sensible heat flux to the soil 


Since the equation expresses the principle of conservation of energy, 
it is clear that sensible heat (H) becomes an increasingly important 
product of the transformation of radiant energy at the earth’s surface 
(Rn) when soil water depletion limits latent heat production (LE). In 
the south Texas study area daytime air temperatures increase inland in 
relation to increasing aridity (annual rainfall decreases going inland 
at about 1 cm/4.5 km over a large portion of the study area). The pre- 
vailing wind is from the Gulf of Mexico toward the mainland and evidently 
helps moderate the inland surface and air temperatures. One objective 
of the synoptic analyses was to document the daytime increase of radio- 
metric surface temperature with distance inland. 


Increased surface temperatures is the mechanism by which reportion- 
ing of transformed energy takes place. Under the influence of midday 
insolation the surface is the warmest feature and heat moves from the 
soil surface into the soil by conduction and into the air by convection 
and radiative transfer. Stressed vegetation, which is above air temper- 
ature, also contributes to heating of the air. Because of vertical 
mixing with cooler air above, the air temperature lags surface tempera- 
ture temporally and in amount. At midday the surface of dry soil is 
typically 20 C above air temperature. Thus the drier the landscape 
condition, the greater is the expected difference between surface 
temperature and bulk air temperature. 


The question to be investigated was: Can the degree of aridity be 

represented on a synoptic scale by comparing Weather Service air tempera 
ture measurements (daily maximum) with HCMM indicated daytime surface 
temperatures in the vicinity of the weather stations? This was inves- 
tigated in the climatic gradient that exists between the Gulf of Mexico 
and weather station locations up to 170-km inland. We coined the term 
"aridity index" for THCMM-Tmax. 


Another objective was to relate nighttime HCMM temperatures to 
ambient air dew point temperatures at the Weslaco Environmental Station 
and to minimum air temperatures at 14 weather stations located from the 
Gulf of Mexico to 170-km inland. Hygrothermographs usually maintained 
at weather stations are not very accurate so there would be measurement 
error in humidity and temperature values taken from their charts. 
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Dew is important in agriculture because it provides a wet surface 
on which pathogens can develop and because it delays the onset of tran- 
spiration after sunrise. Our study area conditions are often favorable 
for nightirae dew formation (high humidity, light wind, clear sky) except 
during the winter months when the air is sometimes very dry. On nights 
of dew formation the landscape surface will usually be a few degrees 
colder than air temperature at shelter height. This temperature in- 
version results from reduced air temperature near the surface associated 
with radiation heat loss and the lower dew point temperature near the 
surface caused by removal of moisture from the air by condensation. 

The minimum air temperature (usually at dawn) may be closely 
correlated with surface temperature at HCMM overpass time. This was 
investigated for a summer night with data collected at the Weslaco 
environmental station. The high moisture content of Gulf of Mexico 
air (vapor pressure of approximately 28-mb during half the year) 
limits the air temperature decline during dew formation in the late 
nocturnal hours (relatively great amount of dew formation, with 
attendant heat release as the ambient air temperature approaches 
the dewpoint temperature) . 

There were other minor objectives including the determination of 
thermal islands in the vicinity of inland water bodies and cities, and 
comparisons of day and night temperature transects across the test area. 

3.3 Methods 

Procedures described in steps a through s of section 2.1.1 were 
applied in meeting the objectives of this chapter. The HCMM indicated 
surface temperatures were corrected by a radiative transfer model (J. C. 
Price, personal communication, 1980; Barnes and Price, 1980). Radio- 
sonde data from Brownsville, Texas (0600 CST or 1300 CST as appropriate) 
were used in determining the corrections which were applied uniformly 
to the entire test area. It is recognized that ideally the corrections 
would be determined at overpass time at several locations and interpolated 
over the area. 

Surface temperature gradients from the coast to 240-km inland were 
determined for the ten daytime HCMM scenes. The profiles were prepared 
from parallel east-west transects 20-km apart. 

To study the aridity index the maximum air temperatures measured 
at 14 weather stations on the dates of HCMM overpasses were paired with 
radiometric surface temperatures derived from the HCMM emissive band for 
various HCMM ground sample sizes, ranging from 0.3 to 29-km , centered 
on weather station coordinates. 

For nighttime HCMM scenes the HCMM radiometric temperatures were 
compared with ambient air dewpoint temperature recorded at weather sta- 
tions and on one night with the minimum air temperature recorded at the 
central environmental station. 
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3.4 Results and discussion 
3.4.1 Thermal contouring 

One of our objectives was to determine if there are locations in 
the LRGV that have consistently different temperatures from surrounding 
areas. One possibility of detecting areas that are consistently dif- 
ferent from surrounding areas is to produce and examine temperature 
contour maps. 

Each of 10 scenes of the LRGV was stored on a computer disc 
registered as described in 1.2.7. A computer program was modified to 
permit drawing of contour lines over any portion of the slightly more 
than 100,000 pixels in each scene. Due to memory limitations of our 
computer, individual pixel data were contoured by breaking the scene 
into sections 32 pixels by 64 pixels, contouring that section, then 
matching adjacent sections until the LRGV was covered. Lesser detail 
was obtained by averaging groups of pixels in the array. 


We found that isotherm lines drawn at 2 C intervals gave adequate 
detail to show the temperature patterns in the Valley. We used three 
colors of ink in the plotter pen to distinguish one isotherm from an- 
other. Temperatures from -11 to 40 C were contoured where present. 
Isotherms representing temperature multiples of 5 C were cross-hatched 
to assist in reading the maps. Most scenes had more than a six degree 
range of temperatures so isotherms drawn by any color could represent 
temperatures six degrees apart. Cross-hatching of the isotherm multi- 
ples of five, and the sequence of colors made it possible to identify 
the temperature represented by each isotherm. 

Temperature isotherms, especially at the single pixel scale, showed 
the location of clouds. Cloud tops were cooler than surface temperatures 
Large cloudy areas were also evident in maps of blocks of pixels, but the 
temperature gradients were not as steep as for the individual pixel con- 
tours; the steepness of the gradient depended upon the size and arrange- 
ment of the clouds within the block of pixels. 

The VIS channel data were also contoured. Rather than converting 
the VIS digital counts to the engineering unit, albedo, the contour 
lines were drawn with a digital count interval of eight between adjacent 
contours, or 24 counts between contours of the same color. Again the 
contours that were multiples of five were cross-hatched to simplify 
reading the resulting maps. Because of the difficulty and expense of 
reproducing colored maps, they are not included in this report. 

3.4.2 Soil association temperatures 

A soil association map, based on Soil Conservation Service soil 
surveys (section 1.2.7) was digitized and stored on computer disc file 
so that each area corresponding to a satellite pixel was given a symbol 
representing the predominant soil in the pixel. This soil map was 
registered with the HCMM data so there was a soil association symbol 
for each pixel location of the reflective and emissive data. Table 3.1 
summarizes the information for the soil association map of the LRGV area. 
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Because the two sets of data (satellite and soil association) were 
registered on a pixel by pixel basis, it was possible to determine the 
albedo and the surface temperature of each soil association pixel. The 
number of pixels of each temperature from each soil association was 
counted and a histogram developed. The temperature histograms of the 
clear scenes showed a sharp peak for each soil association. The tempera 
ture at the peak of the temperature histogram was used to represent the 
temperature of the soil association. For clear scenes this value con- 
tained the median temperature for each soil association. The term 
median temperature was used to refer to the temperature at the peak of 
the temperature histogram for each soil association; only in a few 
cloudy scenes was the true median temperature not in the histogram 
cell at the peak of the histogram. 


Clouds over a scene gave combinations in a scatter diagram of 
reflective and emissive data such that-few pixels fell in the same 
scatter diagram cell. Data from ground surface, on the other hand, 
yielded a high frequency of occurence of reflective and emissive data 
pairs in a tight cluster (section 2. 2. 1.7). 

Figure 3.1 and 3.2 are scatter diagrams of reflective and enissive 
data for a clear and a cloudy scene, respectively. This characteristic 
of clouds and ground surface to group, or to be scattered, was used to 
separate cloud-covered pixels from pixels with a clear view of tie 
ground surface. Any pixel falling in a scatter diagram cell containing 
17 or more pixels was considered to be cloud-free; fewer than 17 data 
pairs in a cell indicated cloud-contaminated pixels. Table 3.2 lists 
the number of cloud-free pixels in each of 36 soil associations for the 
ten dates. According to this criterion 27 May, 8 Jul, and 15 Aug 78 
were cloud-free. 

Because clouds are not uniform in either surface temperature or 
brightness, temperature histograms for cloudy soil association pixels 
peaked at essentially the same temperature as did the cloud-free pixels. 
Table 3.3 shows the minor differences from calculating soil association 
temperatures using only cloud-free pixels and using all pixels. The 
three dates included in Table 3.3 were the three with the fewest cloud- 
free pixels in Table 3.2. The small differences in these three dates 
indicated that other dates with fewer cloudy pixels would have iless 
temperature differences from using all pixels or excluding the cloudy 
pixels. 

The clustering method of distinguishing between cloudy and cloud- 
free pixels requires that data from both reflective and emissive chan- 
nels be available. This limits the procedure to daytime passes of the 
satellite. The median temperature for each soil association can be 
determined from the emissive data from either a daytime or a nighttime 
pass. Because the histogram peak temperature was found to be very 
similar using both all pixels or only cloud-free pixels, the tempera- 
ture histogram procedure can be used when only emissive data are 
available. 


Temperatures determined by the histogram peak method for the 36 
soil associations on the 10 daytime dates are given in Table 3.4. 

Linear correlation coefficients for all combinations of soil association 
temperatures for the 10 dates (n=10) showed that most temperatures were 
significantly correlated statistically. The fact that surface vegeta- 
tion affected emissive measurements may have masked possible differences 
in soil associations. 

Soil temperatures at the nightime HCMM overpass tended to be 
uniform across the test area in summer (Table 3.5). This is illus- 
trated by Figures 3.3 and 3.4 which show the mean temperature of all 
pixels within each of the soil associations. The standard deviations 
related to them are indicated by dots in the Figures. The symbols 
representing the soil associations in the figures are identified in 
Table 3.5. The W's represent water temperatures of the Gulf of Mexico, 
Laguna Madre, and Delta Lake. 


In the LRGV dew frequently forms on vegetation at night, and 
instrument shelter air temperatures approach the dew point within a 
degree or two. The liberation of heat in dew formation tends to re- 
tard further temperature decline and results in relatively uniform 
air and surface temperatures through the area. This is indicated 
for most dates listed in Table 3.5 and is illustrated for the night 
30 Jun 78 by Figure 3.3. 

The mean temperatures and standard deviations of HCMM pixels 
falling in the grouped associations on the 30 Jun 78 night are illus- 
trated by Figure 3.4. The grouped loamy upland soil associations 
(identified as 3, consisting of C, D, E, F, G, H, I, J, and of 
Figure 3.3), and saline clays and loams (identified as 5, consisting 
of U, X, and Y) were about half a degree warmer than the other soils. 

Variable amounts of soil visible to the satellite, climate 
characteristics typical of the various soil associations, and their 
different thermal inertias were among conditions causing these results. 
However, the overriding factor appears to have been nighttime dew 
formation over the entire test area. 

Figure 3.5 shows the soil association temperature means and standard 
deviations for cool season conditions of 21 Oct 78. There was a pro- 
nounced decline of temperature with distance inland from the Gulf of 
Mexico. There was no dew on this night, however, because the dew point 
was much below ambient temperature. 


3.4.3 Daytime temperature increase from the coast 
3. 4. 3.1 Transects 

Surface temperature transects were made across the test area 
(Figure 1.3) using HCMM data for all 10 daytime dates studied. An 
example of summer daytime conditions is shown in Figure 3.6, represent- 
ing 15 Aug 78. The east ends of the transects originated in the Gulf 
of Mexico (right side of figure) and traversed inland to west longitude 
99.5° (left side of figure). These parallel east-west transects were 
separated from each other in north-south direction by 0.2° latitude. 

Each transect is identified (right margin) by the north latitude which 
it represents. 

Figure 3.6 shows that the surface temperature of the Gulf of Mexico 
was uniform with little change from pixel to pixel, or for the entire 
length of the transect across the water on 15 Aug 73. (The temperature 
transect representations of the five latitudes have been offset from one 
another in the figure by 10 C temperature increments to facilitate exam- 
ination of the individual traces. The temperature scale on the left side 
of the graph applies to the temperature trace made at 26,01° north lati- 
tude.) . 

A trend of increasing temperatures with distance inland, from right 
to left, was apparent in all daytime traces that we prepared except for 
27 May 78 and 06 Jun 78, when there was little overall change across the 
test area. 

In Figure 3.6, the relatively cool temperatures at 99.00° W near 
the western (left) end of the 26.21° N transect were caused by the water 
in the Marte Gomez reservoir in Mexico (Figure 1.3). This water was at 
about the same temperature as the Gulf of Mexico. 

Cloudy conditions bounded by 97.7° N to 98.0° N resulted in the 
irregular traces on transects 26.41, 26.61 and 26.81° N. Areas of data 
removed from analyses by automatic cloud screening were "bridged across" 
by the computer in making these traces. The temperature dip at 97.95° W 
on transect 26.41° N was real and represented the water temperature of 
Delta Lake (Figure 1.3). Similarly, the cooler temperatures on transects 
26.61 and 26.81° N between 99.2 and 99.4° W were caused by Falcon Reser- 
voir (Figure 1.3). 

The surface temperature transects of a clear winter day are illus- 
trated in Figure 3.7. These temperatures are 15 to 20 C cooler than the 
summertime temperatures discussed above. On this date, 07 Feb 79, the 
Gulf of Mexico's temperature illustrated by the rightmost portion of 
each trace was colder near shore (Padre Island) than farther out. This 
contrasts with summertime traces which showed rather uniform Gulf tem- 
peratures. 
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The two southernmost transects (26.01 and 26.21° N) showed increasing 
surface temperature with distance inland (right to left). The transects 
through the middle and northern part of the test area showed a positive 
temperature gradient immediately inland from the coast, but generally 
constant or even declining temperatures in the western half or two-thirds 
of the land area represented. 

As in the summertime, water bodies were cooler than daytime tempera- 
tures of land areas. The water bodies (Gulf of Mexico, Delta Lake, 

Marte Gomez, and Falcon Reservoir) (Figure 1.3) are apparent in the 
temperature traces. 

3.4. 3.2 Ground sites 

The data set used is the same one as analyzed in section 2. 3. 1.2. 

For this purpose, however, the linear regression equation for each of 8 
dates that related the temperature observed by HCMM for 9x9 pixel sites 
on the ground (T9, °C) to distance from the coast (DIST, km) was deter- 
mined. The results given in Table 3.6 showed that r2 was significant for 

5 of the 8 dates. In general, the proportion of the variation in tempera- 
ture accounted for by site distance from the coastline was poor for the 
early sampling dates, went through a maximum for 23 June, and then de- 
clined --although remaining highly significant— through 15 Aug. 

The 27 May data were taken near the end of a long rainless period, 
whereas the 06 Jun data were taken following rains on Jun 03-05 that 
should have left the soil surface wet at all sites. Another long, essen- 
tially rainless period lasted through the 15 Aug observation date. If 
the same sites had been used on all dates, a continuously increasing tem- 
perature with time after 06 Jun would be expected. However, only sites 
left after cloud-contaminated data were deleted were used so that number 
of observations and their distance from the coast was unique for each set. 

The intercepts of the regression equations of Table 3.6 are the 
estimated temperatures at the coastline and the coefficients in front of 
DIST give the rate, °C knH, at which temperature increased with distance 
from the coast. The average coefficient for the 8 dates was 0.05 so that 
there was an average 5 C temperature difference between the coast and the 
mean distance inland of the sites (100 km). The coefficients increased 
steadily from 06 Jun to 03 Jul as soil water was depleted and seasonal 
evaporative demand increased. 

Based on the expected steady increase in temperature expected after 

06 Jun, the 17 Jun and 03 Jul intercepts or coastline temperatures are 
notably low (Inland sites have reasonabe temperatures for 03 Jul because 
of the large regression coefficient, 0.09 C km-1.). The coastline temper- 
ature on these two dates was too low to be accounted for by inclusion of 

a disproportionate number of irrigated sites in the sample (the tempera- 
ture difference between means of irrigated and dryland sites was not that 
large; see Table 7.3). Precipitable water in the atmosphere that atten- 
uated the emissions, lower temperatures due to shadowing of the sites by 
cumulus clouds just preceding the HCMM overpass, and undetected cirrus 
clouds are also possible explanations for the low coastline temperatures. 
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The lower part of Table 3.6 provides information on the mean tempera- 
ture and distance from the coast for each date and their standard 
deviations. The mean distance from the coast for sites sampled on 17 Jun 
was only 56.8 km as compared with an average of 107 km for the other seven 
dates. The temperature standard deviation generally increased as the dis- 
tance from the coast increased, probably because the sites farthest from 
the coast were dry rangeland sites with highest temperatures. 

Figures 3.8 and 3.9 show mean and standard deviations of surface 
temperatures at "representative sites" that consisted of 3x3 pixels. 

The letters designate the type of land use at the sites as follows (crop- 
land may be fallow except in summer): 

A Sand 

C Cotton 

G Grass 

K Sugarcane 

M Mixed crops 

R Rangeland 

S Grain sorghum 

T Citrus 

U Urban and industrial 

W Water 

Z Wildlife refuge 

The positive and negative standard deviations are noted by dots above and 
below the letters representing the sites. All of the 94 HCMM study sites 
that were not cloud covered are shown in the Figures. 

Figure 3.8 shows the situation on 17 Jun 78, a day characterized by 
a moderate inland climatic gradient and more than the usual spread of val- 
ues. Figure 3.9 represents 08 Jul 78 when a strong climatic effect 
resulted in a more marked temperature gradient from the ocean. 

It had been 9 days since rainfall in the test area on 17 June and 10 
days on 08 Jul. Factors which may have contributed to the varying results 
were the presence of 26% cloud cover over the test area at HCMM overpass 
on 17 Jun compared with 3% on 08 Jul, and daily atmospheric transmissiv- 
ities (fraction of solar radiation which survives passage through the 
atmosphere to the earth's surface) measured at the Weslaco Experimental 
Station of 0.53 and 0.68 on the respective dates. Also, the mean satur- 
ation deficits were 12.0 and 15.8 mb, and the calculated potential eva- 
potranspiration (Kincaid and Heermann, 1974) were 7.0 and 9.1 rnm/day on 
the respective days. 

Striking temperature differences betweer land uses are absent although 
some land uses did maintain their relative relationships with one another. 
For example, sugarcane (K) tended to be below average temperature, and 
rangeland (R) was generally warmer than cropland. 

Examples of winter temperatures at the representative sites are given 
in section 3.4.7. 
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3.4.4. Nighttime thermal islands in the vicinity of water bodies 

Nighttime surface temperatures in the vicinity of Delta Lake, La Sal 
Vieja and Sal del Rey are shown in Figure 3.10 for 26 Feb 79. The first 
two lakes are located and identified by name in Figure 1.3. Sal del Rey 
is the small lake in the upper left of Figure 3.10. The area represented 
in this surface temperature map measures about 40 by 40 km. The scale 
distortion is caused by the line and character spacing of the line printer 
that made the illustration. 

The nighttime surface temperatures shown in the figure were made to 
fall into three classes: warm, medium, and cold. Each class contains 
approximately one- third of the total pixels as determined from a frequency 
distribution of temperatures. In this case, the warm designation (darkest 
symbol) represents 3.44 C or warmer, cold (lightest symbol) is 2.23 C or 
colder, and medium is the range in between. Although the water bodies were 
the warmest features of the scene their areas have been left blank to 
facilitate identifying their locations. 

The heat island effect surrounding the lakes on this winter night was 
apparent in Figure 3.10. This was the coldest clear night (free of sub- 
visible cirrus clouds) for which we have HCMM data. The warm areas some 
distance to the southeast of the lakes (lower right) are probably unrelated 
to the water bodies. They may have resulted from a high proportion of fal- 
low land and possibly high water tables. Bare soil, especially if moist, 
has warmer nighttime temperatures than low growing vegetation (Geiger, 

1965; Monteith and Szeicz, 1962; and Nixon and Hales, 1975). The night- 
time temperature decline was not impeded by dew formation on this night; 
when HCMM overpassed, the average dew point depression at 1.5 m above 
the ground was 6.7 C at the Weslaco Environmental Station (Station 9, 

Figure 1.3). 

The uniformly cold northern (upper) portion of the illustration is 
rangeland, outside of the irrigated area. A dry vegetated area would be 
expected to be the coldest part of the landscape in contrast with fallow 
soil (especially if moist) or irrigated crops. The area is level to 
gently rolling so topographic relief was not a major factor in establish- 
ing nocturnal temperature patterns resulting from air drainage. 

The heat island effect was also apparent in Figure 3.11 where HCMM 
surface temperature patterns of two cold nights (21 Oct 78 and 26 Feb 
79) are combined. The area represented by Figure 3.10 is located in the 
upper middle portion of the figure. The outlines of the three lakes 
help identify the area. 

The heat island effect was also apparent during a summer night, 30 
Jun 78, illustated for the three lake area by Figure 3.12. The darkest 
symbol represents temperatures of 20.54 C and warmer, the lightest symbol 
19.50 C and colder, and the medium toned symbol the range in between. 

There was a relatively uniform distribution of temperature across the 
scene, except in the vicinity of the lakes. Dew had probably formed in 
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most of the area by the time of the HCMM overpass, and it probably caused 
the uniform temperatures . The dew point depression measured at the 
Weslaco Environmental Station (Section 9, Figure 1.3) during the hour in 
which HCMM overpassed was 0.9 C at 1.5 m above the ground. This is indi- 
cative of dew conditions (Jensen, 1974). 

3.4.5 Synoptic aridity index 

Even at midday, canopy temperature of a fully transpiring crop is 
close to air temperature, unless a strongly advective situation prevails. 
However, as soil water becomes limiting and evapotranspi ration is impeded, 
the balance of energy transformation is changed. Sensible heat becomes 
important and is associated with surface temperatures that are signifi- 
cantly above air temperature. This is especially true at the soil 
surface (at midday dry soil is typically 20 C above air temperature). 

Thus, the drier the condition, the greater is the difference between bulk 
air temperature and surface temperature. The question investigated was: 
Can the deqree of aridity be represented on a synoptic scale by comparing 
Weather Service air temperature measurements (daily maximum) with HCMM 
indicated daytime surface temperatures in the vicinity of the weather 
stations? 

Figures 3.13, 3.14, and 3.15 are given to examine this question. 

They show the difference between the daytime surface temperatures (°C) 
of 9x9 pixel footprints surrounding the weather stations and the maximum 
air temperatures (°C) recorded on the day of the HCMM overpass. We call 
this an aridity index (THCMM-Tmax) . 

The symbols showing the aridity index in the figures are identified 
in Table 3.7 as to their weather station locations. The table also gives 
identifying numbers that show geographic locations in Figure 1.3. 

The HCMM data used in this section have been corrected by a radi- 
ative transfer model (J. C. Price, personal communication) using 1800 
CST radiosonde data from Brownsville, Texas, applied to the entire test 
area. The temperature correction for each date is noted in Table 3.8. 

Figure 3.13 gives the aridity index on 27 May 78 following a 34-day 
dry period. It shows that the aridity index increased with distance 
inland from the coast, as land use changed from irrigated agriculture 
to rangeland, and went from an area of coastal showers to a more arid 
interior. Annual rainfall decreases at about 10 mm/4. 5km from the coast 
over a large portion of our test area. The predominant land uses within 
the HCMM 9x9 pixel footprints surrounding the weather stations are given 
in Table 3.7. 

The meteorological conditions for the dates represented by Figures 
3.13, 3.14, and 3.15 were: 


Date 


Percent Solar Daily Mean Potential 

clouds radiation, atmos. saturation evapotranspi ra 
at overpass ly/day trans- deficit,mb tion, mm/day 

mission 


27 May 78 

12 

651 

0.70 

11.0 

8.6 

06 Jun 78 

20 

628 

0.66 

8.4 

7.7 

23 Jun 78 

25 

616 

0.64 

12.4 

8.4 


The percent cloud cover at HCMM daytime overpass was determined by analy- 
sis of the VIS (reflective) coverage of the test area in conjunction with 
the IP* (emissive) data. The percentages do not include areas, if any, 
that were covered by subvisible cirrus clouds (SCi). The daily solar 
radiation and the other entries were measured at the Meslaco Environmen- 
tal Station (station 9, Figure 1.3). The daily atmospheric transmission 
is the ratio of measured daily solar radiation (horizontal surface) to 
the theoretical extraterrestrial radiation on the day under consideration. 
The mean saturation deficit is the integrated daily difference between 
the actual vapor pressure and the saturation vapor pressure at the exist- 
ing temperature. Daily potential evapotranspi ration was calculated from 
meteorological records using the Penman (1953) method as applied by 
Kincaid and Heermann (1974). 

The condition on 06 Jun 70 is shown in Figure 3.14. Rain had 
occurred at all locations during one to four of the preceding five days, 
averaging 33 inn, but ranging from a minimum of 4 mm at Port Isabel (H in 
Figure 3.14 and station 16 in Figure 1.3) to 73 mm at Raymondville 
(G and 11 in the figures) and Zapata (A and 1) (Point "A" for Zapata is 
missing from these figures because the northwest corner of the test area 
was truncated in processing the scene.). Potential evapotranspi rati on 
on these three dates was above normal (Nixon and Smithey, 1979) for their 
respective times of the year. Recent rainfall was apparently the cause of 
the small aridity index on 06 Jun 78. 

Aridity indexes for 23 Jun 78 (Figures 3.15) resulted from 14 
essentially rainless days (A few light showers occurred at stations 
nearest the coast between 17 and 21 Jun, with the maximum cumulative 
amount of 12 mm recorded at Raymondville (G and 11). The high aridity 
indexes in the inland portions of the test area can be attributed to 
soil moisture deficits that accumulated as a result of successive days 
of high evaporative demand without moisture replenishment. The steep- 
ness of the best fit line of Figure 3.15, compared with Figure 3.13, 
showed that more arid conditions prevailed inland on the June date than 
27 days earlier in May. The best fit lines (dashed) were determined 
using only land-bound weather stations. Stations A, H, L, and <1 were 
excluded from regression calculation on this and other dates because of 
nearby water (Laguna Madre and Falcon Lake; see description of station 
surroundings in Table 3.7). The water surfaces in the 9x9 pixel foot- 
prints surrounding the stations could distort the land representations. 
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Table 3.8 lists the equations of the aridity index best fit lines 
with respect to distance inland (km) for all daytime HCMM data processed. 
On one date, 08 Jul 78, the radiative transfer model was not used be- 
cause of missing radiosonde data. Of necessity, the correction for that 
overpass was estimated and applied to the intercept. 

The regression coefficients of aridity indexes listed in Table 3.8 
reflect the influence of decreasing available soil moisture with time 
in the interior rangeland, and the increased evaporative demand from 
spring to summer. In general, the coefficients increased from May to a 
maximum on 08 Jul. The relatively small regression coefficient calculated 
for 15 Aug was likely caused by subvisible cirrus cloud contamination of 
data, as discussed in Chapter 2 of this report. 

The positive regression coefficients in Table 3.3 showed that the 
aridity indexes were frequently greatest in the interior rangeland. 

This is expected considering the decline of rainfall from the coast 
(annually 10 mm/4.5 km inland), and the fact that irrigated surroundings 
predominated within 90 km of the coast. 

The standard deviations of aridity indexes from the best fit lines, 
listed in Table 3.8, averaged 2.88 C for the 10 dates. The standard 
deviations are in °C because the aridity index is calculated by subtract- 
ing air temperature, °C, from satellite indicated surface temperature, 

°C. The lowest values occurred on 06 Jun 78 following a rainy period 
and on 07 Feb 79 a cool winter day (9.7 C mean). Sy comparison the 
standard deviations of aridity indexes from best fit average were 2.35 C 
for nights of unimpeded radiational cooling (21 Oct 73 and 26 Feb 79) 
and 0.90 C for nights of dew formation (30 Jun 78, 23 Aug 78, and 
10 Oct 78). 

The deviation of daytime aridity indexes from the best fit line with 
respect to distance inland is given in Table 3.7 for individual weather 
station locations. Some locations had consistently high aridity indexes 
while others tended to be below the trend. This is indicated in the 
Table's right column where the number of negative to positive occurrences 
are given. 

Locations that usually had positive aridity index deviations were 
McCook (C), Mission (D), and McAllen (K). These values above the best 
fit line may have resulted from either unusually hot surface temperatures 
in the 9x9 pixel footprint, or cold air temperatures. The urban sur- 
roundings of McAllen would contribute to a positive deviation (White, 
et al. 1978). 

Negative aridity index deviations were observed at Raymondville (G) 
and Las Escobas Ranch (K) . These could have been caused by unusually 
cold surface or hot air temperatures. 
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Other causes of either positive or negative aridity index devia- 
tions could be an inaccurate thermometer, a poorly located instrument 
shelter, or a nonrepresentative microclimate, lie recognize, however, 
that the aridity index we calculate is not the true temperature differ- 
ence between the air at instrument shelter height and the land surface 
at the time of the satellite overpass. 

Analysis was made of 365 consecutive days (10 May 78 through 09 May 
79) of thermograph records from the Weslaco environmental station 
(station 10, Figure 1.3, representative of the irrigated area) and Las 
Escobas Ranch (station 25, Figure 1.3, representative of rangeland). 

In both the surroundings, the maximum air temperature was usually warmer 
than at nominal HCMM overpass time; within 1 C over half the time, and 
within 3 C during 90^ of the time. 


The HCMM data do not measure landscape surface temperature precisely; 
they need correction for atmospheric effects and surface emissivity. We 
corrected the daytime HCMM data used in this section by applying a radi- 
ative transfer model ( J. C. Price, personal communication) with radio- 
sonde measurements from Brownsville, Texas (station 17, Figure 1.3). 

To represent actual conditions the corrections listed in Table 3.8 should 
have been higher because the air temperatures were higher (resulting in 
more correction) at satellite overpass than at 1800 CST when the obser- 
vations we used were obtained. Also, we applied the Brownsville radiative 
transfer corrections to the entire test area, but in reality the correc- 
tions should have varied with distance inland due to increasing air 
temperatures and decreasing humidity in the lower atmosphere. 


Our HCMM data have not been corrected for emissivity. During the 
summer much of the irrigated area (eastern portion of test area) was 
covered by adequately watered vegetation. The rangeland (inland western 
portion of test area) had vegetation that sparcely covered the soil in 
many locations, and the soil surface was generally dry. Thus, the tenper- 
ature corrections for emissivity would have been greatest in the range- 
land (Fuchs and Tanner 1966, 1968; Sutherland and Bartholic, 1977). 

A change in emissivity of 0.01 corresponds to a temperature change of 
about 0.7 C (Phinney and Arp, 1975). 

The factors discussed above (higher actual surface temperature and 
lower actual air temperature at satellite overpass) gave greater differ- 
ences than expresssed by our aridity index. Also, the differences 
between surface and air temperatures would increase more with distance 
inland if emissivity corrections were made. Even though the aridity 
index does not exactly represent the difference between surface temper- 
ature and air temperature at satellite overpass, it deserves further 
examination. 
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3.4.6 Nighttime HCMM temperatures vs ambient dewpoint 
temperatures 

Examples of nighttime surface temperatures in the presence of dew 
on 30 Jun 78 have been presented in the discussion of soil association 
temperatures (Figure 3.3) and thermal islands in the vicinity of v/ater 
bodies (Figure 3.12). Those illustrations show that generally uniform 
surface temperatures prevailed under dew conditions. The dew point de- 
pression averaged 0.9 C at 1.5 m above the ground at the Weslaco 
environmental station during the night hour in which HCMM overpassed 
on 30 Jun 78. 

The differences between HCMM indicated nighttime surface tempera- 
tures surrounding weather stations (average of 9x9 pixels) and the 
minimum air temperatures recorded at the weather stations are illus- 
trated in Figure 3.16 for 30 Jun 78. This representation of THCMM- 
Tmin was computed in the same manner as the aridity index described in 
section 3.4.5. The HCMM temperatures of the test area used here were 
adjusted by the radiative transfer model results obtained with 0600 CST 
radiosonde data from Brownsville, Texas. 

Figure 3,16 shows that THCMM-Tmin was positive at most stations, 
and it tended to increase with distance inland. The weather station 
identifications corresponding to the letters in the Figure are given in 
Table 3.7. 

Several factors are involved in the THCMM-Tmin results. Nocturnal 
cooling under a clear sky causes surface temperatures that are below air 
temperature because of heat loss by radiation. Hence, if it were not 
for other factors, THCMM-Tmin would be negative. Because landscape eiiiis- 
sivities range from 0.95 to 0.99 (bare soil to full canopy), the equivalent 
black body temperatures seen by the satellite are lower than actual sur- 
face temperatures, a further cause of negative THCMM-Tmin. 

Negative tendencies may be counteracted by the fact that minimum air 
temperature is used in calculating THCMM-Tmin, instead of the usually warmer 
air temperature at HCMM overpass time. Analysis of 365 consecutive days 
of thermograph records (10 May 78 through 09 May 79) from the Weslaco 
environmental station (irrigated area) and Las Escobas Ranch (rangeland) 
showed that during 85% of the nights, air temperatures at overpass time 
were within 3 C of the minimum for the night. 

Even though a radiative transfer model was used, the radiosonde 
data were obtained several hours after overpass, hence adjustments were 
not exact. Furthermore, even if the adjustments were exact at one loca- 
tion (perfect model and atmospheric data obtained at overpass time), they 
would not precisely apply to all pixels of the large test area. 
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The radiative adjustment of 3.61 C applied to the 30 Jun 7C night- 
time HCMM data appeared to be too large at the Weslaco Research Farm 
environmental station in an irrigated area. Weather information from 
the hour in which HCMM overpassed and HCMM surface temperature data from 
9x9 pixel footprints surrounding the weather stations are listed below: 



Envi ron . 
Station 

Las 

Escobas 

Avg. All 
Stations 

Dewpoint at overpass, °C 

22.11 

- 

- 

Air temperature at overpass, °C 

23.01 

24.64 

- 

Minimum air temperature, °C 

22.67 

23.19 

22.11 

THCMM, unadjusted, °C 

19.68 

20.20 

19.86 

THCMM, adjusted, °C 

23.29 

23.81 

23.47 

Std. Deviation (9x9 THCMM) ,°C 

0.91 

0.81 

0.78 

Std. Error (9x9 THCMM), °C 

0.10 

0.09 

0.09 


The situation shown for the Weslaco environmental station (first data 
column) is unlikely in that the adjusted THCMM exceeded air temperature 
at the nighttime overpass, 23.29 C vs 23.01 C. The humidity at this 
time was 94% at shelter height. The minimum air temperature of 22.67 C 
(hourly average) occurred between 0300 and 0400 CST. The THCMM-Tmin at 
the environmental station was 0.62 (23.29-22.67). 

The situation shown for Las Escobas Ranch in the rangeland area 
(second data column) is reasonable. The adjusted surface temperature 
was 0.83 C below air temperature. 

At the inland location of Las Escobas Ranch (station 25, Figure 1.3) 
a slightly lower dew point might be expected than at the Weslaco environ- 
mental station. The similarity of nighttime minimum air temperatures in 
summer in the central and western portions of the test area is in agree- 
ment with widespread nocturnal dew formation. The narrow range in 
dewpoint temperatures, in turn, resulted in a rather uniform distribution 
of HCMM nighttime surface temperatures . By inference the dew point (not 
measured) at Las Escobas was less than a degree below the minimum air 
temperature of 23.19 C (Geiger, 1965; Jensen, 1974) and the adjustment 
made to the HCMM surface temperature was within 0.5 degree of being 
correct . 
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The last two lines of the above information give an idea of the 
thermal variability of the nighttime scene. The standard deviation and 
standard error relate to the 81 pixel (9x9) HCMM footprint surrounding 
the weather station. The mean is the THCMM value. 

Averages of values obtained at all test area weather station loca- 
tions are presented in the right column, except that stations are 
excluded whose 9x9 HCMM footprint contained water of Laguna Madre or 
Falcon Reservoir. These data substantiated the earlier statement that 
surface temperatures are quite uniform over the entire test area when 
dew is present. 

The use of temperature data to the second decimal place is unreal- 
istic considering the magnitude of the measurement errors that were 
involved. However, decimal fractions of degrees are important to the 
usefulness of HCMM data, especially where differences among surface 
temperatures rather than absolute amounts are considered. 

The mean daily maximum and minimum air temperatures, °C, for five 
summer months for the period 1931-1962 (Orton et al . , 1967) are given 
in Table 3.9 for stations representing coastal, intermediate, and 
interior locations (see Figure 1.3 for station locations corresponding 
to the identification numbers in the table), tiote the similarity of 
minimum air temperatures at Weslaco and Rio Grande City. The table 
shows that the large inland diurnal temperature range is due to higher 
inland maximum temperatures. The wintertime diurnal temperature 
ranges are also greatest at inland locations, but they result mainly 
from lower minimums (Table 3.9) made possible by low dew points character 
istic of that time of year. 

3.4.7 Day and night temperature transects for 07 Feb 79 and 
26 Feb 79 


For our LRGV test area, it was not possible to obtain day/night 
HCMM coverage 12 hours apart, nor was it possible to obtain cloud-free 
36-hour day/night coverage. We did receive a day-IR scene for 07 Feb 
79 and a night-IR scene for 26 Feb 79. These two were clear winter 
scenes with weak atmospheric absorption, characteristic of that time of 
year. The data in these two profiles were smoothed with a 3x3 moving 
window [%) . Air temperatures averaged 6 C and 21 C for the night and 
day scene, respectively. 


Figure 3.17 shows that the 07 Feb 79 day scene and the 26 Feb 79 
night scene unadjusted temperature profiles across the LRGV at latitude 
26.2° N were almost mirror images; that is, high day land surface 
temperatures on 07 Feb 79, were generally mirrored by the low night 
land surface temperatures on 26 Feb 79. The temperatures of water in 
Falcon Reservoir and the Gulf of Mexico are almost identical on both 


dates, as expected for these deep water bodies. If 12-hour or even 
36 -hour day minus night HCMM coverage were available, the temperature 
difference between day and night land surface temperature profiles may 
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have been more pronounced. The Figure indicates, however, that day/night 
overpasses need not be that frequent if vegetative cover remains essen- 
tially static and little or no rain has fallen in the intervening period. 

The clear and weak atmospheric absorption characteristics of the 
07 Feb 79 and 26 Feb 79 day/night scene pair means that the temperature 
profiles shown are probably due to near surface soil water content and 
to vegetative cover. Surface temperature adjustment calculated by a 
radiative transfer model using Brownsville, Texas radiosonde data for 
1300 CST, 07 Feb 79 was 0.92 C and for 0600 CST, 26 Feb 79 it was -0.21 C. 
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Figure 3.1 Scatter diagram of emissive (X-axis) and reflective (Y-axis) pixels of clear (07 Feb. 79) satellite pass. The 
numbers represent the count of pixels falling in each cell of the diagram. Counts of 10 to 16 are coded A to G 
counts above 16 are designated by X. 
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figure 3.2 


Scatter diagram of emissive (X-axis) and reflective (Y-axis) pixels of cloudy (.27 May 78) satellite pass. The 
numbers represent the count of pixels falling in each, cell of the diagram. Counts of 10 to 16 are coded A to G 
counts above 16 are designated by X. 
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Figure 3.4 HCMM temperatures versus distance inland from the coast for 
six groups of soil associations with similar properties for 
the night of 30 June 78. 
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Figure 3.5 Soil association temperatures from HCMM versus distance from 
the coast for the night of 21 Oct. 78. (Standard deviations 
of temperature are indicated by dots above and below each 
soil association symbol.) 
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WEST LQNGITl'QE 


3.6 Surface temperature transects across the test site at each of 
five latitudes for the daytime overpass on 15 Aug. 78. 
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SURFACE TEMPERATURE CC) 


Figure 



WEST LONGITUDE 


3.7 Surface temperature transects across the test site at each of 
five latitudes for the daytime ot'erpass on 07 Feb. 79. 
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TEMPERATURE 



DISTANCE INLAND (KM) 


Figure 3.8 Surface temperature (3x3 pixel means and standard deviations) 
by 11 land use categories at 94 sites versus distance inland 
from the coast for the daytime overpass on 17 June 78, 
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TEMPERATURE 



DISTANCE INLAND (KM) 


Figure 3.9 Midsummer gradient of temperature inland from the coast by land 
use codes (identified by letter) .Standard deviations of temper- 
ature are indicated by dots above and below each land use symbol 
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Figure 3.10 


Nighttime surface temperature on. 26 Feb. 79 displayed by pixel 
as 1/3 warmest (£3.44 C, darkest symbol) , 1/3 coldest (>2.23 C, 
light symbol) and 1/3 intermediate (2.23-3.44 C, intermediate 
symbo.1)' for a portion of the test site. Water bodies left blank 
and outlined. *^*UOvAL FAGij }.$ 

OF POOR QUAIJTY 
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Areas of greatest emissions during nighttime overpass 
Oct. 78 and 26 Feb. 79. Inland water bodies are outlii 
the centers of cities and towns are identified by dot 









figure 3.12 
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Portion of scene displayed as in Figure 3.10 for the night of 
30 June 78. 
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ARIDITY INDEX 


27 MAY 78 DAY 
3X3 SITE 


B Rio Grande City 
C McCook 

D Mission 

E Research Farm 

F Weslaco 



DISTANCE INLAND CKM) 



;ure 3.13 Aridity index (T^ where T is maximum air temperature 

for the day at weather stations anO^^ is the temperature mean 
of a 9x9 pixel sample surrounding the s tationsv-the dashed line is 
the regression fit (see Table 3 . 8 ) . 
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Figure 3.14 Aridity index as in Figure 3.13 for 06 June. 78 
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Figure 3.15 Aridity index as in Figure 3.13 for 23 Jun© '78 
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Figure 3.16 The difference between HCMM indicated nighttime surf ace temperature 
and daily minimum air temperature versus distance inland for the 
night of 30 June 78. 
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Table 3.1 Summary of number of pixels, rank, percent of total, and number 
of acres represented (based on 600 x 600-m pixels) of each of 
the 36 soil associations in the Rio Grande Valley of Texas. 


Soil 

Code 

No. 

No. 

Pixels 

Rank 

% of 
Total 

Acres 

Soil 

Association 

45 

5483 

1 

IS. 351 

487,775 

MeAllen-Brennan 

54 

2505 

2 

7.013 

222,844 

Nueces-Sarita 

2 

2139 

3 

•5.989 

190,284 

water 

81 

1851 

4 

5.210 

165,554 

Se j ita-Lomaita-Barrada 

21 

1827 

5 

5.115 

162,529 

Catalina-Copita 

31 

1648 

6 

4.614 

146,605 

Delmita-Randado 

11 

1518 

7 

4.250 

135,041 

Laredo-Olmito 

*5 

1345 

8 

3.766 

119,651 

Rayraondville 

46 

1340 

9 

3.752 

119,206 

McAllen-Zapata 

13 

1328 

10 

3.718 

118,138 

Rio Grande -Hat am or os 

48 

1194 

11 

3.343 

106,218 

Willacy-Pharr-Hargill 

92 

1098 

12 

3.074 

97,677 

Willamar 

29 

- 1083 

13 

3.032 

96,343 

W illacy-Raymond vill e 

61 

984 

14 

2.755 

87,536 

Harlingen 

43 

952 

15 

2.665 

84,690 

• Hidalgo-Pharr 

42 

884 

16 

2.475 

78,640 

Pharr-Brennan 

23 

847 

17 

2.371 

75,348 

Hidalgo 

27 

809 

18 

2.265 

71,968 

Raymondville-Lyford-Willamar 

82 

719 

19 

2.013 

63,961 

Mustang- coastal dunes 

47 

704 

20, 

1.971 

62,627 

Willacy-Delfina-Hargill 

62 

645 

21 

1.806 

57,378 

Harlingen-Benito 

28 

626 

22 

1.753 

'5,688 

Willacy-Racombes 

55 

594 

23 

1,663 

52,842 

Sarita-Faifurrias . 

24 

565 

24 

1.582 

50,262 

Hidalgo-Kaymondville 

41 

537 

25 

1.503 

47,772 

Copita 

63 

437 

26 

1,223 

38,875 

Mercedes 

22 

354 

27 

.991 

31,492 

Lyford-Raymondville-Lozano 

71 

340 

28 

.952 

30,241 

Jimenez-Quemado 

12 

306 

29 

.857 

27,222 

Reynosa-Runn 

91 

297 

30 

.1832 

26,421 

Laredo saline-Lomalita 

26 

295 

31 

.832 

26,243 

Raymondville-Hercedes 

52 

231 

32 

}. ‘ 

.647 

20,550 

Comitas -Delmita 

44 

73 

33 

.204 

6,494 

McAllen 

53 

56 

34 

.157 

4,982 

dune land 

1 

51 

47 

46 

35717 

35 

36 

.132 

.129 

4,181 

4,092 

3,105,461 

spoil banks 
Comitas 




Table 3.2 Number of cloud-free pixel? (determined by the scatter diagram 
method) in each of the 36 soil associations on each of the 10 
satellite passes. 


Dates of satellite pass 



Soil 

1979 


— . 

- - 1978 

. _ — 

, 

— — 

.. — — . 



No. 

Code 

2/07 

5/27 

6/06 

6/12 

6/17 

6/23 

7/03 

7/08 

8/15 11/08 

1 

1 

A 7 

47 

34 

36 

44 

38 

33 

47 

47 

44 

2 

2 

2139 

2139 

1022 

703 

1808 

1624 

1738 

2139 

2139 

1825 

3 

11 

1518 

1518 

1076 

1072 

973 

1345 

1106 

1518 

1518 

1491 

4 

12 

306 

306 

303 

302 

56 

305 

302 

306 

306 

305 

5 

13 

1318 

1328 

1108 

1208 

911 

1228 

1176 

1328 

1328 

1314 

6 

21 

1827 

1827 

1785 

1483 

1289 

1501 

1720 

1827 

1827 

1826 

7 

22 

35A 

354 

275 

295 

212 

335 

270 

354 

354 

325 

8 

23 

8A7 

847 

836 

791 

54 

828 

732 

847 

847 

844 

9 

2A 

565 

565 

481 

416 

311 

462 

190 

565 

565 

565 

10 

25 

13A5 

1345 

1191 

1028 

755 

898 

502 

1345 

1345 

1338 

U 

26 

295 

295 

281 

234 

295 

251 

103 

295 

295 

292 

12 

27 

809 

809 

625 

778 

515 

753 

620 

809 

809 

795 

13 

28 

626 

626 

504 

382 

200 

370 

2 

626 

626 

614 

1A 

29 

1083 

1083 

856 

928 

613 

741 

423 

1083 

1083 

1012 

15 

31 

16A8 

1648 

1207 

1287 

1257 

992 

1342 

1648 

1648 

1635 

16 

A1 

537 

537 

517 

483 

443 

431 

515 

537 

537 

536 

17 

A2 

88A 

884 

861 

859 

218 

849 

822 

884 

884 

854 

18 

A3 

952 

952 

875 

847 

505 

874 

641 

952 

952 

952 

19 

44 

73 

73 

73 

71 

53 

70 

73 

73 

73 

63 

20 

45 

5483 

5483 

4945 

4790 

4007 

4300 

4343 

5483 

5483 

5025 

21 

A6 

13A0 

1340 

1072 

1134 

952 

866 

1034 

1340 

1340 

1317 

22 

A7 

704 

704 

571 

535 

625 

502 

245 

704 

704 

696 

23 

A8 

1194 

1194 

1048 

1018 

627 

957 

175 

1194 

1194 

1174 

2A 

51 

46 

46 

35 

32 

31 

23 

30 

46 

46 

46 

25 

52 

231 

231 

191 

210 

193 

171 

153 

231 

231 

228 

26 

53 

56 

56 

55 

56 

54 

56 

56 

56 

56 

55 

27 

5A 

2505 

2505 

1636 

1821 

1890 

1682 

1907 

2502 

2502 

250* 

28 

55 

594 

594 

327 

317 

458 

520 

477 

594 

594 

586 

29 

61 

984 

984 

810 

735 

277 

842 

845 

984 

984 

981 

30 

62 

645 

645 

487 

449 

296 

583 

417 

645 

645 

638 

31 

63 

437 

437 

340 

275 

143 

312 

24 

437 

437 

427 

32 

71 

340 

340 

339 

316 

251 

286 

311 

340 

340 

335 

33 

81 

1861 

1861 

1511 

1294 

1761 

1527 

1497 

1861 

1861 

1751 

3A 

82 

719 

719 

402 

248 

656 

630 

501 

719 

719 

545 

35 

91 

297 

297 

223 

150 

270 

202 

230 

297 

297 

285 

36 

92 

1098 

1098 

865 

1066 

945 

947 

981 

1098 

1098 

1085 

Total 

35707 

35717 

28767 

27649 

23948 

28301 

25536 

35717 

35717 

34505 



Table 3.4 Median temperatures of 36 soil associations in the Rio Grande 
Valley of Texas determined by the histogram peak, method for 
the ten daytime satellite passes (Atmospheric corrections of 
Table 1.3 not applied). 

^979 “1978- ~ — — ~ 

No. Code 2/07 5/27 6/06 6/12 6/17 6/23 7/03 7/08 8/15 11/08 



Table 3.5 Nighttime mean surface temperatures and standard deviations, C, of south Texas soil association 
areas indicated by HCMM during 1978 and 1979. 


Distance 
Inland , km 


Soil Associations 


30 Jun 


Mean S.D. 


Mean Temperature <1 Standard Deviation 
23 Aug 10 Oct 21 Oct 


26 Feb 


Harlingen 19.7 

Harlingen-Benito 19.6 

Mercedes 19.7 

Laredo-Olmito 19.6 

Reynosa-Runn 19.8 

Rio Grande-Matamoros - 

Catarina-Copita - 

Lyford-Raymondville 19.5 

Hidalgo 20.2 

Hidalgo-Raymondville 19.7 

Raymondville 19.6 

Raymondville-Mercedes 19.6 

Raytnondvil le-Ly f ord 19.3 

Willacy-Recombes 19.5 

Willacy-Raymondville 19.4 

Delmita-Randado 20 . 8 

Copita 

Pharr-Brennan 19.9 

Hidalgo -Pharr 20. 1 

McAllen 20.1 

McAllen-Brennan 20.5 

McAllen-Zapata 20.2 

Willacy-Delfina-Hargill 20.1 
Willacy-Pharr-Hargill 19 . 9 

C omit as 20.1 

Comitas-Delmita 20.1 

Neuces-Sarita 19.5 

Sarita-Falfurrias 19.5 


Mean 


S.D. 

(1) s 
.7 
.7 
.6 

( 2 ) 1 
.7 
.7 

1.7 

.5 

.5 

.6 

.5 

.5 

.5 

.5 

.5 

.5 

(3) 1 
.5 
.5 
.5 
.5 
.5 
.6 
.6 
.6 
.6 

(4) S 
.7 
.5 
.6 
.6 


Mean S.D. 
hrink-Swell C 
17.8 .6 

18.1 .6 
17.7 .5 

oams and claS 


Mean 

S.D. 

Mean 

S.D 

£ 

14.8 

.9 

2.3 

1.2 

15.4 

1.0 

3.4 

1.0 

15.5 

.6 

3.6 

1,0 

15.6 

.9 

3.9 

1.3 

14.1 

.9 

2.8 

1.1 

14.7 

1.3 

3.6 

1.4 

13.7 

.9 

3.3 

1.2 

15.5 

.6 

4.1 

1.0 

14.7 

.7 

3.1 

.9 

15.2 

.7 

2.8 

1.0 

14.6 

.8 

2.6 

.8 

15.1 

.7 

2.4 

1.0 

14.8 

.8 

3.1 

1.1 

15.1 

.7 

3.0 

.7 

14.7 

.7 

3.0 

.8 

14.2 

.9 

2.3 

1.0 

14.0 

.8 

3.7 

1.2 

13.2 

2.1 

3.3 

.9 

14.0 

16.1 

3.1 

1.0 

14,3 

.6 

3.7 

.7 

13.9 

1.0 

2.8 

1.3 

14.6 

1.0 

3.7 

1.5 

14.6 

1.0 

3.1 

2.0 

14.5 

.8 

2.9 

1.2 

13.8 

.7 

1.9 

.7 

14.4 

.7 

2.7 

.8 

13.7 

.8 

1.3 

1.0 

13.9 

.7 

1.3 

1.0 







Table' 3.6 Equations expressing dependence of daytime temperature 
(T9) observed by HCMM on distance inland from the coast, 
mean temperatures and distances, and standard deviations for 
eight dates in 1978. 


Rate 

Equation 


No. of Observ. 

r2 

27 May 

T(9)=31. 87 + 0.010 

(DIST) 

18 

0.075 

06 June 

T(9)=26.97 + 0.010 

(DIST) 

24 

0.022 

12 June 

T(9)=26.22 + 0.061 

(DIST) 

17 

0.491* 

17 June 

T(9)=22. 68 + 0.054 

(DIST) 

20 

0.220 

23 June 

T(9)=25. 29 + 0.067 

(DIST) 

10 

0.910** 

03 July 

T(9)=23.45 -f- 0.090 

(DIST) 

15 

0.832** 

08 July 

T(9)=28. 88 + 0.060 

(DIST) 

38 

0.604** 

15 Aug. 

T(9)=32.06 + 0.052 

(DIST) 

26 

0.633** 



Temperature, C 

Distance, km 



Mean 

S. D. 

Mean 

S. D. 

27 May 

32.9 

1.09 

104.3 

29.3 

06 June 

27.9 

1.72 

90.8 

26.1 

12 June 

32. 3 

1.98 

100 .4 

32.5 

17 June 

25.8 

1.23 

56,8 

10.6 

23 June 

32.2 

1.67 

102.8 

23.7 

03 July 

34.8 

2.91 

126.2 

29.6 

08 July 

35.4 

2.55 

109.0 

33.3 

15 Aug . 

38.3 

1.93 

119.0 

29.4 
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Table 3.7 Weather station identification, description of surroundings, and deviation of daytine Aridity Indexes 
from the trend with distance inland. 


Map 
I dent , 


Station name & description 27 Mai 
Zapata (rangeland, lake) 


Rio Grande City (range- 
land; irrigated agric.) 
McCook (dryland grain 
sorghum) 

Mission (irrigated field 
crops & citrus) 

Research Farm (irrigated 
agric . ) 

Weslaco (irrigated agric. 
& urban) 

Raymondville (irrigated 
agric. , urban) 


grassland) 

Brownsville (irrigated 
agric. , airport, urban) 


airport, urban) 
McAllen (irrigai 
airport , urban) 


land) 


06 Jun 12 Jun 17 Jun 23 Jun 03 Jul 08 Jul 15 Au£ 


Sign | ; 
recur-? 
rences 
minu s | . 

08 Nov 07 Feb plus |j 


-7.1 -6.6 


-2. 0 

- 

-4.0 

-1.0 

0.3 

0.8 

-3.1 

-7.2 

1.4 

-0. 3 

6/3 

1.3 

1.9 

4.7 

3.9 

1.0 

0.1 

2.1 

3.3 

4.2 

2.7 

0/10 

0.5 

0.3 

2.4 

0.0 

0.8 

-0.1 

1.1 

0.6 

0.5 

0.3 

1/8 

0.7 

-0.8 

1.4 

1.2 

-1.6 

-3.1 

-2.1 

-0.9 

-1.0 

-0.2 

6/4 

1.1 

1.1 

0.0 

0.6 

-0.3 

1.4 

3.5 

0.8 

-0.9 

-0.8 

3/7 

0.0 

-0.2 

-1.0 

-0.6 

-2.0 

-2.0 

0.7 

-3.0 

2.0 

1 

o 

7/2 

-1.2 

-2.1 

0.5 

2.2 

0.9 

2.6 

0.6 

0.0 

-7.5 

-10.5 

4/5 

-2.7 

-1.4 

0.5 

2.8 

0.9 

1.9 

-1.8 

-1.2 

-5.1 

-1.8 

6/4 

-1.0 

-0.5 

-3.2 

-5.4 

0.1 

- 

-4.3 

-0.6 

2.1 

0.5 

6/3 

3.8 

2.3 

3.7 

- 

2.9 

2.9 

5.9 

3.4 

1.2 

2.1 

0/9 

-0.2 

-1.2 

-0.6 

-2.0 

0.4 

-0.8 

-3.9 

0.7 

3.6 

-1.7 

7/3 

-3.5 

-2.6 

5.3 

-1.6 

-2.3 

-1.6 

-2.3 

2.9 

-4.5 

-2.2 

8/2 

2.6 

0.7 

1.9 

4.0 

4.3 

5.2 

2.8 

-0.2 

-5.4 

-5.9 

3/7 



Table 3.8 Linear best fit of Aridity Indexes (daytime 9x9 pixel HCMM minus 
maximum air temperature, °C) with respect to distance inland (X, 
km), and atmospheric temperature correction used each date. 


Atmospheric 
Correction. C 


27 

May 

78 

06 

Jun 

78 

12 

Jun 

78 

17 

Jun 

78 

23 

Jun 

78 

03 

Jul 

78 

08 

Jul 

78 

15 

Aug 

78 

08 

Nov 

78 

07 

Feb 

79 


Eauation 


-0.93 + 0. 0377 -X 
0.80 - 0.0144-X 
-1.00 + 0.0267.X 
-1.58 + 0.0346-X 
-0.73 + 0.0482 *X 
-5.36 + 0.0548-X 
-2.68 + 0.0755 *X 
2.02 + 0.0276 -X 
8.36 - 0.0173-X 
4.73 + 0.0014‘X 


Standard 

Deviation 


Estimated value, radiosonde data missing for this date. Resulting intercept 
of the equation is consequently an estimate also. 





Table 3.9 Mean daily maximum and minimut summertime and wintertime air 
temperatures (C), by months, for the period 1931 - 1962 (from 
Orton et al. 1967) for - a coastal, an intermediate, and an 
interior location. 


Month 

Brownsville 
Station 17 

Weslaco 
Station 10 

Rio Grande City 
Station 2 

May 

30.3/21.7 

32.2/21.1 

34.4/20.4 

Jun 

32.2/23.9 

33.8/22.9 

35.6/22.8 

Jul 

33.4/24.5 

34.9/23.4 

37.6/23.5 

Aug 

33.6/24.4 

35.1/23.2 

37.9/23.2 

Sep 

31.9/22.9 

33.1/21.8 

34.7/21.4 

Dec 

21.9/12.1 

22.8/11.2 

22.4/8.7 

Jan 

20.7/11.2 

22.0/10.4 

21.7/7.9 

Feb 

22.7/12.8 

24.4/12.2 

24,4/9.3 
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4.0 FREEZE HAZARD 


4.1 Abstract 

Areas of least freeze hazard in our south Texas test area were 
identified by computer processing of HCMM data from two cold nights, 

21 Oct 78 and 26 Feb 73. First, a temperature frequency distribution 
for individual pixels that comprised our test area identified all pixel 
locations where temperature was in the top third (Had usable data been 
available from other cold nights they would have been processed to de- 
termine the locations of warmest pixels on those nights also.). Then 
the locations of warm pixels from both dates were superimposed to form 
a map showing the areas that had least freeze hazard. A bold symbol 
was used to identify locations that were warm on both dates. This 
approach offers promise as a useful method of classifying land areas 
according to freeze safety. Because of variations in temperature pat- 
terns from date to date, information from several freeze or near freeze 
events should be used in determining the safest areas to grow tempera- 
ture sensitive crops. 

HCMM nighttime data from the serious freeze of 03 Jan 79 could not 
be validly analyzed for landscape temperatures because an absorbing- 
emitting atmospheric layer, probably subvisible cirrus (SCi), obscured 
surface emissions from view of the satellite. The SCi layer, which was 
not apparent to ground observers, retarded the loss of nocturnal ter- 
restrial radiation to space and slowed the cooling of the landscape 
and air that night. This ability of satellites to detect an atmospheric 
layer that blocks the loss of terrestrial radiation should have important 
application in real-time freeze forecasting. Knowledge of the presence 
of such a blocking layer is justification for revising upward the minimum 
temperatures forecast for a cold night. This information could save 
growers from undertaking expensive freeze protection measures 
unnecessarily. 

On the coldest SCi-free night for which we have HCMM data (26 Feb 79) 
the HCMM indicated surface temperature surrounding weather stations ap- 
proximated the minimum air temperatures observed at the stations. In 
other words, presuming all calibrations were correct, the differences 
between air temperature and surface temperature at the 0209 CST HCMM 
overpass were approximately equal to the air temperature declines that 
occurred between overpass time and the mini mums (usually at dawn). The 
cold-night temperatures of HCMM footprints surrounding the weather sta- 
tions adjacent to water bodies did not closely match minimum air temper- 
atures because of the mixed scene comprising warm water and cold land. 

Cold nighttime temperatures in valleys and depressions were visible 
in the HCMM thermal image. In two cold air drainage situations investi- 
gated, the minimum air temperatures at dawn were lower than the satellite 
indicated surface temperatures surrounding the stations at overpass time. 
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The nighttime standard deviation from mean of pixels comprising foot- 
prints surrounding weather stations were typically 9.5 C to 0.3 C for 
cropland and rangeland locations, compared with 0.2 C to 0.4 C for 
offshore (Gulf of Mexico) water. 

4.2 Introduction and objectives 

Cold damage to crops usually occurs on clear, calm nights when 
radiational energy loss from the landscape results in critical vegeta- 
tion temperatures. The low atmospheric moisture and lack of clouds 
usually associated with freeze events are ideal for observations by 
satellites. The HCMM overpass time is late enough to well establish 
nighttime surface temperature patterns. 

Protection of vegetable and fruit crops against winter and spring 
freeze damage is becoming more difficult and expensive because of higher 
cost heating fuel. This makes it important to avoid planting tempera- 
ture-sensitive annual and perennial crops in the coldest portions of 
local areas . 

The objective was to test the suitability of HCMM data for identi- 
fying areas of least: freeze hazard in south Texas. "Safe" warm areas 
and hazardous cold spots were to be examined for recurrence by comparing 
satellite data from all cold nights. Techniques developed in this study 
should be applicable to other areas where subtropical crops are grown 
(California, Arizona, Florida). Hie method may also be useful in selec- 
ting sites for deciduous fruits (prevention of bloom damage) and other 
temperature-sensitive crops. 

A further objective was to determine the relationship between 
minimum air temperature and the surface temperature surrounding weather 
stations at the HCMM nighttime overpass. 

4.3 Methods 

The HCMM data were processed according to steps a through u 
(excluding f, g, h, o, p and q) outlined in section 2.1.1. The surface 
temperatures used in this chapter were corrected by a radiative transfer 
model (J. C. Price, personal coinnuni cation) using 0600 CST radiosonde 
data from Brownsville, Texas applied to the entire test area. 

4.4 Results 

4.4.1 Identifying areas of least freeze hazard 

Figure 3.11 shows the portions of the four southmost Texas counties 
that were warmest on two cold nights (21 Oct 78 and 26 Feb 79). Tempera- 
ture patterns of additional cold nights would have been incorporated into 
the figure had more HCMM data been available. The freeze night of 03 Jan 
79 for which HCMM data are on hand could not be used because the surface 
temperatures were obscured by an absorbing-emitting atmospheric layer 
believed to be subvisible cirrus (see section 4.4.2) . 


In preparation of Figure 3.11 a map-grid pixel was considered to 
be warm if its temperatures fell in the top third of all test area 
pixels. The darkest symbol used in the figure represents locations that 
were warmest on both dates-, the lighter symbol indicates warmest on only 
one date. Blank areas ware never in the top third of temperature, 
except that water bodies which were the warmest feature of the landscape 
have been left blank to facilitate identifying their location on the map. 

The figure shows an extensive warm area in the immediate vicinity 
of the coastline (right side of figure). Apparent in these nighttime 
data are heat islands adjacent to inland water bodies (outlined). Other 
heat island?* surround dots that locate the centers of cities and towns. 
The map sfvsws warm areas in unexpected locations such as portions of 
Starr County (left side). 

Such analysis based on several freeze nights could guide in selec- 
ting areas to plant temperature-sensitive crops. This is provided other 
factors such as availability of irrigation water and/or freedom from 
soil salinity are also favorable. The final site based on radiometric 
temperature surveys should consider the effects that vegetation density 
(amount of exposed bare soil) and enissivity may have had on a classifi- 
cation at a particular site. 

4.4.2 HCMM-detected protective absorbing-emitting atmospheric 
layer during freeze night 

A radiation absorbing-emitting layer prevented accurate assessment 
of surface temperatures on the only freeze night for which HCMM data are 
available for our test area, 03 Jan 79. The atmospheric layer is be- 
lieved to have been sub visible cirrus (SCi). This was unexpected to 
ground observers because the sky was visibly clear, as it usually is on 
radiational cooling nights when energy loss from the landscape causes 
critical vegetation temperatures. The hourly weather records from the 
FAA McAllen, near station 21 (see Figure 1.3), also indicated the sky 
was clear. Vie had not observed absorbing-emitting layers in any of the 
five NOAA-2 through UOAA-5 cold night scenes of the test area that we 
had examined over a five year period. 

The ability of the emissive (IR) channel of HCMM to detect essen- 
tially invisible absorbing-emitting atmospheric layers is discussed in 
section 2.2.2 of this report. The HCMM overpass on the freeze night 
occurred at 0205 CST on 03 Jan 75. The portions of the scene covered 
by the absorbing-emitting layer had an indicated temperature everywhere 
the same, or essentially so; much more uni form than landscape features 
outside the influence of the blocking layer. 

Figure 4.1 infers that much of the test area was hidden from satel- 
lite view by a SCi layer during the HCMM nighttime overpass on 03 Jan 79 
In the figure, temperatures at representative sites (means of 3x3 pixel 
footprints) are plotted against temperatures of the same sites in the 
NOAA-2 data of the historic 21 Dec 73 freeze (Mixon, Phinney and Arp, 
1974; Mixon, Phinney and Sautreaux, 1978). 


The figure shows that the representative sites (from which cities 
and sites adjacent to water bodies were excluded) had a narrow tempera- 
ture range (-1.5 C to -2.4 C) on 03 Jan 79. These data were apparently 
dominated by the effect of the SCI layer. In contrast the satellite 
indicated temperatures for the same sites ranged from +2.0 C to -5.3 C 
during the 21 Dec 73 event, giving no suggestion of the presence of a 
SCi layer on that night. 

A comparison of atmospheric-corrected surface temperatures (means 
of 3x3 pixel footprints) experienced at the same representative sites 
on cold, clear, SCi-free nights is made in Figure 4.2. In this case, 

26 Feb 79 is compared with 21 Dec 73. Both dates have considerable 
temperature ranges. Even though several years elapsed between events 
the land use patterns have changed little in the interim. 

A consistent temporal temperature relationship among sites would 
have resulted in a well defined linear grouping of points that plot 
diagonally across the figure from cold to warm. This was not the case, 
although there was a weak trend in this direction. 

A further consideration entering into a comparison of the cold 
nights is that the NOAA-2 overpass was at 2197 CST December 73 in 
contrast to the HCMM overpass at 0209 CST on the later date. Thus the 
observations by satellite were made at different stages of nighttime 
landscape temperature development and this might have influenced the 
relationships among the stations. 

However, practical conclusions related to freeze protection cannot 
be drawn from the figure because 26 Feb 79 was not a serious freeze event 
and because the overpasses were not at the same time of night. In any 
case, verification by study of additional freeze events would be prudent. 
The tentative conclusion to be drawn from Figure 4.2 is that comparative 
freeze damage from site to site varies from one event to another, as in- 
dicated by the wide scatter of points. 

Meteorological data obtained on the night of 02-03 Jan 79 illustrate 
the effectiveness of an absorbing-emitting atmospheric layer in halting 
noctural temperature decline. Typically on clear, calm nights character- 
istic of radiational freezes, the air temperature continues to decline 
until a minimum is reached at dawn. In contrast, thermograph data 
(Figure 4.3) shows what happened during the night under discussion. This 
example is from station 9 (Figure 1.3) where there was a customary decline 
from the afternoon high to 2320. This drop of temperature from 4.4 C to 
-4.5 C in less than 9 hours forbode disaster to the subtropical crops and 
vegetation of the Lower Rio Grande Valley during the remainder of the 
night. The presence of clouds or a high dew point reading would have 
given hope of averting a continuing temperature decline. But visible 
clouds were not present and the dewpoint of -12.7 C was far below the 
prevailing temperature. 
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If the air temperature trend established before 2320 had continued 
throughout the night (dashed line of Figure 4.3), a minimum temperature of 
about -8.0 C would have been reached at dawn (Nixon 1974). However, the 
development of the atmospheric blocking phenomenon, invisible to the 
human eye but detectable in the iiCMM emissive scene, caused a reduction 
of the net- terrestrial radiation during the remainder of the night and 
a consequent air temperature increase. The surface air drift during 
this time was insufficient in itself to account for the marked change 
in temperature trend. Apparently the reduction of net-radiative energy 
loss from the surface resulted in a reapportionment of the soil heat flux 
to heating the air. 


It appears that the absorbing-emitting atmospheric layer started 
to affect air temperatures at some locations between 2200 and midnight 
on the cold night of 02-03 Jan 79. Table 4.1 shows that nighttime 
temperature trends were moderated or reversed at or before midnight at 
stations 3, 4, 9, 13, 17 and 25 (Figure 1.3), These stations have a 
west-northwest east-southeast alignment across the test area. The moder- 
ation was short lived at western stations 3, 4, and 25 as the downward 
decline resumed at these sites. Haddock and Connolly (1979) have 
reported the thermograph of station 4 for this night. At stations 5, 

11 and 21 on either side of the main alignment, the effects of radiation 
blockage were delayed until 0200 or later. The Rio Grande City area 
(station 2) benefited little, if any, from the atmospheric layer. At 
the time of HCMM overpass, the Vaso Marte Gomez, a reservoir in Mexico 
opposite Rio Grande City, was outside the influence of the atmospheric 
layer and showed clearly in the emissive scene. Delta Lake and La Sal 
Vieja in the vicinity of station 11 (Figure 1.3) could also be clearly 
seen on the emissive band imagery. 

Falcon Reservoir was obscured in the HCMM emissive scene. Minimum 
air temperatures near the reservoir were -3.3 C at station 1 and -5,6 C 
at station 23. At station 1 the minimum was 0.6 C -warmer than it had 
been the night before. These relatively warm air temperatures for loca- 
tions so far inland can be partly attributed to their proximity to Falcon 
Reservoir, but they were mainly due to the absorbing-emitting atmospheric 
layer. These stations reached mini mums of -6.7 C and -7.2 C, respective- 
ly, in the December 1973 freeze which was said to have been generally 
less severe than the January 1979 freeze. 

The cold night of 02-03 Jan 79 caused serious freeze damage to 
citrus, sugarcane, vegetable and ornamental plants, but the damage was 
not nearly as disasterous as it would have been had the radiation 
absorbing-emitting atmospheric layer not been present. It can be seen 
from Figure 4.3 that the number of degree-hours of air temperature below 
some critical value would have been greater in the case of uninhibited 
temperature decline (dashed curve). The respective areas under the -3 C 
line suggest the relative extent of damage that would have been sustained 
by crops with that critical temperature. 
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Sugarcane freeze damage studies 

It was hoped that the HCMM coverage on the freeze night of 02-03 
Jan 79 would provide opportunity to verify the applicability of thermal 
satellite data to estimating relative amounts of freeze damage sustained 
by crops. The presence of the absorbing -emitting atmospheric layer 
prevented the test for which preparations were complete. The Rio Grande 
Valley Sugar Growers Association participated in the preparation by pro- 
viding the coordinates of field loading sites so that the HCMM indicated 
surface temperatures of areas surrounding these sites could be correlated 
with sugarcane yields and mill processing records for these sites. 

Surface temperatures observed with the NQAA-2 satellite on the 
occasion of the 20-21 Oec 73 freeze had correlated well with sugar con- 
tents of the crop harvested following that cold night (Nixon, Phinney 
and Gautreaux 1973). These comparisons between freeze-night satellite 
indicated temperatures at field locations and subsequent sugar contents 
of mill processed Sugarcane showed that the freeze damage was least 
where the satellite-indicated surface temperatures were warmest. 

According to t-tests made on data from individual field locations 
the correlations between freeze-night temperatures and percent commer- 
cially recoverable sugar ’was significant at p = 0.10. Since the sugar 
content of freeze-damaged sugarcane declines as it stands unharvested in 
the field following the freeze, the F-ratio was also calculated to take 
into account the multiple correlation of harvest dates and satellite 
sensed temperatures with recoverable sugar. The F-ratio for these test 
conditions was significant at p = 0.05. When satellite temperatures ware 
adjusted for the additional temperature decline between overpass and the 
minimunis at dawn the data were significant at p - 0.01. 

Importance of absorbing-emitting atmospheric layer 

The ability of satellites to detect the presence of an atmospheric 
layer that blocks the loss of terrestrial radiation to outer space should 
have important application in real-time freeze forcasting because as soon 
as the nocturnal radiation loss from the landscape is blocked, the rate 
of air temperature decline lessens. Consequently, the presence of a 
blocking layer justifies revising upward the minimum temperature forecast 
for the cold night. In borderline temperature conditions a more accurate 
low temperature forecast could save growers from undertaking expensive 
freeze protection measures unnecessarily. 


Currently, polar orbiting weather satellites with a thermal band 
pass over the U.S. nightly. Since the thermal band responds to atmospheric 
layers such as developed on the night of 03 Jan 79, the potentita! value 
of real -time processing of their data in conjunction with freeze 'warnings 
and advisories for growers is apparent. The development and operational 
testing of such procedures might be carried out under the multiagency 
AgRI STARS program. Under this program episodic events such as freezes are 
of high interest as an early warning of crop losses . The data would also 
be useful to the sugarcane, citrus and other fresh commodity industries 
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who need to react quickly to freezes. As mentioned above, the sugar 
content of the sugarcane stalks decreases 'fastest in the fields most 
seriously damaged so that such data would be a valuable tool in re- 
prioritizing harvest schedules following freezes to minimize sugar 
yield losses. 

If multiple uses --such as real-time advisories to growers during 
the freeze, estimates of crop losses as they affect the market, and 
commodity industry quick reaction decisions —are all made from essen- 
tially the same data, then there is little doubt the effort can be 
cost effective. For the Lower Rio Srande Valley of Texas, there would 
only be three to five nights of interest in a normal winter. However, 
all parts of the data acquisition, transfer, interpretation and rapid 
dissemination system would need to be ready to function on about one 
day's notice. The data would not have sufficient resolution to make 
decisions about individual fields, but the gradient of temperature from 
the coast, relative temperature differences of production areas several 
kilometers apart and presence of any protective cloud or atmospheric 
vapor layer over any portion of the production area would be apparent. 

4.4.3 Comparison of Tmin at weather stations and THCHM for 
weather station sites on cold nights 

On clear, radiations! freeze nights it is the temperature of vege- 
tation that determines how much crop damage is sustained. Customarily 
the likelihood and extent of crop damage is estimated from air tempera- 
tures . 

Table 4,2 compares the minimum air temperatures (Tmin) measured at 
weather stations with atmospherically corrected HCT1 temperatures (THCMM) 
indicated for the footprints surrounding the stations. The station num- 
bering in the table corresponds to locations in Figure 1.3. The size 
of the footprint varies from one pixel (20-seconds longitude by 20- 
seconds latitude) to 3x2 pixels (3-minutes longitude by 3-minutes lati- 
tude ) . 

A study of this table for 26 Feb 79 shows that the minimum air 
temperatures which occurred several hours after overpass were roughly 
estimated by the HCT! indicated surface temperatures. In other words, 
presuming all calibrations ware correct, the difference between air tem- 
peratures and surface temperatures at liCTi overpass at 9209 CST were 
approximately equal to the air temperature declines that occurred between 
overpass time and the minimums (generally at dawn). 

‘.Jhile "on the average" the above statement is valid, the table 
shows that the relationshio varied from site to site. Especially notable 
is Port Mansfield (26) where the surface temperature exceeded air tem- 
perature by 4.5 C for the single pixel to 2.4 C for the 9x9 footprints. 
The close vicinity of water at this station caused the discrepancy of 
temperatures. As the footprints increased in size there was a continuing 
improvement due to the influence of increased land area (land was colder 
than water). 
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An opposite effect occurred at Las Escobas Ranch (25) where in- 
creasing the footprint size increased the average surface temperature. 
Cold air drainage caused this rangeland site to have the lowest minimum 
air temperature of any station. Examination of the USCS Las Escobas 
Ranch 7. 5 -minute topogranhi c quadrangle sheet shows the site to be on a 
flat area near the confluence of two ephemeral streams. The drainage 
ways apparently facilitated cold air drainage. As the footprints size 
was increased it included more sloping land, most of it higher in eleva- 
tion and outside the cold pocket, thus resulting in a higher average 
temperature. 

Apparently, cold air drainage was also the cause of the low minimum 
air temperature at the Rio Grande City (2) station. The location of 
this station, 5 km west of the city, is on the Rio Grande River flood 
plain near the toe of land that slopes from the upland. Table 4.2 shows 
the minimum air temperature at this site to be about 3 C lower than the 
satellite indicated temperatures at its overpass. 

Several examples of the effects of cold air drainage are visible 
in the HC'i'l photographic image of 25 Feb 70 night for western Starr 
County (left portion of Figure 1.3 and Figure 3.11). This area has much 
more topographic relief than the remainder of the test area. 

The degree of nighttime thermal homogeneity of the landscape is 
indicated by the standard deviations (second data -line of each site) and 
standard errors (third data line). The marked differences of temperature 
between water and land are reflected in the large standard deviations 
(>1.5 C) of footprints that included both water and land: Zapata (1), 

Port Isabel (16), Port Ians field (26) and the largest footprints at 
Falcon Dam (23). Several sites were thermally homogeneous , such as 
McCook (3), and had standard deviations less than 0.3 C. 
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Figure 4.1 Temperature of 3x3 pixel samples (T3, C) of representative sites, 
excluding cities and water bodies, for two cold nights (03 Jan. 79 
and 21 Dec. 73) . 




26 FEB 79 


Figure 4. 



21 DEC 73 


Temperature of 3x3 pixel samples (T3, C) of representative sites, 
excluding cities and water bodies, for two cold nights (26 Feb. 79 
and 21 Dec. 73). 
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Table 4.1 Air temperatures at weather stations in the Lower Rio Grande Valley test area during the freeze of 02 - 03 
. January 1979. 


Station 








Time 







Name 

No. 

1900 


2100 

2200 

2300 


0100 

0200 

0300 


0500 


0700 


Rio Grande City 2 
McCook 3 
Texan Gardens 4 
Mission, Bates Power 5 
Weslaco Environ. Sta. 9 


Raymondville 

n 

Adams Gardens 

13 

Brownsville 

' 17 

McAllen 

21 

Las Escobas Ranch 

25 


-Z.8 

-3.9 

-3.3 

- 2.8 

-4.2 

-3.9 

-4.4 

- 2.8 

- 1.1 

-4.8 


-5.U 

-b . L 

-o. / 

-/ .z 

-4.4 

-5.0 

-5 . 6 

-5.6 

-4.4 

-5.0 

-5 . 6 

-6.1 

-5.0 

-6.1 

-6.7 

-6.7 

-4.8 

-4.7 

-4.0 

-3.9 

-5.6 

-5.6 

-5 . 6 


-3.9 



-3.9 

-3.3 

-2.2 

-2.8 

-3.3 

-3.3 

-4.4 

-3.3 

-3.9 

-3.3 

-4.6 

-4.7 

-5.6 


Data for all stations except 9 and 25 are from Texas Citrus Mutual Freeze Report No. 1, dated 05 Jan. 79. Data were 
reported to the nearest degree F. 


Table 4.2 Minimum air temperatures (T ) and nighttime HCMM radiometric surface 
temperatures (first lines) aPweather station locations on 26 Feb. 79, 
standard deviations (second lines) and standard errors(third lines) of 
nultipixel footprints. 

Sta. Sta. 


No. 

Name 

T 

Min 
( S C) .... 

HCMM Radiometric Surface Temperature, C 
lxl 3x3 5x5 7x7 9x9 

1 

Zapata 

6.67 

7.95 

8.93 

9.30 

9.27 

8.95 





1,75 

2.40 

2.63 

2.85 





.58 

.48 

.38 

.32 

2 

Rio Grande City 

0.00 

2.83 

3.37 

3.54 

3.64 

3.61 





.40 

.61 

.77 

,80 





.13 

.12 

•11 

.09 

3 

McCook 

2.78 

1.20 

2.15 

2.23 

2.29 

2.10 





.50 

.63 

.54 

.61 





.17 

.13 

.08 

.07 

5 

Mission 

1.67 

2.42 

2.74 

2.86 

2.65 

2.68 





.63 

.74 

.75 

.75 





.21 

.15 

.11 

.08 

9 

Research Farm 

2.22 

3.23 

2.60 

2.39 

2.41 

2.38 





.90 

.73 

.69 

,70 





.30 

.15 

.10 

.08 

10 

Weslaco 

3.33 

3.23 

3.28 

2.94 

2.63 

2.72 





. 58 

.75 

.77 

.76 





.20 

.15 

.11 

.08 

11 

Raymondville 

1.67 

4.04 

3.65 

3.36 

3.07 

2.78 





. 98 

.96 

.90 

.85 





.33 

.19 

.13 

.10 

16 

Port Isabel 

11.11 

9.49 

10.86 

11.17 

11.22 

11.19 





1.51 

1.60 

1.86 

2.07 





.50 

.32 

.27 

.23 

17 

Bro&msville 

3.89 

3.23 

3.10 

3.23 

3.53 

3.74 





.60 

.76 

.85 

.80 





.20 

.15 

.12 

.09 

20 

Harlingen 

3.33 

4.44 

3.63 

3.55 

3.28 

3.33 





.85 

.94 

1.10 

147 





.28 

.19 

.16 

43 

21 

McAllen 

4.44 

3.23 

4.04 

3.80 

3.75 

3.74 





.49 

.57 

.62 

.76 





.16 

.12 

.09 

.08 

23 

Falcon Dam 

4.44 

3.64 

3.50 

3.63 

3.93 

4.21 


■ ■ 



.28 

.78 

1.04 

1.55 





.09 

.15 

. 15 

.17 

25 

Las Escobas Ranch 

-0.67 

1.61 

2.42 

2.74 

3.09 

3.25 





.70 

1.05 

1.20 

1.22 





. 23 

•21 

.17 

.14 

26 

Port Mansfield 

8.33 

12.85 

11.77 

11.20 

10.95 

10.77 





1.75 

2.07 

2.31 

2,55 


•:! 


5.0 PLANTING DATE ADVISORY 


5.1 Abstract 

HCMM radiometric temperatures on a clear night early in the 1979 
planting season (26 February), when there was little atmospheric atten- 
uation, were very close to minimum air temperature, 3.1 C. Minimum 
temperature at the 10-cm soil depth was about 10 C higher. The radio- 
metric temperature of the sandy uplands (1.4 C) differed from that of 
the loams and clays (3.3 C). Three variables— average daily air tem- 
perature, insolation at ground level, and extraterrestrial radiation- 
accounted for from 89 to 93% of the variation in maximum daily (Tmax), 
minimum daily (Tmin), and average daily (T5) temperature at the 5-cm 
soil depth in dry plots and for from 88 to 91% of the variation_for 
wet plots. Standard errors of the estimate of Tmax, Tmin, and T5 from 
the three environmental variables ranged from 1.5 to 1.9 C. 

Estimation of planting depth soil temperatures from radiometric 
temperatures appears: feasible, provided the relation between planting 
depth and surface temperature is understood. The estimates should be 
superior to estimates currently made from air temperature observations 
and can probably have coefficients to tailor them to major surface 
texture groupings (soil associations). 

5.2 Introduction and objectives 


The flow of heat in soils follows linear-flow laws, the same as 
for the flow of water and air in soils. The calories of heat Q trans- 
mitted per unit area per unit time through a soil in the x direction, 
where the temperature is T and the temperature gradient is dT/dx, is 
given' by: 



where K is defined as the thermal conductivity. The negative sign is 
used because heat flows in the direction of decreasing temperature. 

In soils K can vary considerably with moisture content and soil type. 

It can be shown (Kirkhum, 1964) that vertical heat flow is expressed by 


k 


a x 


it 


where the thermal diffusivity k=k/pc, and p is density and c is specific 
heat. 


In this manner diurnal and annual soil surface temperature fluctua- 
tions are transmitted to the deeper soil layers. Practical implications 
(Chirkov, 1979) of the above equations are that: (1) At a given depth 

in the thermal lly active soil zone the interval between minima and 
maxima is 24 hours and 12 months; (2) The amplitude of the temperature 
range decreases in geometric progression in relation to the arithmetic 
progression of depth; and (3) There is a lag in occurrences of minimum 
and maximum temperatures with depth. This time lag is directly pro- 
portional to depth. 


Soil temperatures are critical for germination of crop seeds. 

For example, it is recommended that seed of corn, sorghum, and cotton 
not be planted in the spring until average soil temperature at the 
2-inch (5-cm) depth reaches 55, 60, and 65 F (or 13, 15, and 18 C), 
respectively, and a warming trend is forecast for the subsequent five 
days. Coarse-textured soils (sands and sandy loams) are better drained, 
retain less water and warm up earlier in the spring than fine-textured 
soils (clays and clay loams). Thus synoptic soil radiant temperatures 
should relate generally to soil association patterns if antecedent pre- 
cipitation was similar and the emissivities are the same for all soil 
types. Then, because of the correlation between surface temperature 
and temperatures at shallow soil depths, one would expect radiometric 
surface temperatures to give guidance on adequacy of soil temperature 
for planting. 

In this task, we proposed to examine the nighttime soil temperature 
patterns in the HCMM data for one or more nights as spring planting time 
approached and compare them with the soil association patterns as mapped 
by the Soil Conservation Service. Also available were (a) daily maximum 
and minimum soil temperature at the 10-cm depth for four monitoring sites 
during the period 01 Jan through 31 Mar 79, which the NWS, NOAA, used for 
planting advisories, and (b) soil temperatures in wet and dry plots we 
maintained at the Research Farm, along with the environmental station, 
during the HCMM investigation period. The soil temperatures recorded by 
NWS observers could be compared with the HCMM temperatures, and the con- 
tinuous wet and dry plot and environmental station data enabled an examin- 
ation of the dependence of soil temperature at various depths on the am- 
bient environmental variables. 

Toy et al_. (1978) produced and reported linear equations relating 
monthly, seasonal, and annual air temperature (F)— which is widely 
available in climatic summaries--and soil temperature (F) at the 10-cm 
depth for eight widely separated sites in the U.S.. Mean monthly soil 
temperatures, based on data from all stations combined, could be estimated 
to +2.3 C from the air temperature data. Meikle and Treadway (1979) used 
fourth degree polynomials to estimate maximum and minimum temperature at 
the 10-cm depth from air temperature for 114 U.S. locations. For Weslaco, 
Texas (Willacy fine sandy loam soil) the standard error of the estimate 
of soil temperature was 2.4 C for the maximum and 2.0 C for the minimum. 
Bonn (1977) found that the coefficients relating radiometric temperature 
(equivalent blackbody temperature observed with a Barnes PRT) and thermis- 
tor probe temperatures are close to unity for bare soil but vary with the 
seasons for vegetated surfaces. The closest correlations between remotely 
sensed temperature was with surface soil temperature for bare ground and 
with air temperature at the half canopy height over vegetated areas. 
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5.3 Methods 


There were three sources of data for looking at soil temperature 
as it related to planting crops: (a) daily maximum and minimum soil 

temperature at the 10-cm depth for the period 01 Jan through 31 Mar 
79 from four sites within the test area as made by volunteer observers 
for the NOAA NWS office at Brownsvil le, (b) HCMM radiometric temperatures 
for the night of 26 February 1979, and (c) soil temperature for the first 
90 days of 1979 (except for a few days when equipment was down) as meas- 
ured with thermocouples in the small wet and dry plots maintained at the 
Research Farm. Of these data sources "a" and "c" are described somewhat 
more fully in section 1.2.8. 

The data for the 10-cm soil depth at four sites and the maximum 
and minimum daily air temperature for these same sites were tabulated, 
keypunched, and graphed for each site to display the time trends and 
the relations between air and ground temperature. 

The radiometric temperatures from HCMM were displayed in relation 
to soil associations by the same procedures as for other dates under 
Synoptic Analyses (see section 3.4.2). 

For the thermocouple data of the wet and dry plots, the methods of 
acquisition and reduction were as follows: 

The center areas (about 1.5 m dia.) of each of two soil plots 6.1 
m square and 6.1 m apart were instrumented with type T thermocouples at 
the following depths: 0.5, 5, TO, 20, and 100 crn. In order to reduce 

"stem conduction error" at the thermocouple junctions special wire con- 
figurations were incorporated which placed the first several centimeters 
of wire from the junctions in the junction isotherm. Each junction buried 
5-cm and deeper was located in a spiral of 20-gage lead wire totaling about 
50 cm in length (Figure 5.1 A). Because the thermal conductivity of copper 
is about 10 times that of the mating copper-nickel alloy about 45 cm of 
the wire in the spiral was copper and the remainder was the mating alloy 
wire. 

Because the vertical temperature gradient in the vicinity of the 
0.5-cm depth thermocouples was so steep even more precautions were 
taken to minimize stem conduction and to preserve the integrity of the 
thin layer of soil above the thermocouple junctions. For these thermo- 
couples the junction was part of a grid consisting of 0.13-mm dia. wire 
wound on a rectangular frame made of vinyl -insula ted, 20-gage copper 
wire as shown in Figure 5. IB. The 2Q-gage lead wire was joined to the 
grid wire outside the frame. An electrical insulating spray material 
was applied to all uninsulated wire surfaces in both junction designs. 

Each plot had 3 independent thermocouples at the 0.5-cm depth and 
3 thermocouples connected in parallel at the 5-cm depth. One thermo- 
couple was located at each of the remaining depths. In addition, two 
shielded, air-temperature sensing thermocouples were mounted on a nearby 
tower at .5 m and 1.5 m heights. All of these instruments terminated at 
a Brown, 20-point, strip-chart recorder located in an instrument shelter. 
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The Brown recorder was timed to scan the thermocouples once every 
hour. The "wet” plot was watered weekly if rainfall did not exceed 0.7 
of the amount of weekly pan evaporation less the sum of the antecedent 
rainfall as described by Linsley, et jH_. (1958) using a factor of 0.9. 

The "dry" plot received only rainfall. Data were transcribed from the 
strip charts to IBM cards preparatory to computer processing. 

5.4. Results 

5.4.1 Soil association and NWS observer temperatures 

Figure 5.2 summarizes the mean temperature, of all HCMM pixels 
representative of each soil association, and the standard deviation 
for each soil association in the HCMM data for 26 February 1979. The 
symbols used are decoded in Table 3.5. It is evident that the scatter 
of temperatures versus distance (km) inland of soil association centroids 
has a "V" pattern. Sand dunes (V symbol) on Padre Island and saline 
coastal clays have temperatures of about 7 C, At the bottom of the "V" 
pattern at about 60 km inland, the "P" and "0" represent the sandy up- 
land soil associations, Sarita-Falfurrias and Nueces-Sarita which had 
radiometric temperatures of about 1 C. T, I, and D at the upper right 
arm of the "V" with radiometric temperatures of about 4 C identify the 
loamy upland inland soil associations Jimenez-Quemado, McAllen-Zapata, 
and Capita, respectively. When the soil associations were grouped by 
soil surface texture into six categories and displayed as in Figure 5.3, 
it was obvious that the mean temperature (1.4 C) of the sandy uplands 
(category 4) differed from that of the loams and clays (category 2) 
which had a mean temperature of 3.3 C. The two categories are equi- 
distant from the coast and the standard deviations do not overlap. 

The shrink-swell clays (category 1) are a fraction of a degree cooler 
than the loams and clays (category 2) even though the shrink-swell 
clays are closer to the coast. 

By comparison, the trend line in soil temperatures at the 10-cm 
depth, as measured by NWS observers with indicating maximum and minimum 
thermometers was a minimum on 09 Feb at all four sites; Adams Gardens 
(12.5 C), Mercedes (13.9 C), Weslaco (13.3 C), and McCook (12.5 C). 
Thereafter the temperature increased an average of 0.21 C per day 
through the end of March. On 26 Feb, the mean daily 10-cm soil depth 
temperature estimated from the trend line was 16.6, 16.2, 17.3, and 
16.4 C at Adams Gardens, McCook, Mercedes, and Weslaco, respectively 
and the minimum averaged 12.7 C. Thus the soil was warm enough to 
plant corn and sorghum, but it was still about a week too early to 
plant cotton. 

The average minimum air temperature at these four sites on 26 
Feb was 3.1 which is very close to the mean radiometric temperature 
of Figure 5.2. 


In summary, although there are many sources of variation in the 
data (differing vegetative cover conditions, mixture of irrigated and 
nonirrigated sites) the data show that radiometric temperature of the 
sandy uplands (1.4 C) differed from the temperature of the loams and 
clays (3.3 C) on a night of strong radiational cooling and that the 
radiometric temperatures were close to the minimum air temperatures 
that night. 

5.4.2 Wet and dry plot temperatures 


Although soil temperature data was recorded hourly at 5 levels, in 
each plot from late in 1978 until the middle of 1980, only the complete 
data for the three days around the HCMM overpass date of 26 Feb 79 was 
reduced and evaluated in detail together with selected daily data points 
during the period from 01 Jan through 31 Mar 79. 

The 26 Feb, 3-day period data provides an example of how wet and 
dry soil temperatures varied as functions of depth and time, as shown 
in Figure 5.4 and 5.5. At the 100- cm depth the diurnal variation dis- 
appeared from the temperature traces leaving an almost constant average 
daily temperature increase of 0.1 C for the wet plot and 0.23 C for the 
dry plot for this 3-day period. The 20-cm traces exhibit time lags 
ranging from three to eight hours from the 0.5-cm temperature traces. 

Air temperature (shown in Figure 5.6) exhibits a more consistent lag 
of 1.5 to 2.6 hours from the 0.5-cm temperature traces. The variation 
in surface temperature changes seems to be the chief cause of the wide 
variation in the lag times of the 20-cm deep soil temperatures. 

Figures 5.6 and 5.7 illustrate how the wet and dry plot soil tem- 
peratures compared with each other at depths of 0.5, 5, 20, and 100-cm. 
The average difference between the wet and dry soil temperatures for 
this period is 3.7 C at 0.5-cm deep, 3.3 C at 5-cm, 3.3 C at 10-cm (not 
shown), 2.9 C at 20-cm and 0.9 at 100-cm. 

The low air temperature of 26 February, lower than any wet and dry 
plot soil temperature (Figure 5.6), may be lower than we would expect 
if it were measured over a large area of uniformly bare, hard-packed 
soil. Hard-packed, bare soils cool less radiatively than the loose- 
packed and vegetative-covered soils that surround our test area. The 
temperature of the air that we recorded was affected far more by the 
general character of the surrounding, up-wind terrain than it was by 
our small test plot area, 

5.4.3 Prediction of wet and dry plot 5-cm soil 

temperatures from meteorological variables 

Da_ily maximum (Tmax) , daily minimum (Tmin), average daily tempera- 
ture (T5) and daily maximum minus minimum temperature (M-M) at the 5-cm 
depth in each of the wet and dry plots, as measured by thermocouples, 
were regressed on the following ten environmental variables: 
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~X 


Ta 

ICDPD 


TW 

INSOL 


I/W 


IOA 


BP 

MDP 

SDEF 

SDW 


These eight analyses as well as an additional one, the average 
daily temperature difference between the wet and dry plots, are sum- 
marized in Table 5.T. From the simple correlation coefficients, r, in 
the upper part of Table 5.1, it is evident that: 

(1) The 10 "independent" variables are all significantly cor- 
related with maximum daily temperature (Tmax) and mean daily 
soil temperature at the 5-cm depth (T5) (Any r > 0.217 is 
significant at p=.05). Seven of the 10 independent variables 
are significantly correlated with Tmin and 6 are correlated 
with maximum minus minimum (M-M) temperature in both wet and 
dry plots. 

(2) Mean daily air temperature is the most influential environ- 
mental variable affecting Tmax, Tmin, and T5 for both wet and 
dry plots (The simple correlations range from 0.847 to 0.945). 
Variables containing air temperature, such as air temperature 
times wind (TW) , also yield highly significant correlations. 

(3) Variables containing insolation are more highly correlate^ 
with Tmax of the soil and (M-M) than with either Tmin or T5. 

(4) For a given variable pair, the correlations are about the same 
for both wet and dry plots. 


- average daily air temperature, C 

- average of current day and previous day's 
insolation, Ly day”' 

- air temperature times wind movement, °C km day“l 

- insolation, Ly day”^ (same as current day 
insolation in ICDPD) 

- ratio of insolation for current and previous 
day to wind movement current and previous day, 
Ly day” '/km day”' 

- insolation outside the atmosphere incident on a 
surface of unit area parallel to the earth's 
surface— extraterrestrial radiation—, Ly day”' 

- barometric pressure, inches Hg 

- mean dew point, °C 

- saturation deficit of the air, inches Hg 

- saturation deficit times wind movement, inches 
Hg km day -1 
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The means summarized at the bottom of the simple correlations table 
show that maximum and minimum daily temperatures in the dry plots were 
both warmer than in the wet plots and that the average daily temperature 
difference was 1.8 C. The daily excursion in temperature (M-M) was 
greater in the dry plots by 2.1 C, 


At the bottom of Table 5.1 information is presented on the order in 
which the independent variables were added in the step-wise multiple re- 
gressions, on the proportion of the sum of squares (R^) accounted for at 
each step where the new variable added had a regression coefficient that 
significantly deviated from zero (t-test) and on the standard error of 
estimate for each dependent temperature variable in this particular data 
set. 


For estimating Tmax, Tmin, and T5, average daily air temperature was 
always the first variable chosen. The second or third variable chosen 
was always either I/W (ratio of insolation for the current and previous 
day to wind run the current and previous day) or IOA (insolation outside 
the atmosphere or extraterrestrial radiation). The current day's insola- 
tion (INSOL) was the 2nd, 3rd, or 4th variable chosen in 5 out of the 6 
analyses. The first three variables accounted for from 89 to 93% of the 
variation in Tmax, Tmin, arid T5 for the dry plots and for from 88 to 91% 
for the wet plots. 

For estimating the max-mi n temperature at the 5-cm depth, INSOL was 
the first variable chosen, average of current and previous day's insola- 
tion ( ICDPD) was the second variable chosen, and TW was the third variable 
chosen. These three variables could account for 78% of the variation for 
the dry plots but for only 65% for the wet plot. Had we had a soil water 
content measurement daily it undoubtedly would have been included and 
would have probably been more important for the wet than for the dry soil 
plot. 

Still a different set of independent variables accounted best for 
the average daily soil temperature difference between the wet and dry 
plots (Twet-Tdry). The first variable chosen was I/W, the second was 
IOA, and the third was saturation deficit of the air times daily wind 
run, SDW, They accounted for 66 % of the variation in (Twet-Tdry). It 
is noteworthy that two of these involve windspeed. 


The standard errors of estimate of Tmax, Tmin, and T5 were all 
similar and ranged from 1.5 to 1.9 C. From air temperature alone (data 
not presented) Tmax, Tmin, and T5 could be estimated with standard 
errors 2.6, 2.0, and 1.6 C for the wet plots and 3.6, 2.1, and 2.1 C 
for the dry plots, respectively. 

In general these findings indicate that Tmax, Tmin, and T for 
shallow soil temperatures correlate well with average daily temperature, 
as in the literature cited. The variables ICDPD, INSOL, I/W, and IOA 
are all insolation measurements and they were always the 2nd, 3rd, and 
4th variables chosen to estimate soil temperature. Consequently, soil 
temperature estimates would be improved if ground level daily insolation 
data were routinely available from some source. Extraterrestrial radia- 
tion, IOA, is evidently a surrogate for daylength or season in the rela- 
tions, 
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The simple correlations between the independent variables used in 
the multiple regression analyses are summarized in Table 5.2. The in- 
solation variables derived from ground level measurements ICDPD, INSOL, 
and I£W are closely correlated and hence any one would suffice. Like- 
wise Ta and TW are highly correlated (r=0.832) so that TW is dominated 
by the information it_contains about temperature, not wind. IOA cor- 
relates highest with Ta (0.618) and mean dewpoint temperature, MDP 
(0.545), but contains significant additional seasonal or daylength 
information. Barometric pressure, BP, is negatively correlated with 
every other variable. It and SDEF were not included in any of the 
multiple regression equations. Saturation deficit times wind (SDW) 
which correlates most highly with variables containing ground level 
insolation and air temperature variables but contains some other 
information was the third variable chosen in estimating the wet minus 
dry soil temperature difference, Twet-Tdry. 

In summary, of 10 variables we used, air temperature, ground level 
insolation, and extraterrestrial radiation could be used to estimate the 
5-cm soil temperature parameters daily maximum (Tmax), daily minimum 
(Tmin) , and average daily temperature (T5) with average standard errors 
of 1.87, 1.77 and 1.53, respectively. Multiple regression equations 
containing these variables accounted for 88 to 93% of the variation in 
the 5-cm soil depth temperature measured with thermocouples for the 
period (preplant and planting) 01 Jan to 31 Mar 79. 
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INSULATED COPPER- 
NICKEL ALLOY WIRE- 

THERMOCOUPLE JUNCTION - 
INSULATED COPPER WIRE- 


20-GAGE, TYPE T THERMOCOUPLE 

WIRE, VINYL INSULATED 

AND JACKETED ^ 


THERMOCOUPLE JUNCTION- 
COPPER-NICKEL ALLOY WIRE, 

0.13MM DIA. 


TYPE T THERMOCOUPLE WIRE^ 
20-GAGE, VINYL INSULATED 


E \ 


COPPER WIRE, 0.13MM DIA; 


20-GAGE, ALLOY WIRE LOOP* — - 1 CM 

VINYL-INSULATED 

(B) 

Figure 5.1 Thermocouple configurations for 5-cm and deeper depths (A) and 
0.5-cm depth (B) . 
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SOIL association temp 





GROUPED SOIL ASSOCIATION TEMP* (C) 


26 FEBRUARY 1979 NIGHT 


£ SANDS AND DUNES 


5 SALINE CLAYS 


SHRINK-SWELL CLAYS 1 


I 2 LOAMS AND CLAYS 
, ! 3 LOAMY UPLANDS 


4 SANDY UPLANDS 


-40 


40 80 120 160 

DISTANCE INLAND (AM) 


Figure 5.3 Radiometric temperature for soil association categories 
grouped by texture. 
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Correlation coefficients between nine dependent soil temperature 
variables and ten environmental variables (upper part); and sum- 
mary of multiple regression variables chosen, R*, and standard 
errors of estimate of nine soil temperature variables (lower 
part). 


Dependent Variable 


, , , Wet Plots 

Independent ___ 

variable T T . T5 

max mxn 


Dry Plots 


T5 (M-M) T T . 
v ' max mxn 


(H-M) CT ttet -T dcy ) 


(1) T 

a 

(2) ICDPD 

(3) TW 

(4) INSOL 

(5) X/W 

(6) IOA 

(7) BP 

(8) HDP 

(9) SDEF 

(10) SDW 


,875 .920 ,945 .017 .847 .923 .934 .136 


.186 
.819 
.138 
.104 
.624 
-.642 
.800 
. 420 
.476 


.240 

.805 

.168 

.171 

.684 

-.654 

.782 

.438 

.480 


.505 .749 
.758 -.035 
.473 .850 
.448 .770 
.707 .142 
.667 -.106 
.649 -.270 
.633 .581 
.607 .402 


-.553 

-.653 

-.379 

-.631 

-.663 

-.625 

-.408 

-.209 

-.605 

-.507 


Mean: 
S.D. : 



. .2 8.95 

6.67 5.49 5.73 4. 21 1.31 

; r(.01)=0.283 * 

in step-wise multiple regression: 

1114 5 

4 6 2 3 6 

5 4 6 5 10 

6 - - - 


Coefficients of determination, R , for one, two, three,... variables with 
t-value significant at the 0.05 level. 


.766 

.846 

.893 

.524 

.717 

.852 

.872 

.722 

.440 

.867 

.872 

.902 

.613 

.904 

.872 

.912 

.765 

.615 

.876 

.884 

.908 

.648 

.920 

.893 

.927 

.785 

.659 

- 

.885 

- 

- 

.927 

- 

— 

- 

- 


Standard error of the estimate, C(f or maximum no. of terms in regression) 
1.90 1.72 1.50 1.98 1 1.84 1.81 1.56 1.99 .7 






Table 5.2 Simple correlation matrix among input environmental variables 
for estimating 5-cra soil temperature in wet and dry plots* 


T 

ICDPD 

TO 

INSOL 

I/W 

IOA 

BP 

MDP 

SDEF 

SDW 

3i 1 

(1) T " 

.339 

.832 

.310 

- r - 
.278 

.618 

-.730 

.780 

.569 

.582 

(2) ICDPD 

- 

.234 

.827 

.905 

.360 

-.242 

-.125 

.697 

.606 

(3) TO 


- 

.195 

.023 

.506 

-.622 

.725 

.417 

.664 

(4) INSOL 



- 

.772 

.295 

-.203 

-.156 

.717 

.617 

C50 I/W 




- 

.354 

-.174 

-.175 

.670 

.431 

(6) IOA 





- 

-.373 

.545 

.250 

.217 

C7) BP 






- 

-.571 

1 

• 

-o 

o 

h-» 

-.453 

(8) MDP 







- 

-.019 

.087 

(5) SDEF 








- 

.882 

10) SDW 









— 

Mean 16.2 

69.7 

56 S3. 

270.4 

• 78 

667.1 

29.9 

9.68 

0.19 

67.1 

S.D. 6.4 

124.3 

3045. 

149.9 

.45 

97.2 

.22 

7.62 

.12 

44.7 


n=80 r(.05)-0.217 ; r(.01)=0.283 
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6.0 CROP WATER STRESS 


6.1 Abstract 

HCMM temperatures for sorghum, citrus, and sugarcane crop sites 
were related to soil water depletion (DEPL) since the last rain or irri- 
gation calculated from the USDA-ARS Irrigation Scheduling Program (ISP), 
to days since irrigation or rainfall (DSR), and to the ratio of estimated 
actual to potential evapotranspi ration (ET/ETp) the day of and the day 
prior to HCMM overpasses. For a combined crops data set involving three 
HCMM overpasses, temperatures observed by HCMM (means of 3 x 3 pixel 
arrays, T3) were highly significantly correlated with ET/ETp (r=-.704**) 
and positively correlated with DSR (r=.484*). Temperatures were also 
highly significantly correlated with insolation on the day of the pass 
( IDOP) r=.812**, the vegetation greenness index (PVI6) calculated from 
Landsat for the corresponding dates (r=-.642**) and maximum air tempera- 
ture (Tmax) the day of the HCMM pass. The simple correlation between 
HCMM-observed surface temperatures and variables important in evapotran- 
spi rati on for seven dates (27 May through 08 Jul 78) consisting of 51 
observations were: insolation the hour of the HCMM overpass (IHROP), r= 

.668**), IDOP, r=.660**, ET/ETp, r=-.581**; DSR, r=-.477**, Tmax, r=.378** 
and, percent vegetation cover (PC), r=-.277**. For this data set, Tmax 
was highly significantly correlated with DSR, DEPL, and ET/ETp. The data 
demonstrate the strong correlation between variables related to evapo- 
transpi ration and HCMM surface temperatures, as well as the strong 
coupling of maximum daily air temperature to availability of water for 
evapotranspiration. In general, the findings indicate that surface tem- 
peratures observed by satellites, such as HCMM and the NOAA-6 to -9 series 
are responding to shifts in the energy balance of the Earth's surface that 
are strongly affected by vegetative cover and availability of water for 
evapotranspiration. Operational applications of such interpretations 
may be forthcoming under AgRISTARS. 

6.2 Introduction and objectives 

The remote sensing of plant canopy temperatures appears feasible in 
studying a broad array of situations: assessing the need for irrigation, 

conversely, the extent and severity of drought; distribution of rainfall 
as revealed by vegetative thermal patterns; proximity of water tables to 
the soil surface; pi ant- temperature variation with distance from coastal 
or inland water bodies; ecological studies of plant communities; crop 
inventories related to gross thermal patterns; and, estimating evapo- 
transpiration from farmland, forest, and rangeland (Myers, et al_, , 1970; 
Wiegand and Bartholic, 1970, and Wiegand, 1971). 

Plant canopy temperature is a function of the energy balance that 
is dominated by incident radiation (Wiegand and Namken, 1966) and by 
availability of water at the leaves for transpiration (Ehlig and Gardner, 
1963; Byrne, et aK , 1979). When water supply to the leaves is restricted 
by root zone water depletion, soil salinity (Thomas and Wiegand, 1970), 
or restricted rooting, transpiration is reduced and the energy balance 
shifts to a higher temperature equilibrium consistent with increased 
sensible heat loss. The energy balance is also affected by advected 
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energy, water vapor pressure gradient between the crop and the air 
(atmospheric demand) (Jackson et aU , 1981), by crop stomatal resistance 
(physiological), surface roughness (turbulent transfer), and albedo. 

Thus plant, soil, and ambient environmental variables all affect crop 

canopy temperatures. 

Crop canopy temperature is normally very close to ambient air tem- 
perature, a degree or two below at low light intensities and relative 
humidities (Pallas and Harris, 1964) to a few degrees above under severe 
water stress. The magnitude of effects are illustrated by the following 
examples from the literature: Cotton leaf temperatures were 3.5 C higher 

for plants with a leaf relative turgidity* of 59% (severely wilted) than 
for plants with a relative turgidity of 83% (well watered) (Wiegand and 
Namken, 1966). Tanner (1963) reported a 0 to 3 C foliage temperature 
difference between irrigated and nonirrigated potatoes; he estimated 
that a 10% decrease in transpiration by a full cover of alfalfa evapo- 
transpiring at 0.76 mm hr”' would cause a temperature increase of approx- 
imately 1 C when the associated windspeed at the 1-m height was 3.1 m 
sec”T and the net radiation was 0.95 ly min - *. Gates (1964) calculated 
the energy balance for single leaves and concluded that for each 0.10 
mm hr” T of transpiration, the temperature of the transpiring leaf would 
be 5, 2.5, and 1 C cooler than a similarly exposed nontranspiring leaf 
at windspeeds of 1, 5, and 15 mph, respectively. Wiegand and Namken 
(1966) who measured cotton leaf temperatures under partly cloudy condi- 
tions in early afternoon in June found a linear increase in leaf tem- 
perature of 9 C when insolation ranged from 0.6 ly min"' (cloud shadowed 
plants) to 1.6 ly min”' (full sun). 

In this task, we proposed to relate HCMM temperatures to soil water 
depletion and to days since irrigation or rainfall; and, to look at how 
canopy temperatures differed among the crops sorghum, citrus, and sugar- 
cane. 

6.3 Methods 


A computerized irrigation scheduling program developed by Jensen, 
Wright, and Pratt (1971), commonly referred to as the USDA-ARS Irrigation 
Scheduling Program, (ISP) (Wright and Jensen, 1978), provides estimates 
of the timing and amount of irrigation water needed using weather data, 
crop coefficients, and soil data. The necessary crop coefficients for 
local use of the method have been developed for cotton, sugarcane, and 
citrus (Salinas and Namken, 1977). Coefficients for grain sorghum were 
developed in this study. Given the date of the last irrigation or rain- 
fall , daily weather data, applicable crop coefficients, soil water holding 
characteristics, and rooting depth, the method calculates daily water 
depletion. At the time of the proposal submission, we were routinely 
scheduling the irrigations of about 150 fields and experimental plots 


* Relative turgidity (RT) = 100 (FW-DW)/TW-DW) , wherein FW is the fresh 
or field condition weight of leaf samples, TW is the turgid weight 
achieved by floating the leaf samples on distilled water overnight 
under illumination, and DW is the dry weight of the samples after dry- 
ing at 60 C, 






.1 ■ 
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of sugarcane, citrus, and cotton for our own research and farmer- 
cooperators. The method is used to schedule irrigations for about 
400,000 ha in the irrigated west (EPA, 1975). Thus we chose the cumu- 
lative water depletion and actual evapotranspi ration estimate on the 
day of and the day prior to the HCMM overpass dates calculated by this 
procedure as the acceptable standard to relate the HCMM indicated crop 
canopy temperatures to. 

Because of the nominal 40 ha instantaneous field of view of the 
HCMM data, we utilized various precautions in selecting fields for study. 
One precaution was to find fields surrounded by other fields of the same 
crop so that an ample area of crop, even if not all within the same field 
was obtained. Part of the cotton and sorghum fields were also chosen in 
nonirrigated areas where soil water conditions would be uniform among 
adjacent fields. These considerations as well as ground travel and man- 
power resulted in our selecting 22 fields that ranged from 18 to 1240 ha 
in size (only three <40 ha) that were devoted to sorghum, sugarcane, 
cotton, citrus, and improved pastures (grass). Ground truth acquired 
on weekly visits from 30 May to 13 Jul 78 included crop species, plant 
height (PH), percent crop cover (PC), percent weed cover, stage of 
maturity, and dates of irrigation (See section 1.2.2 for more detail). 

Ideally, the crops should have attained about 70% ground cover at 
maximum canopy development to eliminate exposed sunlit soil effects on 
the canopy temperatures and evapotranspiration rates (Tanner and Jury, 
1976). However, the crops did not all reach such ground covers, so we 
include percent cover (PC) and spectral vegetation indexes (Kauth and 
Thomas, 1976; Richardson and Wiegand, 1977) as variables in the analyses. 
The latter were available from two Landsat Satellite overpasses--for 15 
Jun and 12 Jul 78— that could be paired with HCMM data for 12 and 17 
Jun, and 8 Jul 78, respectively. (See section 1.2.4 for additional 
detail). Because ground cover ranged from 20 to >70%, we developed 
crop coefficient curves to cover the ranges, <40 PC, 40-60 PC, and >60 
PC. 

The equation ET = Kc ETp was used to estimate the soil water deple- 
tion where 

ET = estimated evapotranspiration (cm/day) for a given crop 
ETp = the evapotranspiration potential 
= 0.0112 (Tm-20)Rs 

Rs = incident solar radiation converted to equivalent depth 
of water evaporated (cm/day) 

Tm = mean daily temperature (°F) as obtained from the 
closest weather station 

Kc = Kco Ka + Ks 



where Kco is the crop coefficient 


Ka is a soil moisture availability factor 

Ks is. a factor that lets evaporation directly 
from the soil contribute to ET for a short 
time (3 days in summer) when the soil surface 
is wet from irrigation or rainfall. When Kco 
_<.9, Ks does not apply. 

The mean daily temperature from the weather station geographically 
closest to the chosen field was used to calculate ETp. Daily insolation 
as needed for Rs in the irrigation scheduling program (ISP) and hourly 
insolation for relating to the HCMM overpass temperatures were taken from 
the central weather station at Weslaco. (See section 1.2.5 for more 
details) . 

Rooting depth and plant available water in the root zone were es- 
timated for each crop on each soil type as shown in Table 6.1. ET was 
calculated for the dates 21 May through 15 Jul 78 inclusive. The book- 
keeping of soil water took rainfall, irrigations, and drainage in excess 
of field capacity into account. 

We acquired digital magnetic tapes for the reflective (0.55 to 1.1 
urn) and emissive (10.5 to 12.5 urn) bands of HCMM for afternoon overpasses 
on 27 May; 6, 12, 17 and 23 Jun; and 3 and 8 Jul 78. The overpass on 27 
May was 1416 1ST; the overpass time became slightly earlier on succeeding 
passes and was 1400 on 8 Jul. The scenes were subjected to a long series 
of processing steps including registration to the ground scene, assign- 
ment of latitude and longitude to each pixel, etc. (See report section 
2.1.1 for more detail). From these procedures it was possible to extract 
the single pixel that centered on the individual site of interest or to 
use the mean of 3x3, 5x5, 7x7 or 9x9 arrays of pixels centered on the 
central pixel as representative of the site. For the analyses of this 
task, we used only the central pixel (1*1 matrix of pixels) and the mean 
of the 3x3 array of pixels for analysis. We used the "original" HCMM 
calibration to calculate the HCMM temperatures fro'' ! ’he digital counts 
not the "recalibration" (Barnes and Price, 1980). i he original calibra- 
tion yields temperatures 5.5 C higher than the recalibration. 

We merged the data from HCMM (albedo and surface temperatures for 
1x1 (T-j , R-j) and 3x3 (T^, R^) pixel matrices) with ground truth (crop 
species; days since irrigation or rainfall (DSR); plant height (PH); 
percent cover (PC)), weather station (daily maximum (TMAX) and minimum 
temperatures (TMIN); rainfall), Landsat (vegetation indices KVI, BRT, 

PVI , PVI6), and central environmental station (insolation the day of the 
pass (IDOP) and insolation the hour of the pass (IHROP)) with outputs 
from the irrigation scheduling program (cumulative depletion (DEPL) 
since rainfall or irrigation), ratio of actual evapotranspi ration to 
potential evapotranspi ration for the day of and the day preceding the 
HCMM overpass (ET/ETp), and calculated site distance from the coast 
(OIST) into three data set combinations: 
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(1) HCMM data for 3 dates paired with Landsat data for 2 dates 
plus the ground truth, weather station, irrigation scheduling 
and central environmental station data for the appropriate 
dates. This set was available with cloud and water vapor 
contaminated pixels included and with the affected data 

dele ted. 

(2) HCMM data for 7 dates (27 May through 8 Jul) paired with ground 
truth, weather station, irrigation scheduling, and central 
weather station data for the same period (15 Aug 78 data (see 
Table 1.2) not included for lack of ground truth and irrigation 
scheduling data past 15 Jul). 

(3) T-j , R] , and T 3 , R 3 segregated by the crops sugarcane, citrus, 
and sorghum for 7 dates. (Did not include cotton or grass 
because of insufficient number of observations). 

The raw and derived data for the 22 intensive sites are summarized in 
Appendix A. The cloud and vapor contaminated data have been deleted. 

Our approach in analyzing the data was to use the most complete data 
set we could initially, then delete sources of variation such as cloud 
and vapor contaminated pixels and then highly correlated variables. We 
used a statistical program W1106 included by IBM with its basic statis- 
tical package for our IBM 1800 computer. This step-wise multiple re- 
gression program calculates the simple correlation among all the variables 
in the data set before performing a step-wise multiple regression analysis 

6.4 Results 

We were extremely fortunate that 7 usable overpasses of HCMM occurred 
while we were ground truthing and obtaining the data for the irrigation 
scheduling computations, and that Landsat data were available for two 
dates to augment the ground truth taken. 

Results from the various data sets may be summarized as follows: 

6.4.1 Correlation of crop water stress parameters and 

vegetation indexes with HCMM temperatures and albedos 

6.4.1 .1 Before and after removal of atmospherically 
contaminated data 

Data set 1 consisting of HCMM data for 3 dates paired with corre- 
sponding LANDSAT, ground truth, weather station, irrigation scheduling, 
and central weather station data were used for these analyses. 

The simple correlation of the albedo (R) and surface temperature 
(T) and the other variables are summarized in Table 6.2. Results are 
presented for the central pixel for each site (field) and a 3x3 matrix 
of pixels centered on the central pixel. In the upper portion of Table 
6.2 results are presented for the data that included cloud and vapor 
contaminated pixels (n = number of observations = 46) and in the lower 
part with the cloud contaminated pixels and some of the variables deleted 
(n=22). 
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Conclusions reached from considering the content of Table 6.2 are: 

(a) The simple correlations among variables is higher for the 
sample of 9 pixels (R,» T,) than it is for the central pixel 
(R], T-i) in both the upper and lower parts of Table 6.2. 

(b) The simple correlation of other variables with temperature are 
improved by deleting the cloud and vapor contaminated pixels, 
but the correlations of albedo are weakened. (This occurs for 
albedo even though the standard deviation about the mean is 
lessened whereas it remains about the same for temperature.) 

(c) Albedo is significantly correlated with water depletion (DEPL) 
and insolation on the day of the pass (IDOP) before contaminated 
pixels are removed but only with IDOP after removal of con- 
taminated data. 

(d) Before deletion of contaminated data, T] and T 3 are most 
closely correlated with distance from the coast ( D 1ST) but 
they are also significantly correlated with the ratio of 
actual to potential evapotranspi ration the day of and the 
day prior to the pass (ET/ETp), maximum air temperature the 
day of the pass (Tmax) and the vegetation indices PVI, PVI 6 , 
and GVI. After deletion of the variable DIST (to be dis- 
cussed under synoptic analyses), all but PVI 6 of the vegeta- 
tion indices (because they are highly intercorrelated and PVI 
is least affected by atmospheric conditions) and the contami- 
nated data, HCMM temperatures are significantly correlated 
with days since rainfall (or irrigation) (DSR), ET/ETp, Tmax, 
IDOP, PVI 6 , and albedo. 

The multiple regression analyses are useful for indicating the 
amount of variation explained by the environmental variables. Table 
6.3 lists the variables in the order they were included in the step- 
wise regression, gives the multiple correlation coefficients, and 
provides the standard error of estimate of 3x3 pixel temperature (T 3 ) 
for the first eight variables included in the analyses for data set 1 
before and after removing cloud contaminated pixels and DIST. 

Table 6.3 shows that eight variables explained 66 % of the variation 
in T 3 before removal of cloud and vapor contaminated pixels and 96% of 
the variation after their removal. These results point out the over- 
whelming variability imposed on the HCMM thermal data by atmospheric 
effects. The standard error of estimate of mean temperature, ST 3 , de- 
creased as additional variables were added to the regression analysis. 

It was 1.1 C for the noncontaminated data and 2.8 C for the contaminated 
data when eight variables were included in the respective regression 
equations. 

The corresponding R^ and ST for the central pixel (data not pre- 
sented) for eight variables in the analyses were .601 and 3.3 before 
contaminated data were deleted and .941 and 1.4 after they were deleted. 
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Generally these results show the strong effect of an unclean at- 
mosphere on the observations, a strong correlation between variables 
related to evapotranspiration and temperatures observed by HCMM, and 
that well over SOS of the variation in observed temperature is ex- 
plained by the variables included in the analyses when samples from 
cloud or vapor contaminated atmosphere are deleted. 

6.4.1 .2 Within the crops sugarcane, citrus, and sorghum 


Data set 3 consisting of and Rg data for 7 dates segregated by 
crops were used for this analysis. 

Table 6.4, which summarizes the statistically significant simple 
correlation coefficients between each pair of variables by crops for 7 
HCMM overpass dates, is an attempt to get a further feel for the vari- 
ables that were consistently related. In interpreting the data, it is 
worthwhile to keep in mind that the sorghum fields were mainly from non- 
irrigated areas, the fields were large and most of the fields had been 
harvested by 08 Jul . The sugarcane fields were making rapid growth and 
both plant height and ground cover were increasing rapidly; to maintain, 
the growth they had to be irrigated regularly after the rainfall in 
early June. The citrus orchards were the smallest sites used, the age 
and condition of groves within the test sites varied considerably, and 
citrus is moderately drought tolerant. 

From Table 6.4 we see that To for the sorghum sites is significantly 
correlated with 8 other variables (PC, DEPL, ET/ETp, IDOP, IHROP, DSR, 
Tmax, and PH) whereas it is significantly correlated with only two for 
sugarcane (IDOP, IHROP) and one for citrus (Tmax). Plant height (PH) is 
significantly correlated for all 3 crops with PC, DEPL, and Tmax. 


Regarding variables that affect water use, we see that calculated 
depletion (DEPL) is significantly related to PC as in the real world; 
the two could be significantly related because the crop coefficients 
used to estimate evapotranpsi ration are ground cover dependent. But 
crop coefficients used here were either taken from the literature or 
were derived from independent data. Days since rainfall or irrigation 
(DSR) and ET/ETp, two variables that are not methodologically dependent 
are very highly correlated for each of the crops (SC=~.945**, CT=-.778**, 
and SR=-.910**). Such findings support the validity of the evapotran- 
spiration estimation procedures and their utility for irrigation sched- 
uling. 


Figure 6.1 presents the cumulative percent of the variation in T -j 
and T^ as variables were added to the multiple regression for lxl and 
3x3 pixel samples representing sugarcane, citrus and sorghum fields. 

For every crop by the time 6 variables had been included in the regres- 
sion equations for the 3x3 pixel sample, >95% of the variation in T3 
was explained. For the lxl pixel sample, from 85% (sugarcane) to 92% 
(sorghum) of the variation was explained by 6 variables. For citrus and 
sorghum 3 variables explain over 80% of the variation in T-| and T-j, but 
not for sugarcane. 


3 * 
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From a look at Figure 6.1, it appears that 6 variables are desirable 
to describe the behavior of T 3 . Numerical values of R 2 and ST 3 for 3 and 
6 variables, and the variables included were (where variables with sig- 
nificant regression coefficients are underscored) : 

6 . -.2 Correlation of PC, DEPL, I00P, IHROP, TMAX, ET/ETP, 
and PH with HCMM temperatures and albedos for seven 
observation dates 

Data set 2 consisting of HCMM data for 7 dates (27 May through 8 Jul 
78) during which ground truth and evapotranspi ration calculations were 
being made were used for the analyses. The simple correlation coeffi- 
cients for 8 variables with the central pixel albedo and temperature 
(R]» Ti) and with the 3x3 pixel mean of albedo and temperature (R 3 , T 3 ) 
for 7 HCMM overpass dates are summarized in Table 6.5. The same variables 
that correlated with R and T for the three dates of data (Table 6.2) again 
correlate. However, the simple correlation coefficients are larger; they 
are, therefore, significant at a higher probability level because the 
number of observations is greater. Also as in Table 6.2, the simple 
correlations are slightly higher for the 3x3 pixel mean data (R 3 , T 3 ) 
than for the central pixel (Ri, T-|) data. The only exception is the 
correlation with percent vegetative cover (PC). 

The simple correlation between albedo and temperature is highly 
significant, whereas it was not for the three dates of Table 6.2. This 
seems to be related to the facts that the additional dates in the data 
set were associated with a greater range in the variables such as in- 
solation the day of (IDOP) and the hour of the HCMM pass (IHROP), and 
in other influencing variables such as percent vegetative cover (PC) 
and plant height (PH). IHROP was not significantly correlated with any 
other variable in the data set in Table 6.2 because those dates had 
quite similar insolation the hour of the HCMM overpass. 

The finding that insolation (IDOP, IHROP) strongly affects surface 
temperature is as it should be because the sun's energy is the "driver" 
of the earth's energy balance. The highly significant negative corre- 
lation between the ratio of actual to potential evaporation (ET/ETp) 
for the day of and day prior to the pass indicates that as water becomes 
less available for transpiration, the canopy temperature increases. 

This finding, again, agrees with reality and the physical and physio- 
logical processes involved. The significant positive correlation between 
maximum air temperature (Tmax) the day of the HCMM overpass and the site 
surface temperature observed by HCMM is also in agreement with reality; 
that is, the two are coupled. 

In the energy balance equation, 

Rn = G + H 4- LE, 

if incident solar radiation remains constant for a series of summer days 
and heat flux into the ground G is a small, essentially fixed proportion 
of net radiation Rn for daily periods, it is obvious that a decrease in 
water availability which reduces the latent heat of evapotranspi ration 
LE has to be compensated for by an increase in sensible heat flux from 
the canopy to the air, H. 
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Jackson et al_. (1981) and Monteith (1973) further show that H 
and LE, respectively, can be expressed by 

H = pep (Tc-TA)/ra 


arc LE = pep (e* c -eA)/[>r(ra+rc)] 


where p is the density of the air, cp the heat capacity of air, Tc 
canopy temperature, TA air temperature, e* the saturated vapor 
pressure at Tc, eA the vapor pressure of tne air, 2r"the psychrometric 
constant, ra the aerodynamic resistance of air to vapor transport, 
and rc is the canopy resistance to water vapor transfer to the air. 


In our studies ET = LE became smaller the longer it had been since 
rain had fallen or irrigation water had been applied. Consequently, 
(Tc-TA) which is directly proportional to H had to increase. If ground 
cover was essentially complete, THCMM was sensing Tc. From May through 
July when the data for these studies were being obtained, TA was also 
increasing seasonally. 

Jackson et al_. (1981) also defined a crop water stress index 
(CHS I) as CSWI = 1-ET/ETP. In Appendix Table A. 2 we list ET/ETP for 
7 sorghum, 6 sugarcane, 4 cotton, 3 citrus, and 2 grass fields for 
each of 3 dates. The very low ET/ETP for a number of the dryland 
sites on JD 189 means they were under severe stress. 

Jackson et aj_. (1981) also document that there is a linear relation 
between Tc-TA and the vapor pressure deficit of the air, e*^-e A. When 
we performed the step-wise multiple regression analyses for several dates 
(see Figure 2.5) we found that the saturation deficit of the air (SDEF) 
was the first variable selected and that it alone explained over 50% 
of the variation in THCMM. Consequently, we believe interpretation of 
canopy or "Earth surface" temperatures between geographical locations 
that differ widely in water vapor pressure’of the air will have to 
include the effect of differing saturation deficits of the air in in- 
terpreting Tc and Tc-TA observations. 


In Table 6.6 the simple correlation matrix among the "independent" 
variables of Table 6.5 are presented. Tmax is highly significantly 
positively correlated with millimeters water depletion (DEPL) (.486**) 
and days since rainfall or irrigation (DSR) (.567**) and highly signifi- 
cantly negatively correlated with ET/ETp (-.474**). These correlations 
indicate clearly that availability of water for evapotranspiration is 
affecting partitioning in the energy balance. Consequently, a decrease 
in availability of water for evapotranspiration is associated with an 
increase in the daily maximum air temperature. The actual to potential 
evapotranspiration ratio (ET/ETp) is significantly correlated with 
vegetative cover (PC), and days since rainfall or irrigation (DSR) as 
well as Tmax. DSR and DEPL are calculated for the same time interval 
and are highly correlated. Plant height (PH) is weakly correlated with 
only one variable, PC. PH is not a sensitive variable when the crops 
are pooled because only cotton and sugarcane were changing in plant 
height; sorghum had headed by 27 May, citrus is a perennial whose height 
remained constant for the two month period, and grass height changed 


little during the observation period. 
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The multiple regression equations for estimating the 3x3 pixel 
HCMM temperatures (through the first 4 variables), and the R 2 and the 
standard error of the estimate, S? 3 (through 5 variables) were: 

No. of 




d2 

st 3 

Variables 




(°C) 


T 3 = 24.67 

+ .090 (IHROP) 

.446** 

2.9 

1 

= 29.61 

+ .079 (IHROP) - 6.863 (ET/ETp) 

.672** 

2.3 

2 

= 23.59 

+ .070 (IHROP) - 6.830 (ET/ETp) 

.701** 

2.2 

3 


+ .386 (R 3 ) 




= - 26.! 

93 - .162 (IHROP) - 7.273 (ET/ETp) 

.781** 

1.9 

4 


+ .989 (R 3 ) - .162 ( ID0P) 






.803** 

1.8 

5 


A t-test of the regression coefficients was significant for each of the 
first five variables (R 3 , ET/ETp, I00P, IHROP, Tmax) included in the 
multiple regression equations; they accounted for 80% of the variation 
in T 3 . The standard error of the estimated temperature decreases from 
2.9 C for inclusion of the first variable to 1.9 C by the time the 
fourth variable is included. The corresponding multiple correlation 
an_d STx for the central pixel data were: for 4 variables, R 2 = .722, 

ST -j = 2.3 C and for 8 variables R 2 = .764, and ST^ = 2.2 C. 
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VARIABLE NUMBER 


Figure 6.1 Percent of variation explained as variables were added to the 
multiple regression for lxl and 3x3 pixel samples of sugarcane 
citrus v and sorghum. 
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Table, 6.1. Soil types, crops, rooting depth, and water characteristics of fields 
and sites for which evapotranspiration was deteroined. 


Soil 

type 


Swelling 
clay ] 


Loamy & 
clayey 
(clay 
) loam) 


Loamy 

upland 

(sandy 

clay 

loam) 


Fine j 
sandy j 
loam • 


Crop^ 


! Water holding j Plan 
■Rooting capacity of availa 
: depth soil ; vate 


cm >in.H 0/ cm.H 2 0/ * 

ft. soil cm. soil . cm 


Sorghum / 
Sugarcane / 
Cotton J 
Citrus 


90 | 2.34 
90 2.34 
90 2.34 


Sorghum / j 120 

Sugarcane r j 120 

Cotton / i 120 


Citrus 

Grass 

Sorghum 

Sugarcane 

Cotton 

Citrus 

Grass 


/ i 120 

i 

j 90 

j 120 

i 

j 120 

/ ■ 120 

! 

! 150 

l ■■ ■ 

/ j 90 

/ ; 120 


Sorghum / 150 • 

Sugarcane / ' 120 

■■ 1 ; 

Cotton J ! 150 

’ 

Citrus 120 

Grass / 150 


.195 I 17.8 




Table 6.2. Linear correlation coefficients for all variables with the central pixel 

albedo and temperature (R^ and with the 3 x 3 pixel mean value of 

albedo ar.d temperature (R^, T^) for 22 intensive sites. Top portion of 

table is for cloud and vapor contaminated pixels in the analysis, and 

a/ 


Variable 

Acronym 


the lower part is with them deleted.— 


I Mean + SD 


PC(%) 

PH (cm) 
DSR(days) 
DEPL. (mm) 
ET/ETp 
Tmax(°C) 
IDOP(Ly) 
DIST (Km) 
PVI 

PVI6 • 
GVI 
BRT 
R(%) 

T(°C) 


PC 

DSR 

DEPL 

ET/ETp 

Tmax 

IDOP 

PVI6 

R 

T 


50.0 + 
136. + 

14.6 + 
59 + 

.56 + 

36.1 + 
415. + 

79.4 + 

14.2 + 
19.8 + 

23.4 + 

76.3 + 

17.6 + 

29.7 + 


22.4 

116 

10.6 

37.6 

.29 

I. 5 
50 
22.4 

7.1 

10.1 

II. 3 
9.3 
2.1 
4.8 


49.2 + 23.2 

14.9 + 9.9 
55.7 + 30.5 

.51 + .27 

36.0 + 1.6 
408 + 54.4 

19.0 + 9.9 

17.2 + 1.6 

31.9 + 4.7 


R 1 

T 1 

R 3 : 

• T 3 1 

Correlation coeff . , r 

.095 

-.156 

.145 

-.203 

-.015 

.145 

.068 

.125 

-.276 

.266 

** 

-.388 

.286 

* 

-.299 

-.116 

* 

-.359 

-.123 

.102 

. ** 
-.441 

.275 

aa 

-.505 

-.151 

* 

.350 

-.258 

** 

.375 

* 

.367 

.278 

** 

.454 

* 

.317 

.221 

A A 

.696 

.107 

** 

.761 


»u 


JL.ru 

-.058 

-.343 

.113 

-.387'' 

-.030 

** 

-.404 

.152 

-.450 

-.039 

j.jt . 

O A 

-.374 

.139 

a* 

-.419 

-.037 

.217 

; -.086 

-.221 

- 

.080 

- 

-.017 

.080 

“ 

i -.017 

1 

— 

.107 

-.255 

.270 

-.269 

-.030 

.448 

-.173 

.484* 

-.220 

-.126 

-.242 

(N 

CO 

O 

• 

I 

-.190 

** 

-.681 

-.069 

-.704** 


** 


** 

.081 

.583 

.043 

.619 

.614** 

A* 

.773 

aa 

.679 

kk 

.812 

-.206 

aa 

-.635 

-.106 

aa 

-.642 

- 

.537" 

- 

■ A* 

.492 

JL 

.SB?” 

- 

* 

.492 



Mean + SD 


Same 

as 

for 

central 

pixel 


* 

17.6 + 1.6 

29.7 + 4.4 


Same 
as . 
for 

central 

pixel 


17.2 +1.3 
31.9 + 4.5' 


a/ 

— For upper portion of table, number of observations of each variable, n 
whereas in the lower part n = 22. 

_ 128 


= 46 




Table 6.3. Effect of cloud and vapor contamination on proportion of the 

sum of squares explained by regression (R“) and the standard 

error of estimate (S= ) of mean temperature (T_) . 

3 


v 

Contaminated pixels included 
(tr« 46) 

Contaminated pixels 
(n = 22) 

deleted 

. 

Variable 

No. 

f 

Acronym 

** 

S T 
1 3 

(°c) 

■ ' 

Acronym 

R 2 

% 

°C 

1 

1 

ET/ETp 

.255 

. 

3.9 

IDOP 

.659 

2.7 

2 

DEPL 

.332 

3.7 

PVI6 

.824 

2.0 

; 3 

! 

Tmax 

.457 

3.4 

ET/ETp 

.884 

1.6 

1 4 

l 

DSR 

.496 

3.3 

R 3 

.887 

1.7 

i 5 

! BRT 

.587 

3.2 

Tmax 

.895 

1.7 

\ 6 

PH 

.587 

3.1 

DEPL 

.923 

1.5 

7 

PVI 6 

.635 

2.9 

DSR 

.944 

1.3 

8 

GVI 

.661 

2.8 

PC 

.960 

1.1 


t-test of regression coefficients! 
significant for each of first 3 
variables and for DEPL, PH, and 
BRT among 8 terms. 


t-test of regression coefficients 
significant for each of first 3 
variables and for Tmax, DEPL, 

DSR, R 3 , ET/ETp and PC for 8 
terms . 





fable 6. 4 Summary of statistically significant linear correlation coefficients 
between each variable measured for sugarcane (SC, n»12), citrus (CT, 
n«ll), and sorghum (SR, n<*18) for seven overpass dates from 27 May 
to 08 July, 1978. 



* designates statistical significance at the 0.05 level 

* designates statistical significance at the 0.01 level 









































7.0 PLANT COVER 


7.1 Abstract 

HCMM temperatures 9 for five crops (sorghum, citrus, sugarcane, 
cotton and grass os s cures) as one data set were significantly linearly 
correlated with the vegetation indexes GVI, PVI, and PVI6 but not with 
percent cover (PC) or plant height (PH). For sorghum, PC (as well as 
the vegetation indexes) was significantly correlated with HCMM tempera- 
tures. Weather (insolation, saturation deficit of the air, and maximum 
daily air temperature on HCMM overpass dates), geographic (distance of 
sample sites from the coast) and time (days since irrigation or rainfall) 
variables influence^ HCMM temperatures as much as or more than vegetative 
conditions when crops that differed in stage of maturity, planting con- 
figuration, and PC vs PH relation were combined in one data set. For 
cotton and sugarcane which were increasing in ground cover during the 
observation period (27 May to 15 Aug 78), the cooling effect of in- 
creasing ground cover was about equally offset by the seasonal increase 
in ambient air and agricultural landscape temperature. In general, the 
spread in HCMM temperatures among crops was narrow and the temperatures 
were below ground-observed air temperatures at the time of the HCMM data 
acquisitions, 1400 to 1416 hours 1ST. 

From these studies we anticipate that surface temperatures observed 
by the polar orbiting operational meteorological satellites (NOAA-S to 
-9 series) will correlate with vegetation indexes calculated from their 
visible and near-infrared bands and be interpretable in the AgRI STARS 
program for synoptic seasonal crop development and drought occurrence 
appl ica tons. 

7.2 Introduction and objectives 

We anticipated that the surface temperature of bare fields would be 
approximately 20 C higher than those with a dense crop canopy (Wiegand, 
et al_. 1968) and that there would be an approximately linear relation be- 
tween plant cover, LAI or spectral vegetation indexes (Kauth and Thomas 
1976; Richardson and Wiegand 1977) and temperature. 

Success in identifying fields that differ in amount of ground cover 
by living or dead vegetation would have many applications in protecting 
the soil against water and wind erosion. 

Also, the thermal data could be used for the Synoptic monitoring of 
drought. For this appl i cation, the more severe the drought the higher 
the observed surface temperatures should be; that is, the energy balance 
(Byrne et al . 1979; Soer 1980) would shift to more of the energy going 
into sensible heat associated with increased soil heating and relatively 
less into the latent heat flux associated with evapotranspiration. 

Drought would also limit plant growth, hence ground cover, leaf area 
indexes, or crop canopy greenness. Al though we coul d not anticipate a 
drought, we selected fields in both irrigated and dryland fanning areas. 


a The HCMM temperatures reported were not adjusted for atmospheric 
attenuation. 



This comparison simulated a drought vs nondrought situation in that 
the annual free water evaporation is approximately 180-cm for the study 
area compared with the normal 60- cm of rainfall, so that dry cropland 
and rangeland sites experience plant growth limiting water deficits 
periodically throughout the year. 

In this study we wanted to document the effect of crop cover on 
HCMM surface temperature. The extent to which this objective could be 
met using HCMM data depended upon selection of adequately sized homo- 
genous areas, degree of atmospheric attenuation of emitted energy, 
registration of HCMM data to the ground scene and other considerations 
that impacted the quality of data available and on the ability to 
implement the analysis, and not on the technical soundness of the 
objectives per se in ground level or low altitude aircraft investi- 
gations. ' i; i {, v . 

For bare soil, Idso et a]_. (1975) have shown that the difference 
between daytime maximum and nightime minimum surface temperature (ATs) 
was inversely related to soil water content in the surface 2 to 4-cm 
of soil. Cihlar _et al_. (1979) expressed water content of fallow fields 
as percent of field capacity ( PFC) and confirmed the inverse linear Ts 
versus PFC for the 0 to 2- and 0 to 4-cm depths. 

The relationships are more complicated when plant canopies and 
crop residues partially cover the ground, but shallow soil moisture 
apparently continues to dominate. For example, Cihlar et al_. (1979) 
who measured percent bare ground, plant height, percent plant cover, 
and percent straw cover in cropped and post harvest barley fields were 
able to explain from 5 to 15% more of the variation in ATs by including 
percent straw cover and percent bare ground in the regression analyses. 
Wiegand e_t aJL (1968) found that the equivalent blackbody temperature 
from aircraft scanner date of grain sorghum in small thinned plots at 
1400 hours LST was a linear inverse function of leaf area index (LAI); 
the temperature difference between LAI 1 s of 0.8 and 4.2 was 17 C, 
Heilman and Moore (1980) found that 1330 hour soil temperature, beneath 
a barley canopy that changed from 30 to 90% ground cover (PC) over the 
45 day study period, as measured by thermocouples 1-mm below the soil 
surface could be predicted (r =.86, S.D.=2.6 C) from canopy temperature 
measured by a Barnes infrared thermometer and an exponential function 
of PC. There was no improvement in estimating soil temperature by 
including LAI, insolation, or maximum daily air temperature. Cihlar 
(1980) studied the day minus night temperature di fference, ATs, as a 
function of plant available water in the 0 to 2 and the 0 to 20-cm 
surface soil layer for barley canopies that ranged from 60 to 98% 
ground cover. The correlation between the two decreased as plant 
cover decreased, and he concluded that plant water stress would be 
accurately indicated by infrared thermometry only for sufficiently 
dense canopies. Sutherland and BartholJc (1977) considered the 
emissivity appropriate for correcting aircraft scanner apparent tem- 
peratures to true temperatures and concluded that if the crop height- 
to-spacing ratio is greater than 1, the composite emissivity of 
agricultural areas is between 0.98 and 1.00 and temperature errors 


are <2 C. Byrne et aU (1979) urged that canopy temperature data be 
suppTemented by data from other sources that document atmospheric 
demand, root density, and plant and soil characteristics since leaf 
water status— hence canopy temperature— is affected by these as well 
as mean moisture content in the root zone. 


7.3 Methods 


The data for this task are described in report sections 1.2.2 
through 1.2.4, and they were processed as described in report sections 
2.1.2, 2.1.3, and 2.1. 5. Other information on methods used that are 
necessary for understanding the data and how they were used are given 
in the particular results section where it applies. The variables 
used are of general information to all results subsections, however, 
and are dev i tied as follows: 


R3 - average albedo relative to the solar constant at the top of 
the atmosphere from HCMM reflective band (0.55 to 1.1 urn) 

for 3x3 arrays of ground site pixels {%). The digital counts 

given in the HCMM data were ratioed to 255, the digital count 
for a perfect reflector normal to the sun's rays, and multi- 
plied by 100 to convert to percentage. 

T3 - average temperature from HCMM thermal band (10.5 to 12.5 um) 

for 3x3 arrays of ground site pixels (°C). The equation used 

that applies to the temperature range 260 K (0 digital counts, 
DC) to 340 K (255 digital counts) was 



1251.6 
14421.6 
n \DC + 118.21 


+ 1 . 0 ) 


PC - plant cover estimated visually from a nadir view as fraction 
of ground area obscured by crop plants and weeds {%) . It is 
an integrated value for the field that incorporates skips in 
rows into the estimate. V.V'-.-.'l 

PH - plant height as measured from the soil surface to the upper- 
most plant part as determined by sighting across the tops of 
plants (cm). 

PVI - perpendicular vegetation index calculated from Landsat bands 
5 (0.6 to 0.7 um) and 7 (0.8 to 1.1 um) digital counts (unit- 
less). 

specifically, 

PVI = ( (Rgg 5 - Rp5) 2 + (Rgg7 - Rp7) 2 )^ 

where Rp5 and Rp7 are the digital counts for the 
candidate pixels in Landsat data and Rgg5 = 0.851 
Rp5 + 0.355 Rp7 and Rgg7 = 0.355 Rp5 + 0.143 Rp7 
are the points of intersection on the soil back- 
ground line of perpendiculars from Rp5 and Rp7 
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PVI6 - perpendicular vegetation index calculated from digital counts 
of Landsat bands 5 and 6 (0.7 to 0.8 urn) (unitless). 

GV I - Kauth's green vegetation index = -0.290 (4) - 0.562 (5) + 

0.600 (5) + 0.491 (7) where 4, 5, 6, and 7 denote digital 

count* ■ c - ? Landsat multi spectral scanner (MSS) bands 4, 5, 

6, and 7, respectively (unitless). 

BRT - Kauth's brightness index * 0.433 (4) + 0.632 (5) + 0.586 (6) + 
0.264 (7) wherein; 4, 5, 6, and 7 denote the digital count of 
Landsat MSS bands 4, 5, 6, and 7, respectively (unitless). 

IDOP - insolation on the day of the pass up to the hour of the over- 
pass (langleys = cal cm" 2 ). 

IHROP - insolation the hour of the overpass, measured from 1300 to 
1400 hour (langleys). 

SDEF saturation deficit of the air for the hour of the overpass (mb) 

Note: IDOP, IHROP and SDEF were measured at one site, the 

central weather station at Weslaco, TX. 

Tmax - air temperature maximum the day of the pass from the observation 
station, of 22 in the study area, nearest the ground site of 
interest (°C) . 

DSR - days since rainfall (or irrigation) from rainfall reporting 
station nearest the ground site of interest and ground obser- 
vations of irrigations in process (days). 

For readers who may not be familiar with the vegetation indexes, 
experience indicate that PVI, PVI6, and GVI are most responsive to the 
amount of living (green) vegetation present whereas BRT responds to 
the amount of exposed soil in the sensor field of view. As mentioned 
in Chapter 5, percent cover (PC) may be living or nonliving vegetation 
(plant residue or senescent plant parts). When the vegetation present 
is healthy and green, PC and the greenness indicating vegetation indexes 
are correlated; however, if the vegetation present is a mixture of living 
and nonliving plant parts, PC may not be correlated with the vegetation 
indexes or BRT. 

The albedo (or reflectance) from the 0.55 to 1.1 um band of HCMM 
straddles the visible wavelengths from 0.55 to 0.70 um wherein soils 
are more reflective then plants and the 0.70 to 1 .05 um range where 
plants are more reflective than soil. Consequently, albedos over this 
wavelength interval are nondescript and may not relate to any of the 
vegetation indicators, t.t.'.. 
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7.4 Results 


7.4.1 Vegetation indexes, percent cover (PC) and plant height 
(PH) in relation to HCMM temperatures and albedos 

Wa used the data for the 22 intensive sites where PC and PH had 
been measured weekly from 30 Hay to 13 Jul. These data were supple- 
mented with Landsat data for two dates (15 dun and 12 Jul 73) that 
could be paired with HCMM albedo and temperature data for three dates 
(12 and 17 Jun and 08 Jul 78). Because clouds obscured some of the 
fields on Landsat overpass dates and some of the HCMM data were deleted 
due to cloud and vapor contamination, we were left with 32 data points 
from the three dates for the intensive sites. There were 30 Landsat 
vegetation index data points from the representative and weather sta- 
tion sites, but ground observations of PC and PH were not taken. These 
data sets were analyzed separately and then were combined into one set. 
The linear correlation coefficient matrices for the HCMM albedo (R3) 
and temperature (T3) with the vegetation data are presented in Table 7.1. 

The linear correlation coefficient matrices of Table 7.1 show that 
HCMM-observed temperatures correlated significantly and positively with 
albedo and highly significantly and negatively with PVI, PVI6 and GVI 
(part C of Table 7.1) the same as for the intensive sites only (part A 
of Table 7.1). Thus, temperature decreased as amount of green vegetation 
increased. The positive linear correlation between R3 and T3 indicated 
that exposed soil dominated the albedo, probably because the sun's energy 
is concentrated in the visible wavelengths where soil is more reflective 
than vegetation. 

The data for the intensive sites where PC and PH were observed 
(part A of Table 7.1) indicated that percent cover and plant height were 
linearly correlated (r=0.422*). However, the correlation was weak due 
to the variety of crops and planting configurations involved. As shown 
at the bottom of Part A, the 22 observations were composed of 9 sorghum, 

5 sugarcane, 4 citrus, 1 cotton, and 3 grass pasture sites. The noncor- 
relation between PH and T3 is probably due to the same reasons. For 
example, the citrus averaged 370 cm in height as opposed to Bermudagrass 
sod that was only 10 cm tall. The citrus trees are typically planted on 

6 or 7 in spacing, and in spite of their height, they averaged only 533 
ground cover. The row spacing of sugarcane is typically 1.5 m compared 
with 90-cm for sorghum and cotton. Sugarcane was growing rapidly during 
the experimental period; it averaged 146 cm in height on 27 May and 201 
cm in height on 08 Jul. Sorghum was at the grain filling stage by the 
time the measurements were begun so its maximum height and greenness 
had already occurred; it was, in fact, senescing (lower leaves drying 
up) during the observations period and two fields had been harvested by 
08 Jul. There was little rainfall during the growing season except the 
first week in June; consequently, dryland cotton and sorghum were water 
stressed on 27 May and after mid- June. 
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The vegetation index GVI is very highly correlated (r= 0.97** to 
0.99**) with the vegetation indexes PVI and PVI6. Whereas GVI is de- 
rived by rotation and translation of the coordinate axes in Landsat 
data space so that one of the axes contains the maximum amount of 
information about green or living vegetation, PVI and PVI 6 are based 
on the perpendicular distance from the: soil background line in Landsat 
data for pixels containing vegetation. For PVI and PVI6 the greater 
the displacement from the soil background line (PVI and PVI 6=0 for the 
soil background line) the greater the vegetation density or green leaf 
area index (LAI) (Wiegand et a]_. 1979). The high interrelation among 
these independently deriveH""vegetation indexes adds confidence that 
they describe vegetation conditions. The PVI6 is less affected by 
atmospheric conditions than PVI; consequently it gives higher linear 

correlation coefficients with R3 and T3 than does PVI.! 

■■ ' . 1 5 

The Kauth brightness (BRT) is negatively correlated with the green 
vegetation indexes and with PC and PH. This is evidence that BRT is 
dominated by the exposed soil sensed by Landsat. In fact, BRT is com- 
monly called the soil brightness index. The linear correlation (r= 
-0.739**) between PH and BRT in Table 7.1, part A, was the highest 
correlation coefficient observed in the data set, excluding the green- 
ness vegetation index intercorrelations. 

In summary, the various analyses presented show that temperature 
as observed by HCMM is negatively correlated with the greenness vegeta- 
tion indexes. Consequently, as the amount of living vegetation increased 
surface temperature decreased. On the other hand, the Kauth brightness 
(BRT) was negatively linearly correlated with the ground estimates of 
PC (r=-0.501*) and PH (r=~0.739**) , but BRT was not significantly cor- 
related with HCMM measured temperature nor with the greenness vegetation 
indexes. Therefore, the amount of green vegetation apparently had a 
stronger effect on surface temperature than did the amount of exposed 
(bare) soil. Perhaps this is so because vegetation actively affects 
several aspects of the energy balance--aerodynamic roughness, light 
energy partitioning per se , and transpiration, a dominant process among 
the energy dissipation mechanisms. 

The linear correlation coefficient matrix is given in Table 7.2 
for the 9 sorghum data points from the pooled crops data set of Table 
7.1 . Percent cover was significantly negatively correlated (r=-0.682**) 
with HCMM temperature and both PC and PH were significantly positively 
correlated with each of the vegetation indexes and significantly nega- 
tively correlated with BRT. The greenness vegetation indexes were 
negatively correlated (r=-0.706* to -0.719*) with HCMM temperatures, 
and brightness was positively correlated (r -0.835**) with HCMM tempera- 
ture. 

The relation between PC and the vegetation indexes PVI and GVI is 
displayed in the upper part of Figure 7.1 , and the relation between the 
vegetation indexes and HCMM indicated temperatures is shown in the 
lower part of the figure. The data were pooled for three dates and 
it is evident that the 03 Jul ground covers are low (due to harvest of 
the fields) and the temperatures are high. 


. . . . .. - - ^av'-2S&> • — '• . . - 
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7.4.2 Influence of other variables relative to that of 
vegetation on HCMM temperatures and albedos 


The data set of part C, Table 7.1, was segregated into irrigated 
and dryland test sites for the analyses of this section. The objective 

was to look at variables other than vegetation that affected tempera- 
ture and albedo. 

The results are summarized in Table 7.3 where the analysis of the 
irrigated sites in Part A of the table and the nonirrigated sites are 
Part B. For the irrigated sites, seven variables were highly signifi- 
cantly correlated with HCMM temperature for a 3x3 matrix of pixels; 
the variables, in order of decreasing correlation coefficients, were: 
insolation on the day of the pass (IDOP), maximum air temperature (Tmax), 
distance from the coast (DIST), saturation deficit of the air (SDEF), 
perpendicular vegetation index from Landsat bands 5 and 6 (PVI6), 

Kauth greenness (GVI) and days since rainfall or irrigation (DSR). 

For the dryland sites, five variables were highly significantly cor- 
related with the HCMM temperature. They were, in order of decreasing 
magnitude of the correlation coefficients, SDEF, Tmax, DIST, IDOP, and 
DSR. 

The step-wise multiple regression equations that included six 
variables are given at the bottom of each part of Table 7.3. For both 
the irrigated and dryland sites, six variables explain 90% of the 
variation in the HCMM observed temperatures. The standard error of 
the estimate of temperature is 1.4 C for the irrigated sites and 1,5 
C for the dryland sites. The variables included in the regression 
analyses are predominantly weather related (insolation, saturation 
deficit of the air, and maximum daily air temperature), geographic 
(distance from the coast) and related to evapotranspi ration (days 
since rainfall or irrigation). For the mixture of crops and planting 
configurations represented in the data set, the spectral vegetation 
index, PVI6, was the third variable chosen in the regression equation 
for the dryland sites but no vegetation index was among the first 5 
variables chosen for the irrigated sites, although PVI and GVI were 
the 7th and 8th variables chosen. 


From these analyses we conclude that when several crops are 
pooled, variables such as insolation which relates directly to energy 
available, air temperature which responds to the sensible heat com- 
ponent, saturation deficit of the air which is a measure of the 
capacity of air to take on moisture, and days since irrigation or 
rainfall which is indicating availability of water for evapotran- 
spiration relate more directly to the HCIIM temperatures than do in- 
dicators of vegetation greenness. When crops are pooled that differ, 
in stage of maturity (or senescence), planting configurations, and 
that vary crop to crop in relation between height and percent ground 
cover (or its inverse, exposed soil) then the vegetation indexes are 
not as closely correlated with observed HCMM temperatures as are the 
weather variables. These observations are borne out for the nine 
sorghum fields of Table 7.2 The multiple regression equation for 
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estimating T3 using 5 variables was T3=-5.70 - 1.075(SDEF) - 2.107 
(PVI) + 1.990(PVI6) + 0. 030( PC) - 0.005(PH) for which R z =0.998 and the 
standard error of the estimate of T3 was 0.13 C. Within this single 
crop, 4 of the 5 variables included were vegetation parameters and 
only one was weather related. The SDEF alone explained 76% of the 
variation in surface temperature of the grain sorghun fields, SDEF + 
pyi6 explained 91%, and SDEF + PVI 6 +• PVI explained 99.2% of the varia- 
tion. These findings suggest the desirability of looking at relation- 
ships within crops. We attempt that in the following sections. 

7.4.3 Percent cover versus plant height within crops 

The data for the analyses are the ground truth for the intensive 
sites corresponding to the eight HCMM overpass dates (27 May, 06, 12, 

17 and 23 Jun; 03 and 08 Jul ; and 15 Aug 78). Table 7.4 describes the 
data in terms of number of observations per crop, summarizes the means 
and standard deviations in both PC and PH, and presents the equations 
that relate PC to PH. The intercept values differed considerably among 
crops and so did the coefficients for the PH term. The growing season 
was Well along before plant measurements were begun in this study so 
that an intercept value greatly different from zero (the predicted PC 
at zero PH) is not alarming. The quadratic equation for sorghum re- 
sulted from head extension above the leaf canopy and the fact that 
some plant heights following combining and even tillage are included 
in the data set (Consequently, there are observations at low PC and 
PH.). Within row crops, plant population (plants per unit area) was 
variable; this introduced variability in PC without affecting PH. In 
plotting the data by soil type for cotton and sugarcane for the water 
use studies (Chapter 6.0), it was evident that on clay soils the PC 
was greater for a given plant height than for the nonclay soil types. 
The slopes were the same, but the intercepts differed. Thus soil 
type caused some of the variability in the relation between PC and PH. 

In summary, there was a highly Significant relation between PH 
and PC for sorghum, sugarcane, cotton, and citrus. However, no two 
crops gave similar percent cover prediction equations due to differing 
slopes and intercepts. 

7.4.4 Temperature versus percent cover within crops 

The data for this analysis are the ground truth for the intensive 
sites corresponding to the eight HCMM overpass dates from 27 May to 
15 Aug 78. Table 7.5 describes the data in terms of number of obser- 
vations per crop, summarizes the means and standard deviations in both 
PC and T3, and presents the equations that relate PC to T3. 

Grain sorghum was the only crop for which the relationship was 
significant (r z =0.522**). The regression equation predicted a surface 
temperature of 36.5 C for 0% ground cover and a decrease in temperature 
of 1.3 C for each 10% increase in ground cover. This result appears 
physical ly meaningful until one considers that sorghum was senesci ng 
during the period Jun to Aug so that the ground cover was decreasing 
at the same time air temperature was increasing. Sugarcane increased 
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in ground cover during the 27 May to 15 Aug time period, and it had 
no correlation between ground cover and crop canopy temperature. Thus 
we believe that for sugarcane the cooling effect of an increase in 
ground cover was offset by the seasonal increase in ambient air and 
agricultural landscape temperature during the observation period. The 
citrus sites all rad 54 ± 2 % ground cover so that changes in their tem- 
perature were a function of time and other variables besides ground 
cover. There were only 8 observations for cotton and they seemed to 
behave like sugarcane. Cotton had a lower temperature than the other 
crops because the temperature observations were limited to the first 
three dates 27 May, 06 Jun and 12 Jun 78. 

Within the range of ground covers of the four crops studied, ground 
cover was not very strongly related to observed surface temperature. 

For the crops that were increasing in ground cover during the observa- 
tion period (sugarcane, cotton), the cooling effect of increasing ground 
cover was about equally offset by the seasonal increase in ambient air 
and agricultural landscape temperature during the observation period,- 

7,4.5 Temperature comparison among crops 

One of the objectives of this study was to determine whether surface 
or canopy temperature differed among crops, HCMM resolution permitting. 
Although we had eight HCMM overpasses in a 2h month period from the end 
of May to mid-August, we deleted approximately one-third of the tempera- 
ture observations due to cloud or vapor contamination. Therefore, no 
ground site was consistently represented in the data. In addition, one 
date (17 Jun 78) yielded low temperatures - probably due to high cirrus 
cloud contamination - and was eliminated completely. However, by scan- 
ning the summaries of the temperature data by Individual ground site, we 
found three sorghum sites, two sugarcane sites and two citrus sites were 
represented on five of the eight dates. 


The comparison of temperature among crop surfaces as well as rela- 
tive to that of a deep water body and maximum air temperature at the 
central weather station are presented in Figure 7.2. The data show the 
average temperature for the vegetated sites is 4.4 C lower on 06 Jun 
than on 27 May. These observations were made preceding and just after 
heavy rains on 04 to 06 Jun, respectively. It was essentially rainless 
again for the remainder of the study period. 

The sorghum grain was maturing in June, the crop was harvested by 
mid-July, and the land was fallow by 15 Aug. By contrast, sugarcane 
(an irrigated crop) was growing rapidly from June to August. The 
perennial citrus trees, also irrigated, were essentially invariant 
in ground cover during this fraction of a growing season. 


The temperature of all the crops increased after the 06 Jun 
observation date. The temperature of sugarcane and citrus tends to 
be similar and both are cooler than the grain sorghum fields. On 15 
Aug the citrus was about 3 C cooler than the fallow sorghum, and the 
Lake water is about 10 C cooler than the cropped sites. For these 
uncorrected temperature data, maximum daily air temperature was above 
crop surface temperature, except on 15 Aug 73 when the fallow sorghum 
was above air temperature, ■ 
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Generally, the data show that the temperature spread among crops 

is narrow and that the uncorrected crop surface temperature was below 

daily maximum air temperature when the sites were vegetated. 
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Table 7.1 Linear correlation coefficient matrices among HCMM albedo (113), HCMM 

temperature (T3) , and plant parameters for mixed crops. 


A. Intensive sites (HCMM, Landsat and ground truth(GT) observations) 

252 PVT PVI 6 BRT GVI PC 


PVI6 BRT GVI PC 

- dimensionless (%) 



R 

PVI 

PVI 6 

BRT 

GVI 

PC 

PH 

T3 


Mean! 17.3 


n » 22 r (0.05) = 0.423 ; r (0.01) « 0.537 

Observations of each crop: sorghum, 9; sugarcane, 5; citrus, 4; 

cotton, 1 and grass, 3 

Representative and weather station sites (HCMM and Landsat observations) 
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• 500 ** 
.967 
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Mean j 17.2 


n = 30 r (0.05) = 0.361 ; r (0.01) - 0.463 

Observations of each crop: sorghum, 4; sugarcane, 1; citrus, 5; 
mixed crops , 1 2; range , 4 ; wildlife refuge , 2 and mud f la t , 2 






Above data sets combined 
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Mean 17.2 
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Table 7.2 Linear correlation coefficient matrix among HCMM albedo, HCMM temperature 
and plant parameters for grain sorghum, r 



R3 

PVI 

PVI 6 

BRT 

GVI 

PC 

PH 

T3 


«> 


—dimensionless— 

— 

(%) 

(cm) 

(°C) 

R 

- 

.265 

.286 

.197 

.282 

.349 

.387 

,122 

PVI 


- 

.996** 

-.651* 

.998* 

. 740** 

.670* 

-.719* 

PVI6 




632* 

-.998** 

. 750** 

. 666* 

-.718* 

BRT 




- 

-.622* 

-.733* 

-.704* 

' 

.835** 

GVI 


r 



- 

.737** 

.653* 

- . 706* 

PC 






- 

.814** 

-.682* 

PH 

! 






- 

-.501 

T3 

• . i 



L:- 




- 

Mean ■ 


1 . I 

16-.5 

81.7 

19.9 

31,0 

64.2 

33.2 

S.D. 


warn 

10. i 

9.0 

11.0 

22.2 

35.6 

1 

3.. 7 







Table 


7.3 Linear correlation coefficient matrix for irrigated sites (part A) and dryland sites (part B) 
Data are from HCMM, Landsat, and weather observations for JD 163, 168, and 189. 


Irriga ted si te s (n - 2 6) 
IDOP T' SDEF 


|R3 

i;. 

IDOP 

SDEF 

DIST 

PVI 

PVI6 

BRT 

CVI 

T.. 


0.276 -0.218 


0.702 
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(km) 


Pyi6 BRT 

-dimensionless- — 


0.610 0.079 ! 0.140 


0.161 


0.580 -0.342 j -0.511 -0.166 

• ' i ' . * ' 


0.138 -0.163 


30.6 80.7 


-0.409 


0.055 


-0.356 1 -0.374 j -0.104 

i i 

I JLX J 

f-0.963 ! 0.274 

j. ! 

i - ! 0.267 


Max 

(®C) (days) 



it* 


19.9 i 76.4 23.5 


0.091 

-0.058 

0.301 

** 

0.876 

k k 

0.602 

** 

0.863 

it* 

0.841 

■A* 

0.580 

** 

0.664 

0.387 

0.182 

** 

0.699 

-0.183 

-0.417* 

* 

-0.454 

-0.394* 

! ** 

[-0.548 

! 

** 

-0.592 

-0.023 

-0.105 

-0.131 i 

-0.340 

-0,517 ' 

** 

-0.555 

j 

i 

0.430" 

A. A. 

0.722" ‘ 

•kiz 

0.541 

36.1 

17.1 

31.1 

1.5 

10.5 

4.0 
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Multiple regression equation for 6 variables: 


X3 = -5.200 - 0. 248(SDEF) + 0.044 (DIST) - 0.244(PVI6) + 1.684(Tmax) + 1.396(DSR) - 10.423((DSR) ) 

2 

R = 0.897;; standard error of the estimate of temperature = 1.5 C. 




ID OP 

SDEF ' 

DIST 
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- ' 
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T3 
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76.6 
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54.3 
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finai] 


p 


0.411* 0.391* 
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0.579** 0.319 
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0.977** -0.028 

0.000 


0.470* 0.141 
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0.038 0.237 
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-0.591** -0.378 
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Table 7.4 Relation between percent cover (PC) and plant height (PH) for 

four different crops. The data is from ground truth observations 

corresponding to 8 HCMH overpass dates between 27 May and 15 Aug. 
78, for the 22 intensive sites. 


Attribute Crop 


No. of observations: 

Sorghum 

47 

Sugarcane 

35 

Cotton 

28 

Citrus 

22 

PH Mean (cm): 

60.1 

169.9 

57.6 

. •_ 

367.5 

PH Stnd. dev. (cm): 

33.8 

23.3 

10.2 

57.7 

PC Mean (%) : 

30.2 

. ' ■ ' 1 

54.4 ; 

52.0 

53.4 

1 

PC Stnd. dev, (%) : 

..... ■ , 

21.9 

16.1 

17.3 

2.4 


Crop 

- ■ ■■ ^ 

Equation 

2 o 2 

r or R 

• ■ ■ 

Stnd. error 
of PC estimate 
(%) 

Sorghum 

PC - 4.375 - 0.302 (PH) + 0.009(PH) 2 

** 

0.872' 

8.1 

Sugarcane 

PC = -50.3 + 0.891(PH) 

A* 

0.794 

I 

7.4 

Cotton 

PC = -33.2 + 1 .47 8 (PH) 

1 ** 
0.764 

8.6 

Citrus 

1 

PC = 68.6 - 0.041 (PH) 

** 

0.989 

0.04 
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Table 7.5 Relation between plant cover and temperature among crops. The data 
are from ground observations corresponding to 8 HCMM overpass dates 

between 27 May an.d 15 Aug. 78, for the 22 intensive test sites. 


Attribute 



Sorghum 

Sugarcane 

Cotton 

Citrus 

No. of observations: 

33 

18 

8 

16 

PC Mean (%): 

28.0 

56.7 

41.0 

54.1 

PC stnd, dev* (%): 

20.2 

15.4 

14.7 

2.0 

T3 Mean(°C) 

. 

32.7 

30.6 

29.7 

32.2 

T3 stnd. dev. 

3.8 

2.8 

2.3 

2.4 


Stnd. error 


Equation 


Sorghum 1 

T3= 36.51-0. 135(PC) 

.522** 

2.6 

Sugarcane 

T3= 29.85+. 012 (PC) 

.005 

2.9 

Cotton 

T3= 28. 17+. 037 (PC) 

.058 

2.3 

Citrus 

T3= 3.58+. 530 (PC) 

.192 

2.3 
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Table A. 1 


Composition of intensive sites by fields, field sizes, and 
crops within larger sample segments. 


Site No. 

Typifying 

field 

Hectares 

Crop 2 ' 

Sample 

segment 

Hectares 

No. of 
fields 

10101 

T 

76 

SR 

1010A 

392 

6 

10102 

2 

99 

CO 




30901 

1 

1538 

SR 

3090A 

1538 

1 

31001 

1 

1237 

SR 

3100 

1237 

1 

31011 

1 

145 

SR 

3101 

145 

1 

41341 

1 

379 

SR 

4134 

379 

1 

41381 

1 

312 

SR 

4138 

482 

3 

41421 

2 

19 

SR 

4142 

123 

8 

41422 

2 

33 

SC 



6 

10194 

4 

118 

SC 

1019 

395 

7 

10196 

6 

73 

SC 




30951 

1 

79 

sc 

3095 

187 

6 

30955 

5 

33 

CO 




41665 

5 

125 

sc 

£166 

202 

6 

51782 

2 

50 

sc 

5178 

116 

3 

10381 

1,2, 6, 9 

35,40,17,26 

CO 

1038 

290 

10 

20611 

1 

231 

CO 

2061 

252 

2 

20681 

11,8 

101,15 

CT 

2068A 

282 

16 

20754 

4 

82 

CT 

2075 

280 

6 

31166 

6, 4, 5,7 

18,29,9,4 

CT 

3116 

108 

10 

20612 

1 

272 

GR^ 

2061A 

272 

1 

20651 

1 

480 

GR^ 

2965A 

480 

1 


-^R, sorghum; CO, cotton; SC, sugarcane; CT, citrus; GR, grass 
^20612, Bermudagrass; 20651, buf felgrass 


Table A. 2. 

Ground truth 

summary for 

intensive sites for 

three sampling 

dates. 


Julian 

PH 

PC 

DSR 

TMAX 

TMIN 

DEPL 

ET/ETp 

DIST 

WCAP 

Site No. 

Date 

(cm) 

(%) 

(days) 

(°c) 

(°C) 

(mm H 2 0) (-) 

(km) 

(cm) 

SORGHUM 











10101 

163 

82 

40 

5 

35.0 

22.2 

25.4 

.58 

50.61 

17.5 


168 

83 

40 

10 

33.9 

21.1 

41.6 

.46 




189 

40 

2 

31 

37.8 

21.7 

89.2 

.17 



30901 

163 

70 

28 

3 

34.4 

23.9 

12.7 

.52 

85.61 

15.5 


168 

70 

28 

8 

35.6 

22.8 

27.1 

.40 




189 

0(BS) 

5 

29 

37.2 

24.4 

65.0 

.13 



31001 

163 

95 

56 

3 

34.4 

23.9 

14.0 

.59 

99.14 

15.5 


168 

95 

56 

8 

35.6 

22.8 

30.8 

.46 




189 

10 (BS) 

5 

29 

37.2 

24.4 

75.4 

.16 



31011 

163 

90 

56 

3 

34.4 

23.9 

14.0 

.59 

99.59 

15.5 


168 

90 -nc 

. 57 

8 

35.6 

22.8 

30.8 

.46 




189 

15 (weeds) 10 

29 

37.2 

24.4 

75.4 

.16 



41341 

163 

60 

21 

8 

36.1 

25.0 

25.9 

.39 

116.51 

15.5 


168 

61 

20 

13 

36.1 

23.3 

37 .9 

.32 




189 

62 

18 

34 

37.8 

24.4 

64.8 

.08 



41381 

163 

93 

65 

5 

35.0 

23.9 

28.7 

.74 

48.4 

17.5 


168 

93 

63 

10 

35.6 

22.2 

51.9 

.65 




189 

0 (BS ) 

0 

31 

37.8 

23.9 

114.2 

.24 



41421 

163 

84 

19 

5 

34.4 

23.9 

22.4 

.52 

88.9 

15.5 


168 

84 

18 

10 

35.6 

22.8 

36.4 

.58 




189 

88 

20 

31 

37.2 

24.4 

77.7 

.14 



SUGARCANE 











10194 

163 

ISO 

71 

4 

35.0 

22.2 

26.7 

.94 

55.0 

17.5 


168 

185 

75 

9 

33.9 

21.1 

57.2 

.88 




189 

220 

83 

17 

37.8 

21.7 

114.2 

.82 



10196 

163 

150 

45 

4 

35.0 

22.2 

26.5 

.92 

54.1 

17.5 


168 

158 

49 

9 

33.9 

21.1 

56.8 

.88 




189 

203 

72 

19 

37.8 

21.7 

126,1 

,82 



3093 1 

163 

150 

95 

5 

35.0 

23.9 

34.7 

.96 

51,9 

17.5 


168 

160 

95 

10 

35.6 

22.2 

6 6. 5 

.92 




189 

205 

94 

10 

37.8 

23.9 

72.7 

.91 



41422 

163 

151 

48 

0 

34.4 

23.9 

0.0 

.89 

89.4 

15.5 


168 

158 

51 

5 

35.6 

22.8 

32.5 

.95 




189 

204 

75 

11 

37.2 

24.4 

78.2 

.87 



41665 

163 

155 

33 

4 

35.0 

22.2 

24.2 

.84 

66.8 

19.5 


168 

155 

35 

9 

r rvey 

33.9 

21.1 

51.1 

.76 




189 

160 70 (weeds) 6 

37.8 

21.7 

45.5 

, 95 
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Table A. 2 Continued 



Julian 

PH 

PC 

DSR TMAX 

TMIN 

DEPL 

ET/ETp 

DIST 

WCAP 

Site No. 

Date 

(c~) 

(%) 

(days) (*C) 

(°C) 

(ram ^0) 

(-) 

(km) 

(cm) 


51782 

163 

155 

47 

9 

35.0 

22.2 

53.7 

.89 

91.3 

19.5 

168 

170 

49 

14 

33.9 

21.1 

82.2 

.82 




189 

215 

72 

18 

37.8 

21.7 

90.0 

.80 ' 



COTTON 

10102 

163 

45 

26 

5 

35.0 

22.2 

30.2 

.80 

61.0 

17.5 


168 

45 

28 

10 

33.9 

21.1 

56.1 

.75 




189 

49 

38 

31 

37.8 

21.7 

145.0 

.46 



10381 

163 

52 

62 

9 

35.0 

22.2 

51.2 

.95 

69.9 

19.5 


168 

58 

70 

14 

33.9 

21.1 

81.8 

.88 




189 

66 

81 

1 

37.8 

21 . 7 

6.8 

.86 



20611 

163 

70 

64 

5 

35.0 

22.2 

30.9 

.84 

63.2 

15.5 


168 

70 

68 

10 

33.9 

21.1 

57.3 

.75 




189 

70 

66 

31 

37.8 

21.7 

111.2 

.00 



30955 

163 

60 

50 

5 

35.0 

23.9 

33.1 

.91 

51.8 

15.5 


168 

62 

54 

10 

35.6 

22.2 

62.6 

,84 




189 

70 

62 

31 

37.8 

23.9 

133.5 

.31 



CITRUS 

20681 

163 

450 

50 

5 

35.0 

22.2 

19.5 

.41 

81.8 

17.5 


168 

450 

50 

2 

33.9 

21.1 

8.9 

.75 




189 

450 

50 

18 

37.8 

21.7 

77.8 

.41 



20754 

163 

330 

55 

19 

36.1 

25.0 

55 . 0 

.43 

110.7 

15.5 


168 

330 

55 

24 

36.1 

23.3 

69.0 

.42 




189 

330 

55 

14 

37.8 

24.4 

54.1 

.45 



31166 

163 

328 

55 

9 

36.1 

25.0 

31.8 

.40 

97.2 

17.5 


168 

328 

55 

14 

36.1 

23.3 

44.8 

.38 




189 

328 

55 

17 

37.8 

24.4 

64.2 

.43 



BERMUDAGRASS 










20612 

163 

11 

60 

5 

35.0 

22 . 2 

27.2 

.74 

61.8 

15.5 


168 

11 

60 

10 

33.9 

21 . 1 

51.0 

.68 




189 

12 

70 

31 

37.8 

21.7 

119.6 

,17 



BUFFEIGRASS 

29651 

163 

22 

75 

3 

34.4 

23.9 

14.7 

.67 

109.7 

15.5 


168 

22 

75 

8 

35.6 

22.8 

34.1 

.52 




189 

20 

70 

29 

37.2 

24.4 

48.4 

. 00 
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Table A. 3 Landsat digital count sample means and vegetation indexes for 
two dates. 



Landsac 

Digital 

Counts 


Vegetation 

Indexes 


Site r 

MSS 4 

MSS 5 

MSS6 

MSS7 ‘ r 

PVI 

PVI6 

BRT 

GVI 1 

A, June 

15, 1978 








41381 

25 

54 

55 

27 

15.4 

21.5 

65.9 

25.1 

30901 

27 

29 

65 

32 

18.2 

26.2 

77.7 

30.5 

31001 

27 

28 

66 

33 

19.8 

27.7 

77.8 

32.5 

31011 

31 

35 

66 

30 

14.3 

22.2 

82.8 

25 .4 

41341 

31 

39 

58 

26 

9.6 

14.5 

79.5 

17.2 

41422 

33 

34 

85 

40 

23.9 

36.9 

97.3 

41.5 

10196 

24 

21 

73 

38 

27.4 

37.6 

77.3 

43.9 

41665 

25 

24 

61 

29 

17.5 

26.4 

69.8 

29.9 

51782 

28 

27 

60 

28 

16.0 

24.0 

72.5 

27 .0 

10102 

33 

37 

66 

28 

12.0 

21.4 

84.2 

23.0 

10381 

26 

25 

76 

35 

23.0 

37.0 

82.0 

41.5 

20681 

23 

25 

48 

22 

10.9 

16.8 

60.6 

19.3 j\ 

31166 

21 

21 

52 

24 

14.7 

21.7 

60.3 

25.1 

20754 

25 

28 

49 

22 

9.9 

14.8 

63.6 

17.3 

10194 

25 

24 

61 

29 

17.5 

26.4 

69.8 

30.0 

2061 2 

29 

35 

51 

22 

7.5 

12.0 

71.6 

13.8 

20651 

33 

43 

56 

24 

5.6 

10.5 

81.6 

11.9 


Table A. 3 Continued 


Landsat Digital Counts Vegetation Indexes 


Site 

MSS 4 

ks 35 

MSS 6 

MSS7 

PVI 

PVI 6 

BRT 

GVI 

B. July 

12, 1978 







41381 

32 

40 

42 

19 

1.9 

2.0 

69.8 

2.4 

30901 

35 

50 

54 

24 

3.2 

3.6 

85.5 

5.8 

31001 

39 

58 

64 

29 

4.6 

5.7 

99.2 

8.5 

31011 

38 

56 

66 

31 

7.2 

8.6 

99.7 

12.4 

41341 

33 

47 

60 

28 

8.0 

10.1 

87.1 

13.6 

41421 

34 

46 

52 

24 

4.8 

5.2 

81.3 

7.5 

41422 

29 

28 

68 

38 

24.2 

29.0 

80.9 

35.2 

10194 

35 

37 

72 

37 

20.0 

25.5 

91.6 

30.1 

10196 

28 

24 

70 

40 

27.4 

33.0 

80.1 

39.7 

41665 

29 

31 

52 

26 

12.0 

15.1 

70.4 

18.1 

51782 

29. 

27 

66 

36 

22.9 ' 

28.0 

79.1 

33.6 

30955 

30 

36 

57 

28 

12.4 

15.7 

77.0 

19.3 

10102 

32 

38 

52 

24 

8.0 

10.4 

75.9 

12.6 

20611 

28 

33 

51 

26 

11.4 

13.3 

70.7 

17.0 

10381 

24 

22 

69 

37 

25.7 

34.2 

75.4 

40.4 

20681 

28 

31 

55 

28 

14.4 

17.5 

72.3 

21.7 

31166 

25 

26 

52 

27 

15.3 

18.7 

66.1 

23.0 

20754 

27 

31 

49 

25 

11.5 

13.2 

67.4 

16.7 

20651 

33 

48 

52 

23 

3.4 

3.8 

82.4 

6.1 

20612 

30 

35 

45 

22 

6.7 

7.2 

67.8 

9.4 


Table A. 4 Temperatures (Tl, T3, T9) and albedos (Rl, R3, R9) and their 
standard deviations for 1x1, 3x3, and 9x9 pixel samples by 

site number for 12 June 1978 (J*D. 163). 


Site No. 

Tl 

T3 

T9 

R1 

R3 

R9 

10101 

***** 

***** 

***** 

***** 

***** 

***** 


*** 

*** 

*** 

*** 

*** 

*** 

30901 

31.71 

31.82 

31.39 

18.03 

17.55 

18.37 



0.4 

1.6 


0.8 

1.8 

31001 

32.34 

32.45 

32.55 

19.21 

19.25 

19.34 



0.4 

0.7 


0.6 

0.7 

31011 

33.90 

33.03 

31.98 

18.43 

19.08 

19.08 



0.9 

0.9 


0.4 

0.8 

41341 

19.84 

23.83 

31.16 

20.39 

17.73 

18.27 



4.2 

4.9 


2.9 

2.8 

41381 

(20.88) 

(22.77) 

(26.31) 

(20.39) 

(20.94) 

(21.31) 



(2.3) 

(2.7) 


(2.1) 

(1.6) 

41421 

30.45 

30.24 

30.40 

18.43 

18.43 

19.30 



0.6 

0.9 


1.3 

1.6 

10194 

(26.58) 

(25.28) 

(25.18) 

(18.82) 

(19.82) 

(20.27) 



(2.2) 

(1.6) 


(1.5) 

(1.4) 

10196 

(25.26) 

(23.72) 

(25.18) 

(20.00) 

(21.17) 

(20.62) 



(2.3) 

(1.7) 


(1.4) 

(1.4) 

30951 

(26.91) 

(27.51) 

(27.28) 

(19.60) 

(20.61) 

(20.61) 



(1.0) 

(2.0) 


(1.4) 

(2 . 1) 

41422 

30.77 

30.55 

30.52 

19.21 

18.73 

19.36 



0.7 

0.9 


1.0 

1.5 

41665 

(28.21) 

(27.14) 

(25.81) 

(17.25) 

(18.87) 

(19.08) 



(0.8) 

(3.0) 


(1.0) 

(2.1) 

51782 

(30.13) 

(29.45) 

(27.70) 

(21.96) 

(22.10) 

(20.21) 



(1.0) 

(2.7) 


(1.3) 

( 2.3) 

10102 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 

*** 

*** 

*** 

10381 

(22.58) 

(23.38) 

(25.93) 

(13.33) 

(16.66) 

(18.68) 



(3.8) 

(3.1) 


(1.8) 

(2.2) 

20611 

29.49 

29.21 

29.34 

24.31 

20.39 

19.11 



1.2 

2.7 


2.4 

2.3 

30955 

(20.88) 

(23.57) 

(26.56) 

(18.43) 

(19.60) 

(21 . 21) 



(2.7) 

(2.3) 


(1.3) 

(2.0) 

20681 

29.17 

29.60 

28.75 

19.21 

18.56 

18.74 



0.6 

2.3 


0.4 

1.3 

20754 

33.28 

33.62 

33.47 

19.60 

18.64 

18.49 



0.8 

0.9 


0.7 

0.7 

31166 

31.71 

33.17 

33.59 

17.64 

18.16 

18.16 



1.0 

1.0 


0.5 

0.8 

20612 

24.26 

25.67 

28.36 

18.03 

19.34 

19.60 



3.2 

2.7 


2.9 

2.3 

20651 

35.44 

34.89 

34. 10 

20.00 

19.08 

19.22 



0.6 

1.1 


0.7 

1.0 


^ Asterisks denote deleted cloud contaminated data: Parentheses denote sub-visible 
cirrus contaminated data (also deleted). 
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Table A. 4 continued. Data for 17 June 1978 (J.D. 168) 


Site No. 

T1 

T3 

T9 

Rl 

i : 

R3 

R9 

10101 

23.93 

24.30 

24.78 

16.86 

15.77 

16.15 



0.4 

0.5 


1.0 

1.5 

30901 

29.17 

29.03 

28.82 

16.47 

16.07 

17.24 



0.5 

1.2 


1.5 

1.7 

31001 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 

31011 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 

41341 

31.08 

30.58 

32.89 

18.82 

17.74 

16.65 



2.8 

2.4 


3.0 

2.8 

41381 

27.23 

26.69 

26.08 

16.07 

15.86 

16.48 



0.6 

0.8 


0.6 

1.1 

41421 

*;* tick* 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 

10194 

25.26 

24.71 

24.79 

14.50 

15.64 

16.15 



0.4 

0.5 


1.2 

1.1 

10196 

24.26 

24.49 

24.67 

14.90 

15.68 

16.33 



0.3 

0.5 


1.1 

1.4 

30951 

24.26 

24.93 

25.25 

16.86 

17.08 

16.00 



0.4 

0.8 


0.9 

1.1 

41422 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


‘ *** 

*** 

41665 

(23.60) 

(.25.00) 

(24.97) 

(14.50) 

(16.07) 

(17.52) 



(0.9) 

(1.8) 


(2.1) 

(2.2) 

51782 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 

10102 

24.26 

24.49 

24.79 

15.68 

15.55 

15.98 



0.5 

0.4 


1.4 

1.5 

10381 

(26.25) 

(25.70) 

(26.55) 

(18.43) 

(17.47) 

(17.72) 



(0.5) 

(2.0) 


(0.9) 

(1.3) 

20611 

26.25 

27.05 

26.34 

16.86 

16.47 

15.78 



0.8 

1.1 


1.3 

l.i 

30955 

24.94 

24.97 

25.20 

15.68 

16.12 

16.19 



0.3 

0.6 


0.7 

1.1 

20681 

***** 

***** 

***** 

***** 

***** 

***** 



*** - ; 

*** 


*** 

*** 

20754 

(34.21) 

(32.23) 

(29.79) 

(16.47) 

(17.29) 

(17.41) 



(1.2) 

(2.3) 


(1.3) 

(3.5) 

31166 

***** 

***** 

***** 

***** 

***** 

***** 




*** 


*** 

*** 

20612 

25-26 

25.59 

25.75 

17.25 

16.64 ' 

15.84 



0.5 

0. 9 


1.0 

1.3 

20651 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 
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Table A. 4 
Site No. 

continued . 
Tl 

Data are 
T3 

for 08 July 
T9 

1978 (J.D. 
R1 

189). 

R3 

R9 

10101 

(23.26) 

(23.70) 

(25.02) 

(17.25) 

(16.94) 

(17.06) 



(1.1) 

(3.0) 


(0.5) 

(0.6) 

30901 

34.82 

35.30 

35.89 

14.50 

15.20 

15.91 



1.1 

1.0 


0.5 ! 

1.1 

31001 

36.35 

36.62 

36.30 

18.43 

18.16 

18.06 



0.5 

0.7 

0.3 

0.3 

0.8 

31011 

35.44 

35.54 

35.95 

15.29 

16.60 

17.85 



1.0 

1.1 


1.1 

1.2 

41341 

38.76 

38.49 

36.88 

18.03 

16.03 

15.47 



0.7 

1.2 


1.5 

1.1 

41381 

(27.23) 

(27.05) 

(26.97) 

(16.07) 

(17.08) 

(17.08) 



(0.8) 

(1.3) 


(0.6) I 

(0.9) 

41421 

36.05 

35.22 

33.55 

18.43 

17.90 ' 

16.88 



1 . 7 

2.0 


0.7 

1.3 

10194 

(26.91) 

(25.13) 

(25.56) 

(16.47) 

(16.81) 

(17.02) 



(2.7) 

(2.7) 


(0.8) 

(0.8) 

10196 

25.59 

26.35. 

26.24 

16.07 

16.86 

16.83 



2.0 

2.8 


0.9 

0.9 

30951 

***** 

***** 

***** 

***** 

***** 

***** 



*** 

*** 


*** 

*** 

41422 

36.96 

36 . 01 

33.50 

18.43 

18.30 

16.82 



1.1 

2.0 


0.5 

1.3 

41665 

(31.40) 

(31.11) 

(30.33). 

(15.68) 

(15.94) 

(16.21) 



(1-5) 

(1.7) 


(0.8) 

(1.1) 

51782 

29. SI 

31.28 

32.15 

18.43 

17.82 

17.85 



1.9 

2.1 


0.9 

0.9 

10102 

(22.92) 

(23.37) 

(25.75) 

(18.03) 

(17.08) 

(17.18) 



(1.1) 

(3.5) 

i 

(0.6) 

(0.9) 

10381 

(30.45) 

(30.41) 

' (30.10) 

! (18.43) 

(17.64) 

(16.99) 



(1.6) 

(1.9) 

i . 

1. ■ • 

(0.9) 

(1.3) 

20611 

(28.53) 

(28.43) 

(28.77) 

(16.47) 

(16.20) 

(15.62) 



(2.3) 

(2.4) 


(0.3) 

(0.9) 

30955 

(27.23) 

(27.05) 

• (26.99) 

(18.43) 

(17.86) 

(17.25) 



(1.2) 

(2.1) 


(0.6) 

(0.9) 

20681 

(31.40) 

(31.78) 

(31.61) 

(16.47) 

(16.38) 

(16.31) 



(0.5) 

(1.2) 


(0.7) 

(1.2) 

20754 

36.96 ... 

36.15 

35.80 

16.07 

16.60 

16.64 



0.8 

1 .0 


0.9 

1.5 

31166 

36.05 

35.53 

36. 

16.47 

16.90 

16.89 



1.1 

0 9 


2.2 

1.4 

20612 

(26.25) 

(27.29) 

(2- 

(16. 47) 

(16.51) 

(15.89) 



(2.0) 

u 


(0.6) 

(1.1) 

20651 

36.05 

36.08 

35. b5 

17.64 

17.73 

16.98 



0.9 

1.7 


0.7 

1.3 
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Table A. 5 Rainfall summary for 22 intensive sites, mid-May to Mid-July 1978.-^ 


Rainfall Rainfall 


Site 

Date 

(mn) 

Site 

Date 

(mm) 

10101 

5/22 

4 

30901 

5/22 

16 

10102 

6/1 

71 

31001 

6/2 

4 


6/3 

54 

31011 

6/3 

29 


6/7 

32 

41341 

6/4 

11 


6/19 

5 


6/8 

25 


7/11 

11 




Source : 

Site, Monte 

Alto, Raymondville 

Source: 

McCook 


41381 

5/22 

4 

41421 

5/21 

14 

30951 

6/3 ; : V 

54 

41422 

6/1 

3 

30955 

6/7 

32 

20681 

6/3 

20 


6/19 ; 

5 


6/7 

23 


7/11 

11 




Source: 

Monte Alto, 

Raymondville 

Source: 

Edinburg 


20611 

5/21 

5 

10194 

6/2 

57 

20612 

6/3 

36 

10196 

6/3 

29 


6/7 

42 


6/4 

37 . 

Source: 

Monte Alto 



6/7 

5 





6/8 

17 

41665 

5/22 

3 


6/19 

7 


6/4 

21 


6/29 

4 


6/27 

4 

Source: 

Santa Rosa 


Source: 

Mercedes (6mi. SSE) 




51782 

5/22 

7 

20754 

5/22 

18 


6/2 

8 


6/4 

15 


6/3 

20 


6/8 

3 

Source: 

Schuster Farm 


7/24 

3 




Source: 

Bates Power 


10381 

6/3 

26 

31166 

6/1-3 

26 


6/7 

8 


7/23 

7 


6/29 

6 

Source: 

McAllen 



20651 5/22 18 

6/2-4 29 

6/8 13 

Source: Av. Bates Power and McCook 


^ Amounts less than 2mm not listed. 
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